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Phases of matter are sharply defined in the thermodynamic limit. One major challenge of accu¬ 
rately simulating quantum phase diagrams of interacting quantum systems is due to the fact that 
numerical simulations usually deal with the energy density, a local property of quantum wavefunc- 
tions, while identifying different quantum phases generally relies on long-range physics. In this 
paper we construct generic fully symmetric quantum wavefunctions under certain assumptions us¬ 
ing a type of tensor networks: projected entangled pair states, and provide practical simulation 
algorithms based on them. We find that quantum phases can be organized into crude classes distin¬ 
guished by short-range physics, which is related to the fractionalization of both on-site symmetries 
and space-group symmetries. Consequently, our simulation algorithms, which are useful to study 
long-range physics as well, are expected to be able to sharply determine crude classes in interact¬ 
ing quantum systems efficiently. Examples of these crude classes are demonstrated in half-integer 
quantum spin systems on the kagome lattice. Limitations and generalizations of our results are 
discussed. 
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I. INTRODUCTION 

Reliably simulating quantum phase diagrams of realis¬ 
tic interacting systems has been one of the central issues 
in condensed matter physics. A number of numerical 
methods have been developed in the past decades, in¬ 
cluding exact diagonalization, quantum Monte Carlo(for 
a review, see RefH]), variational Monte CarlcPH, the den¬ 
sity matrix renormalization group method (DMRGjP^, 
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and methods based on tensor network representations of 
quantum wavefunctionP^^^. Although with advantages 
and disadvantages, these methods have been demon¬ 
strated to be able to successfully simulate various inter¬ 
acting quantum models. For instance, an exotic quantum 
spin liquid phase have been recently identihed in the spin- 
1/2 Heisenberg model on the kagome latticdi^HlSI using 
DMRG methods. 

One major source of the challenges of accurately simu¬ 
lating realistic quantum models is the following fact. The 
full many-body quantum Hamiltonian cannot be exactly 
diagonalized as long as the sample size is not very small. 
Therefore even for intermediate sample size, except for 
systems that do not suffer from the sign problem in quan¬ 
tum Monte Carlo, one has to come up with variational 
wavefunctions, using which to search for the true ground 
states of quantum systems. The guiding principle of all 
variational simulations is simply to minimize the energy 
density, a local property of quantum states, of a given 
sample. On the other hand, generally distinguishing dif¬ 
ferent quantum phases relies on the long range physics. 
Consequently in these variational methods we are trying 
to determine long range physics based on local physics. 
However, competing quantum phases could have similar 
energy densities. In fact, it can be shown that differ¬ 
ent quantum phases could give arbitrarily close energy 

densitiePESl. 

To concretely demonstrate this challenge let’s consider 
frustrated quantum spin systems, for instance, nearest 
neighbor spin-1/2 Heisenberg models on the triangular 
lattice and the kagome lattice. In the triangular lattice 
case a good understanding of the ground state is known, 
based on results from various numerical simulation^i^lb^ 
which show that the system has a long range 120° copla- 
nar magnetic order. To establish this long range mag¬ 
netic order, a statement about the long-range physics, 
it is important to perform finite-size scaling since most 
numerical simulations study samples with small to inter¬ 
mediate sizes. The successful identification of the long 
range order in the triangular lattice model, to a large 
extent, is a consequence of the fact that the magnetic 
ordering in this system is quite stronji^ Even the finite 
size scalings performed quite some time agcP^ on small 
to intermediate sized samples give clear evidences of the 
order. 

The situation for the kagome lattice model is drasti¬ 
cally different. In the past it was known that even if a 
long range order does exist in this system, it is very weak. 
Thus in order to identify the presence or absence of a 
long range order, which represent two different quantum 
phases: a symmetry-breaking phase and a quantum spin 
liquid phase, one needs to perform finite size scaling in 
samples with larger sizes. The simulations on these sam¬ 
ples become possible only recently due to the progresses 
in DMRC methods. 

Cenerally speaking, in order to fully determine the 
quantum phase diagrams of correlated systems in numer¬ 
ical simulations, one cannot avoid studying samples with 


large sizes, simply because general quantum phases are 
sharply defined by the long range physics. But practi¬ 
cally the larger the system size is, the more challenging 
the simulation is. 

But are all quantum phases only distinguished by long- 
range physics? Before we provide an answer to this ques¬ 
tion, it is better to elaborate the question in a slightly 
sharper way. First we emphasize that the a phase is de¬ 
fined only when the global symmetry of the system is 
specified, which may or may not be spontaneously bro¬ 
ken. When limited to finite size samples, the ground state 
wavefunctions necessarily form (generally, irreducible) 
representations of the global symmetry which is usually a 
combination of on-site symmetries like spin rotations and 
space-group symmetries like translations. This statement 
is true even when the global symmetry is spontaneously 
broken in the long range physics. 

We again demonstrate the above statement in the con¬ 
text of frustrated spin-1/2 models. In this context, quan¬ 
tum spin liquids (QSL) are states of matter that do not 
break translation and spin rotation symmetries. In par¬ 
ticular, evidences of a fully gapped Z 2 QSL were reported 
in the kagome lattice model mentioned above. Recent 
theoretical worlPSl supports that this Z 2 QSL has a topo¬ 
logical order which can be described as a usual Z 2 gauge 
theory (i.e., the same topological order as in Kitaev’s 
toric code modeP^^). In such a Z 2 QSL, quasiparticle 
excitations include bosonic spin-1/2 spinon-e, bosonic 
vison-m, and their fermionic bound state / = em. Sup¬ 
pose that we can tune certain parameters in the spin 
model, it is possible that either the spinon e boson con¬ 
denses, which gives rise to certain long-range magneti¬ 
cally ordered (MO) phase, or the vison m boson con¬ 
denses, which gives rise to certain valence bond solid 
(VBS) phase since visons transform nontrivially under 
lattice symmetries. 

However, in this context, the boson condensations of e 
or m quasiparticles are only sharply defined in the long 
range physics. For instance, imagine one does a numeri¬ 
cal simulation for a phase transition between the Z 2 QSL 
and a nearby MO phase (VBS phase) via e (m) conden¬ 
sation. To avoid possible subtlety due to open boundary 
conditions, let’s consider a finite size torus sample. The 
ground state wavefunctions on both sides of the phase 
transition must share the same quantum number in the 
vicinity of the phase transition. Basically the quantum 
phase transition in the long-range physics is not visible 
on a hnite size sample unless a careful finite size scaling is 
performed in large system sizes.^For this reason, we say 
that the Z 2 QSL and the nearby MO phase (VBS phase) 
share the same short-range physics but are distinguished 
in the long-range physics. 

On the other hand, there are lots of examples in which 
ground states of different candidate quantum phases give 
distinct symmetry representations on sequences of finite 
size samples, which persist all the way to the thermo¬ 
dynamic limilP2H2l] Trivial examples include ferromag¬ 
netic phases and paramagnetic phases in spin systems. 
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As a somewhat nontrivial example, in a recent investi¬ 
gation of correlated electronic models on the honeycomb 
lattice, two candidate quantum phases: the chiral spin 
density wave phase and the d-fid superconductor phase, 
are found to host distinct lattice quantum numbers on 
47V X 4A^ X 2 symmetric samples^. In these cases, at 
least on these sequences of samples, clearly these candi¬ 
date phases really give completely different ground state 
wavefunctions which cannot be smoothly tuned from one 
to another. These quantum phases must be distinguished 
by short range physics. Note that the energy density, 
minimizing which is the guideline of all variational meth¬ 
ods, is also a short-range property of the wavefunction. 
One should have the hope of generically being able to 
sharply identify candidate phases distinguishable by short 
range physics even on small or intermediate sized sam¬ 
ples, without worrying about finite size scalings in larger 
system sizes. 

Note that we have made statements on “short-range 
physics” and “long-range physics” without sharply defin¬ 
ing their meanings. Now it is a good moment to comment 
on the sharp meanings of these terms used in this paper. 
By “long-range physics”, we really mean the long-range 
behavior of correlators measured in ground state wave- 
functions. Such long-range correlators, e.g., spin-spin 
correlation functions in a spin model, can be interpreted 
as the conventional Ginzburg-Landau order parameters 
of quantum phases. 

The meaning of “short-range physics” in this paper 
is more unconventional, by which we really mean how 
global symmetries are implemented locally in a quantum 
wavefunction. We will provide a sharper definition of 
this term later since we firstly need to introduce some 
tools to diagnose a local patch of the whole quantum 
wavefunction. But it is important to mention that this 
“short-range physics” is directly related to the quantum 
numbers of ground state wavefunctions on finite size sam¬ 
ples. In addition, both “short-range physics” and “long- 
range physics” in this paper are referred to properties of 
quantum wavefunctions, even in the absence of specific 
quantum Hamiltonians. 

As an interesting example, let’s consider candidate Z 2 
QSL that may be realized in the kagome lattice Heisen¬ 
berg model. Previous studies showed that there exist 
many time-reversal symmetric Z 2 QSL phases respecting 
the full space group symmetry of the kagome lattice^SHSll 
All these Z 2 QSL phases, by definition, are featureless in 
long-range correlators. So their distinctions completely 
lie in the short-range physics. 

The above discussions lead to the following intuitive 
picture. Different quantum phases may be organized into 
crude classes according to short range physics. In each 
class, there may be multiple member phases that are dis¬ 
tinguishable by long range physics. Although identifying 
a particular quantum phase in a correlated model gen¬ 
erally requires careful and challenging finite size scaling, 
identifying a crude class should be easier, even without 
finite-size scaling in large samples. In addition, doing 


the latter is still very useful. First, it would give us 
sharp, although incomplete, information about the quan¬ 
tum phase diagram. Second, determining the crude class 
allows us to focus only on the candidate member phases 
within one class, which helps identifying the complete 
phase diagram significantly. 

This picture motivates the us to separate the task of 
simulating the quantum phase diagrams into a short- 
range part and a long range part, and brings up the fol¬ 
lowing questions. How to systematically, and hopefully 
completely, characterize these crude classes distinguish¬ 
able by short-range physics? Can one construct generic 
variational wavefunctions for each given crude class and 
provide simulation algorithms based on them? The an¬ 
swers to these questions would lead to an efficient nu¬ 
merical method to completely solve the short-range part 
of the simulation task, which is very useful for the long- 
range part of the task as well. We will comment fur¬ 
ther on the sharp information on long-range physics (i.e. 
spontaneous symmetry breaking) that can be obtained 
from short-range physics in Sec |VIB' 

This paper is an attempt to address these questions to 
a certain level. Here we rely on a recently developed lan¬ 
guage to construct quite generic and physically relevant 
quantum wa vefunct ions: the projected entangled pair 
states (PEPSjP^Mol jg version of tensor networks. 
PEPS has been viewed as a powerful and efficient method 
to represent generic quantum states whose entanglement 
entropies do not violate the boundary law(see Ref.|3lj for 
a recent review). In addition, in two spatial dimensions, 
PEPS provides a set of concrete numerical algorithms 
for practical simulations (for instance, Ref|32] discusses 
details of many PEPS algorithms). In this work we con¬ 
struct generic symmetric wavefunctions using PEPS un¬ 
der certain assumptions. We find that there are classes 
of symmetric PEPS which are sharply distinguished by 
short range physics. More precisely, the symmetry re¬ 
quirements on PEPS lead to discrete number of solutions. 
Each solution corresponds to one crude class mentioned 
above, and constrains a sub-Hilbert space that the ten¬ 
sors in the PEPS must live within. 

We find that these classes are related to, but not lim¬ 
ited to, fractionalizations of both the on-site sy mmetr ies 
and the space group symmetries of the systenP^lMHH 
These classes are generally characterized by three sets 
of algebraic data, which are denoted as 0 ’s, y’s and 77 ’s 
in this paper. The first set of data (0’s) represents the 
direct contribution to the symmetry quantum numbers 
of quantum wavefunctions from each local tensor. The 
second set of data (y’s) is related to projective repre¬ 
sentations of the global symmetry, or the second coho¬ 
mology group H'^{SG,U{1)) in mathematics, where SG 
is the symmetry group of the system that is generally 
a combination of on-site symmetries and lattice symme¬ 
tries. The third set of data (r?’s) is related to the so- 
called projective symmetry group (PSGjP^ characterizing 
symmetry fractionalizations of topological quasiparticles. 
Mathematically, p's are related to the second cohomolgy 
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group H‘^{SG, IGG), where IGG is some invariant gauge 
group. We will explain the origin and constraints on IGG 
in detail later. Different possible IGG actually gives a hi¬ 
erarchical structure of the crude classes. As an example, 
half-integer spin systems on the square or kagome lattices 
have /GG’s which at least contain a Z 2 subgroup. 

Moreover, we provide concrete simulation algorithms 
based on these symmetric PEPS wavefunctions in two 
spatial dimensions (2d) and comment on possible algo¬ 
rithms in higher dimensions. We demonstrate the pro¬ 
cedure of crude classifying and constructing symmetric 
PEPS wavefunctions for the half-integer spin system on 
the kagome lattice, in which case 32 distinct classes are 
found under the assumption that IGG = Z 2 . Although 
we mainly consider 2d systems in this paper, the major¬ 
ity of our discussions can be easily generalized to other 
spatial dimensions except for the algorithms specific for 
2d. 

Not surprisingly, the choice of the kagome lattice spin 
system as the main example in this paper is motivated by 
the re cent r eports of a Z 2 QSL in the spin-1/2 Heisenberg 
modeP^*^. It remains an open question that which one of 
many candidate QSL may be realized in the kagome lat¬ 
tice modeP ^. An d very recently there have been a num¬ 
ber of work^^^EH describing how to identify these distinct 
Z 2 QSL in numerical simulations, based on careful quan¬ 
tum number analysis. In our work, when IGG = Z 2 in a 
half-integer spin system, every crude class contains a dis¬ 
tinct Z 2 QSL as a member phase. Therefore part of our 
results can be viewed as a classification and construction 
of Z 2 QSL for half-integer spin systems on the kagome 
lattice, which is somewhat fi ner th an the prev ious classifi¬ 
cations for the spin-I/2 cas^^ZMl^see Sec |VIl| for details), 
and is generally applicable for other half-integer spins. 
In addition, the simulation algorithms proposed here can 
be used to identify the nature of the Z 2 QSL realized in 
the kagome lattice spin-1/2 Heisenberg model efficiently. 

For each given crude class, the other member phases 
can be viewed as ordered phases in the vicinity of the Z 2 
QSL member phase, but with a spontaneous symmetry 
breaking only sharply defined in the long range physics, 
e.g. MO phases (via e-condensations) or VBS phases (via 
TO-condensations). The nonvanishing symmetry breaking 
long range order parameters in these phases are expected 
to be captured in the present symmetric PEPS contruc- 
tion after performing a scaling with repect to both the 
virtual bond dimension D (see Sec^^for definition) and 
system sizes. 

Note that the concepts of invariant gauge groups and 
projective symmetry groups(PSG) have been used to 
study and classify symmetry fractionalizations in topo¬ 
logically ordered phases. In this sense it is not surpris¬ 
ing that we find many non-symmetry-breaking Z 2 QSL 
phases distinct by short-range physics. But in a conven¬ 
tional symmetry breaking phase, such as the MO phases 
or VBS phases mentioned above, there is no topologi¬ 
cal order and the long-range gauge dynamics is confined. 
However, due to the generality of the PEPS language. 


this work suggests that the concepts of invariant gauge 
groups and projective symmetry groups are useful even 
in these conventional ordered phases. 

This interesting question raised by the present work 
can be rephrased in the following way. Do the neighor- 
ing conventional symmetry-breaking phases still “remem¬ 
ber” their parent non-symmetry-breaking liquid phase? 
In many situations the answer to this question is known 
to be positive. For example, consider the two parent 
Z 2 QSL, Sa chdev’s Qi = Q 2 state and Qi = —Q 2 
stat^^MUlI! After the spinon-e condensation, they lead 
distinct long-range MO phases, e.g. so-called q = 0 
MO (for the Qi = Q 2 QSL) and x -s/S MO (for the 
Qi = —Q 2 QSL). However, in some other situations, the 
answer to this question is expected to be negative. For 
instance, the vison-m condensation in these two Z 2 QSL 
could lead to the same VBS phas^^. This phenomenon 
is related to following fact: for the Qi = Q 2 state and 
Qi = —Q 2 state, the PSGs for the spinon-e are different, 
but the PSGs for the visons are the same. 

Therefore, within the framework proposed in this pa¬ 
per, although one non-symmetry-breaking phase only ap¬ 
pears in a single crude class, we cannot rule out the 
situation that certain special symmetry-breaking phase 
appears as member phases in multiple crude classes. 
Namely, it seems possible that two distinct short-range 
implementations of global symmetry lead to the same 
symmetry breaking phase in the thermodynamic limit. 
We will come back to this issue in Sec IVI Hi 

This work may be also useful regarding continuous 
quantum phase transitions. We have mentioned the 
phase transitions between member phases within one 
crude class, e.g. the transition between a Z 2 QSL and 
a nearby MO (VBS) phase. One may wonder whether 
it is possible to have a continuous phase transition be¬ 
tween two phases belonging to distinct crude classes. We 
believe that this is possible and is related to the hierar¬ 
chical structure of crude classes due to different JGGs. 
For instance, one may consider a parent crude class with 
IGG = U{1) that has two distinct descendent IGG = Z 2 
crude classes. Two phases belonging to these two distinct 
IGG = Z 2 crude classes, as a matter of principle, may 
be connected by a critical point described by the parent 
IGG = U{1) crude class. We leave further discussions 
on this topic in Sec |VH| 

Before moving on to the main body of the paper, we 
comment on the limitations of this work. First, due to the 
fact that we use PEPS to construct ground state wave- 
functions, the discussions in this paper is limited to those 
quantum phases whose entanglement entropies do not vi¬ 
olate the boundary law. For instance quantum phases 
with Fermi surfaces are beyond the scope of the current 
work. Even within the PEPS language our work makes a 
nontrivial basic assumption: the on-site symmetries are 
implemented as representations or projective representa¬ 
tions on the virtual degrees of freedom in PEPS. This 
assumption, although appears natural on the superficial 
level, is nontrivial and gives rise to limitations. 
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This problem is related to the recently developed un¬ 
derstandings on symmetry protected topological (SPT) 
phased. SPT phases are gapped quantum phases with¬ 
out anyon excitations and protected by various global 
symmetries. They are generalizations of the topolog¬ 
ical insulators(see Ref[49] and |50] for reviews) in non¬ 
interacting fermion systems. It is known that when at¬ 
tempting to represent SPT phases using PEPS, the con¬ 
straint that the on-site symmetry transforms as represen¬ 
tations or projective representations on virtual degrees 
of freedoms leads to problems, at least in the long range 
physics. For example a fermion state with nonzero Chern 
number constructed using PEPS with a fixed bond di¬ 
mension D under the above constraint is found to have 
power-law correlations in real spacd^i*!^. 


This paper is organized as follows. In Sec|TTj we in¬ 
troduce some basics of PEPS. In particular, We discuss 
gauge redundancy as well as the implementation of sym¬ 
metries in PEPS. We introduce a special kind of gauge 
transformation named as invariant gauge group {IGG). 
In phases with no symmetry breaking, IGG leads to 
low-energy gauge dynamics. Further, for fractional hlled 
systems, there are minimal required nontrivial IGGs for 
any symmetric PEPS under our basic assumption. This 
phenomenon is consistent with the Hast ings -Oshikawa- 
Lieb-Schultz-Mattis theorerrPlHSl! jn Sec III we classify 
symmetric PEPS according to their distinct short-range 
physics, which is characterized by algebraic data 0’s, y’s 
and ry’s. Relations of the data y’s and 77 ’s to second co¬ 
homology are discussed. And an introduction of relevant 
mathematics is given in AppendixAs a main example, 
we give the classification result for symmetric PEPS on 
the kagome lattice with a half-integer spin per site and 
IGG = Z 2 , and obtain the constraints on the sub-Hilbert 
spaces for local tensors for each given class. The detailed 
calculation is presented in Appendix A simpler and 
pedagogical example on the square lattice can be found 
in Appendix]^ We also give efficient algorithms for min¬ 
imization of energy density for a given class of PEPS, 
which can be used to identify these crude classes in in¬ 
teracting quantum systems. We give the physical inter¬ 
pretation of the algebraic data in Sec |IV| Particularly, we 
construct fractionalized symmetry operators to explicitly 
show that ry’s are describing the symmetry fractionaliza- 
tion of spinons in the Z 2 QSL member phase. Detectable 
signatures of the data 0’s, x’s and rj^s are discussed. In 
Sec|Vj we construct a decorated version for symmetric 
PEPS, which serves as a more convenient tool to study 
properties of visons in the Z 2 QSL member phase and 
the properties of the vison-condensed phases. Algebraic 
methods to extract the information of the symmetry frac- 
tionalization on visons are given. In Sec. |VI| we discuss 
symmetry-breaking phases in the symmetric PEPS for¬ 
mulation, and study the effects of the symmetry-breaking 
orders in finite-size scaling on torus samples. In Sec |VII| 
we consider generalizations and limitations of our study, 
comment on relations with previous works, and conclude. 


II. SYMMETRY, GAUGE AND PEPS 

In this section, we will give a brief introduction to 
PEPS. As we will see later, even for the same many- 
body wavefunction, the PEPS representations are not 
unique, and different representations are connected by 
gauge transformations. Further, we will study the im¬ 
plementation of symmetry on PEPS as well as the gauge 
dynamics in the PEPS language. Particularly, for certain 
systems, gauge structures will naturally emerge. 


A. Introduction to PEPS 

Projected Entangled Pair States (PEPS) is a type of 
tensor networks (TN). The basic ingredients of TN are 
“legs”, and every leg is associated with a Hilbert space, as 
seen in Fig.([^). In the following, we will use “leg” to de¬ 
note the associated Hilbert space. As shown in Fig([^), 
tensors formed by several legs simply describe quantum 
states living in the tensor product of these legs, 

j,a&c... (1) 

where V® labels Hilbert space associated with leg i. If 
two legs are the bra space and the ket space of the same 
set of quantum states, they are named as dual space to 
each other. New tensors can be obtained by contracting 
states in dual spaces, or by tracing out states in dual 
spaces, as shown in Fig. 

A TN representation of many-body wavefunction can 
be viewed as a large tensor, which is obtained by con¬ 
tracting small building block tensors. Thus, a TN is 
formed by uncontracted legs (physical legs) and con¬ 
tracted legs (virtual legs). From another point of view, 
we can also treat a TN as a combination of a linear map 
from the virtual Hilbert space (the tensor product of all 
virtual legs) to the physical Hilbert space, together with 
an “input” virtual state. 

Let us construct a PEPS on a two dimensional lattice. 
We first put tensors at both sites and bonds, named as 
site tensors (T®) and bond tensors (i?b) respectively, see 
Every site tensor can be viewed as a linear 
map from several virtual legs to one physical leg, while a 
bond tensor, which is in fact a matrix, labels a quantum 
state (bond state) in the tensor product space of two 
virtual legs. Thus, as shown in Fig. by contracting 
virtual legs of site tensors with bond tensors, we get a 
PEPS as a combination of a linear map from the virtual 
Hilbert space to the physical Hilbert space together with 
an input virtual state, where the map is given by the 
tensor product of all site tensors and the input state is 
the tensor product of all bond states. We can express the 
PEPS representation of the wavefunction as 

1^) = ^ tTr ((ri)'®F..(T^=)'=«=Hi...i?A,J |fci... 

{fcH 

( 2 ) 
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FIG. 1. (a): The leg a is associated with the Hilbert space 
V“. (b): The site tensor (left) and the bond tensor (right) 
label quantum states on Hilbert spaces of tensor products of 
corresponding legs, (c): A new tensor can be obtained by 
contraction of the leg b on T® and the leg a' on Bh, which can 
be expressed as ■ Note that we require 

leg b and leg a' to be dual spaces, (d): The whole PEPS 
wavefunction is obtained by contracting all virtual legs of site 
tensors and bond tensors. 


where 1,2,... Ns{N\^) label sites (bonds), while fcg is the 
physical index. tTr means tensor trace, namely, contrac¬ 
tion of all virtual legs. 

We define that a bond tensor (matrix) is a maximal en¬ 
tangled state, iff singular values of this matrix all equal 
some nonzero constant. By multiplying some constant, 
we can simply set singular values of maximal entangled 
states to be 1. When performing numerical simulations, 
it is more convenient to use maximal entangled bond 
states, or even set bond tensors to be identity matrices. 
As we will see later, by using the gauge redundancy of 
PEPS, it is always possible to do so. 

In the following, we will assume that all virtual legs 
label Hilbert spaces with the same dimension D, while a 
physical leg is associated with a d—dimensional Hilbert 
space. 


B. Gauge transformation on PEPS 

The representation of a many-body wavefunction on 
PEPS is far from unique. Particularly, as shown in 
Fig.§, we are always allowed to multiply W and W ^ 
to two connected virtual legs respectively. This action 
will change the connected small tensors while leaving the 
contracted tensor invariant, 

{nLcdSba'{B^)a'b' = [{ntbcdWbl]6u-[{W-^)l-a>iB^)a'b'] 

(3) 

Every contracted pair of virtual legs will contribute a 
gauge redundancy GL(Z),C). All such gauge transfor¬ 
mations form a group [GL(II, which we call the 

gauge transformation group of the PEPS {Nf, is the num¬ 
ber of bond tensors in the TN). The meaning of the gauge 
transformation can be understood as a change of basis on 
virtual legs. 

From another point of view, in general, for two PEPS 
whose tensors differ at most by gauge transformations de- 



Bi2 




W{l,b) 



-IV-y2.a) 


FIG. 2. Two PEPS describe the same quantum state, iff they 
are differ by gauge transformation together with U(l) phase 
factor. The origin of the gauge transformation is that we 
can multiply identity matrix I = W ■W~^ between connected 
legs, which changes site tensors and bond tensors, but leave 
the whole wavefunction invariant. We can also view TN on 
the left as PEPS transformed by symmetry operation. Thus, 
this figure also express the condition for PEPS wavefunction 
to be symmetric. 


fined above together with overall U(l) phase factors, as 
shown in Fig.(l^, the two PEPS must describe the same 
physical state (up a U(l) phase). In principal, these over¬ 
all U(l) phase factors can occur in gauge transformations 
on both site tensors and bond tensors. But it is straight¬ 
forward to redefine the gauge transformations such that 
the phase factors only appear on site tensors. Mathe¬ 
matically, two PEPS denoted by {T^,Bb} and {T^,Bb} 
respectively describe the same physical state if there exist 
gauge transformations {IT(s,i)} and U(l) phase factors 
{eif'(s)} (g labels a site and i labels a virtual leg on the 
site.), such that 

• [1F(S, l)].„^[IE(s, 2)]pp,... 

(Sb)a/3 = [bF(b, !)]„„'[VF(b, 2)]^^. (i3b)a'/3'. 

(4) 

Here W(h,j) represents a gauge transformation on the 
leg j of the bond tensor Hb, and if a site leg (s, i) 
and a bond leg (b,j) are connected, then IE(s, z) = 
[lE(b, j)“^]*. (The superscript-t stands for the matrix 
transpose.) 


C. Symmetric PEPS 

The purpose of this section is to introduce a generic 
way to implement both on-site symmetrie^^^^^ and lat¬ 
tice space group symmetries on PEPS. We firstly dis¬ 
cuss the finite size symmetric quantum state that can be 
represented by a single PEPS; i.e., such a state would 
form a one-dimensional representation of the symmetry 
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group. Then we define the symmetric PEPS on an in¬ 
finite lattice, which is the main object to be (partially) 
classified in the current study. 


1. On-site unitary symmetries 

The action of a global on-site unitary symmetry S on 
a finite size PEPS wavefunction is defined as 

S\^P) = 1^) = E tTr ... Bm,) 

{k.} 

Us®Us...\kik2...kN,), (5) 

Us is the representation of S on Hilbert space of physical 
leg. These local actions of an on-site symmetry give a new 
TN, with site tensors T® and bond tensors B\, defined as, 

i 

Hb = 5 o Hb = Hb (6) 

We focus on those PEPS that are invariant under the 
global symmetry up to an overall U(l) phase factor. Fol¬ 
lowing the discussion in the previous section, we consider 
the PEPS lip) that differs from the transformed PEPS 
\tjj) only by gauge transformations together with overall 
phase factors, as shown in Fig.([^: 

= 0sWs5 o 

Bb = WsS o Hb (7) 

Here, gauge transformation Ws and phase factor ©g as¬ 
sociated with symmetry S is defined as 

WsoT^ = [Ws(s, l)]o„-[fEs(s, 2)]pp ,... 

WsoHb = [Ws{h,\)]^AWsi>,‘^)hp'{B^)<.'0'- 

( 8 ) 

According to the definition of a gauge transformation, 
if site virtual leg (s, i) and bond leg (b,j) are connected, 
thenVEs(s,i) = [lTs(b, Further, we always choose 

Ws such that only site tensors transform with extra U(l) 
phase factors. Note that so far we do not require matrices 
on the leg (s,i) lTs(s,i) to form a representation of the 
on-site symmetry group when S is tuned. We will come 
back to this shortly. 


2. Time reversal symmetry 

The representation of the global time reversal symme¬ 
try T on a many-body wavefunction is Uj- ® Uj -... AT, 
where K denotes the complex conjugation and Uj- is a 


unitary matrix acting on local physical Hilbert space. Its 
action on PEPS is defined as 

nV-) = 51 tTr . .. BnS 

iU) 

Uj- ® Uj- . . . \k1k2 ■ ■ ■ kuP), ( 9 ) 

Namely, the local actions on a single site or a bond tensor 
read 

= r oT^ = Y,{UT)ki{Ty^ 

I 

i?b=roHb = H; (10) 

We consider the PEPS that is symmetric under T. Sim¬ 
ilar to the previous discussion, we consider a PEPS sat¬ 
isfying: 

T" = ©rlErToT" 

By, = W-jT o By, 

(11) 

where Wj- belongs to the gauge transformation group of 
the PEPS. 

3. Lattice symmetry 

The definition of a lattice space group symmetry R on 
PEPS is 

a/ 9 ... 

Hb = i? o Hb = ( 12 ) 

OL^ 

The action of R on site and bond tensor follows the nat¬ 
ural definition of lattice symmetries. For instance, for a 
square lattice, after a translation along the right direc¬ 
tion by one lattice spacing, the transformed site tensor at 
a given position equals the original site tensor on the left 
neighboring site. Note that the symmetry R not only 
acts on site and bond indices; it may also act nontriv- 
ially on virtual legs. For example, the 90° rotation of a 
site tensor on the square lattice permute the four virtual 
legs. Again, we consider those PEPS symmetric under R 
satisfying the following conditions: 

= QrWrRoT^ 

By, = WrR o By, 

(13) 

where Wr belongs to the gauge transformation group of 
the PEPS. 


4- Symmetric PEPS on infinite lattices 


D. Invariant gauge group and gauge structure 


Space groups of lattices are usually defined for infinite 
lattices. This is because for a finite size sample, the lat¬ 
tice symmetry group is a finite group whose group struc¬ 
ture is non-generic. In this paper, we will focus on PEPS 
on infinite lattices satisfying Eq. ( |7|11|13D under symme¬ 
try transformations. And we define such PEPS as sym¬ 
metric PEPS on infinite lattices, or simply as symmetric 
PEPS. They form the main object to be (partially) clas¬ 
sified in the current investigation. 


A natural question that arises at this point is: are 
symmetric PEPS defined above general enough to cap¬ 
ture ground states of quantum phases? Let us limit our 
discussion within those quantum phases whose entangle¬ 
ment entropies do not violate the boundary law so that 
in principle they may be represented as PEPS. 


Basically, we expect that the symmetric PEPS on infi¬ 
nite lattices defined above are capable to capture all non¬ 
symmetry-breaking liquid phases. After putting on finite 
lattices and performing a scaling with respect to both the 
bond dimension D and lattice sizes, we expect the sym¬ 
metric PEPS are also capable to capture the neighboring 
ordered phases of the liquid phases. Here by “neighbor¬ 
ing” (or “in the vicinity below), we mean that the sym¬ 
metry breaking in these phases is only sharply defined 
in the thermodynamic limit (namely, in the long-range 
physics). Note that we do not have a proof support¬ 
ing the statement above. Nevertheless we are not aware 
of any counterexamples, so at least it is a reasonable 
conjecture.^. 


Sometimes one is forced to use more than one PEPS 
to represent ground state quantum wavefunctions. Eor 
instance, in a quantum spin system with SU(2) spin 
rotation symmetry, this happens for the ferromagnetic 
phase, whose ground states form a large spin representa¬ 
tion. However, such ferromagnetic phases are not in the 
vicinity of any non-symmetry-breaking liquid phases. 


So far, we do not require the transformations matrices 
W’s on the virtual legs to form representations or even 
projective representations for the on-site unitary symme¬ 
tries and the time-reversal symmetry. As mentioned be¬ 
fore, such a requirement leads to difficulties to represent 
SPT phases in two and higher spatial dimensions. (And 
SPT phases are non-symmetry-breaking liquid phases.) 
Indeed, if one translates the ground states of the exact 
solvable models of SPT phases with on-site symmetries 
into the language of PEPS, one still finds PEPS satisfy¬ 
ing Eq.([7p^. But the transformation matrices VP’s form 
neither representations nor projective representations of 
the on-site symmetry group. 


Generally classifying symmetric PEPS defined here is 
a difficult task and we currently do not know how to 
solve. Next we will introduce the invariant gauge group 
for PEPS and will make further assumptions so that we 
could make progress on this difficult task. 


Among the gauge transformations, there is a special 
subgroup which we call the invariant gauge group {IGG). 
Note that generally a gauge transformation will leave the 
physical wavefunction invariant while transforming the 
site tensors and bond tensors nontrivially in a PEPS. 
However, by definition, the action of IGG elements on 
PEPS even leaves all site tensors invariant up to overall 
U(l) phases and all bond tensors completely invariant!^. 
So IGG can be viewed as the “symmetry” of the building 
block tensors with actions only on virtual legs. In the fol¬ 
lowing, we wil l see that IGG is directly related to gauge 
dynamic d^^ l ^^l l^. We will also give examples where non¬ 
trivial IGG’s emerge naturally in fractional filled systems 
under a basic assumption. 

Note that the collection of all gauge transformations 
that leave all site tensors invariant up to overall U(l) 
phases and bond tensors completely invariant forms an 
infinite group, which we denote as IGG. These gauge 
transformations satisfy Eq.Q with T® = = B^. 

Namely, a gauge transformation {IT(s,j)} is in the IGG 
of a PEPS formed by {T®, B^} iff it satisfies: 

intp... = • [fA(s, !)]„„. [W(s, 2)]^g,... 

(Sb)a/3 = [IE(b, !)]„„. [IE(b, 2)]^^, (5b)a'/3', 

(14) 

for certain U(l) phase factors Here W{b,j) rep¬ 

resents a gauge transformation on the leg j of the bond 
tensor B^, and if a site leg (s, i) and a bond leg (b, j) are 
connected, then W{s,i) = [IT(b, j)“^]*. 

Clearly, if certain gauge transformation {IE(s,i)} be¬ 
longs to IGG, then one can straightforwardly multiply 
U(l) phases x(s,*) to the IE(s, z)-matrices: {IT(s,j)} — >■ 
{IE(s, i) = x(s,f)IT(s,j)} and obtain another element in 
IGG, if xi^fi) = X*(s^*0 when (s, *) and (s',z') are the 
two virtual legs connected by one bond tensor. If we view 
the U(l) phase factors {x(s,i)} leaving the bond tensors 
completely invariant as a special kind of gauge transfor¬ 
mations, they form an infinite abelian subgroup in the 
center of IGG, which we denote as the x ~ group, since 
they commute with any gauge transformations. 

In general one should work with the infinite group 
IGG. In this paper, for simplicity, we define IGG as 
the quotient group: 


IGG = . (15) 

X — group 

In addition, we will mainly focus on the cases in which 
IGG is a simple finite abelian group In this situation, 
it is straightforward to show that IGG = IGG x x ~ 
group, indicating IGG is just a simpler way to express 
IGG. This also means that we could equally view IGG 
as a Zn subgroup of IGG. In particular, there exist a 
generator g G IGG, but g ^ x ~ group and g satisfies 
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= 1 where I is the identity gauge transformation — 
the do-nothing gauge transformation. 

Note that if IGG is a more complicated group, since 
the center extension with respect to y — group can be 
nontrivial, it is possible that IGG ^ IGG x y — group. 
In this situation it is better to directly work with IGG. 


1. IGG and gauge dynamics 

Here we will discuss the physical meaning of IGG. We 
use IGG = Z 2 as an example. The following discussion 
can be easily generalized to other IGG groups. 

First, let us clarify the action of Z 2 IGG on PEPS. Ev¬ 
ery virtual leg accommodates a representation of Z 2 = 
{!,(;}. Note that we do not require representations on 
different legs to be the same. However, we require two 
connected legs accommodate representations dual to each 
other, so that applying the g actions on connected legs 
is just a special gauge transformation. The nontrivial Z 2 
IGG element is an action of g on all virtual legs. Follow¬ 
ing the definition of /GG, all site tensors are invariant 
up to ±1 and all bond tensors are completely invariant 
under this action, as shown in Fig.([^). Further, it is 
straightforward to derive that any patch cut from PEPS 
is invariant up to ±1 under the g actions on boundary 
virtual legs, as shown in Pig.(|^). 

The physical meaning of IGG is related to the gauge 
dynamics. To see this, let us first review the Z 2 gauge 
theory. There are two phases in the Z 2 gauge theory: 
the deconfined phase and the confined phase. In the de- 
confined phase, the Z 2 gauge theory describes Z 2 topo¬ 
logical order (toric code). The low-energy excitations 
include four types of quasiparticles: the trivial particle 
I, the chargon e, the fluxon m and the bound state of 
chargon and fluxon / = em. e, m and / can only be 
created in pairs. Each particle is its own anti-particle, 
= /^ = I. e, m are bosons while / is a fermion. 
The braiding statistics of the three nontrivial particles 
are mutually fermionic. In the confined phase, topologi¬ 
cally nontrivial quasiparticles are confined. 

To see the connection between IGG and the gauge the¬ 
ory, let us create nontrivial excitations on PEPS with Z 2 
IGG. We can define e particles living on sites while m 
particles living on plaquettes. As shown in Fig.([^), to 
create two m particles in neighboring plaquettes, we sim¬ 
ply multiply the nontrivial Z 2 element g on one of two 
contracted virtual legs shared by the two plaquettes. The 
insertion of g only on one side of contracted legs is not 
a gauge transformation, and in general will change the 
wavefunction. One can also create a pair of m parti¬ 
cles spatially separated from each other by applying the 
single-sided g-actions over a string of bonds. The fluxons 
are located at the end of the string. Note that although 
the positions of fluxons are physical, the position of the 
string connecting them are not physical since one can per¬ 
form Z 2 gauge transformations on site tensors to move 
the string around while leaving the physical wavefunction 


invariant. 

Now, let us turn to e particles. Let us first define Z 2 
even/odd tensors. The action of g on boundary virtual 
legs of a tensor generally gives a phase factor ±1. If the 
phase factor is -1-1/—1, we call it Z 2 even/odd. The Z 2 
parities of tensors depend on the representations of g on 
virtual legs. If we do not worry about the lattice sym¬ 
metry for the moment, for a Z 2 even/odd tensor, we can 
simply redefine g on one virtual leg by —1, thus this ten¬ 
sor becomes Z 2 odd/even. So we can assume all tensors 
are Z 2 even for the remaining discussion in this subsec¬ 
tion. Creating an e particle on a single site corresponds 
to changing the site tensor from Z 2 even to Z 2 odd, as 
seen in Fig.([^). To detect the number of chargons on a 
patch of PEPS, we simply apply g on all boundary vir¬ 
tual legs; namely, we create an m loop on the boundary. 
If there is an odd number of chargons on that patch, this 
patch tensor should be Z 2 odd and the g action on the 
boundary picks up a —I, see Pig.([^). This —1 can be 
understood as the Berry phase from braiding e and m. 
One can easily convince oneself that an odd number of 
chargons cannot be created on a closed manifold. 

If IGG = Z 2 PEPS describe deconfined phases, then 
separating topological quasiparticles is expected to cost 
zero tension. Consequently one can insert m loops wrap¬ 
ping around torus holes to construct the four-fold degen¬ 
erate ground states on a torus. However if IGG = Z 2 
PEPS describe confined phases, which we expect to be 
possible after a scaling with both bond dimension D and 
system sizes, this is no longer true. We will comment 
further on this in Sec IVIl 

As a final remark, there turns out to be two dis¬ 
tinct types of Z 2 gauge theo ries: th e toric code theory 
and the double-semion theorjEnZnEI]_ They have distinct 
topological orders; e.g., the topological spins (the ex¬ 
change statistics phases) of quasiparticles are [1,1,1, — I] 
i]) for the [I,e,m,em] particles in a toric code 
(double-semion) topological order. We emphasize that 
the IGG = Z 2 PEPS discussed here, when describing 
a deconfined phase, hosts the toric code topological or¬ 
der. The simplest way to see this is to realize the self 
braiding statistics phases of both the e and the m in the 
IGG = Z 2 PEPS are trivial, so they cannot be semions. 

Indeed, when moving an e chargon around a loop by 
a sequence of hoppings, one realizes the Berry’s phase 
is independent of whether there are other e chargons 
inside the loop. Similarly, when moving an m fluxon 
around a loop (giving rise to an m loop), the topological 
Berry’s phase is simply ±1 depending on the Z 2 parity of 
the PEPS patch inside the loop, independent of whether 
there are other m fluxons inside the loop. 


2. Natural emergence of nontrivial IGG 

We will show that, under a basic assumption, the sym¬ 
metric PEPS for certain quantum systems must have 
nontrivial /GG’s. This basic assumption is that the W 
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FIG. 3. (a): Site tensor and bond tensor are both invariant 
under Z 2 action on all virtual legs of tensors, (b): Tensors 
obtained by contracting Z 2 invariant tensors are also Z 2 in¬ 
variant. (c): Acting g on one virtual leg of single bond tensor 
creates two fluxons (m) in plaquettes sharing the bond, (d): 
Z 2 odd tensor indicates there sitting a chargon. (e): By ap¬ 
plying g (or creating fluxon loop) on the boundary of a region, 
we are able to determine chargon number is even or odd inside 
this region. 


matrices on every virtual leg form (generally reducible) 
representations or projective representations for the on¬ 
site symmetries (see Eq|7[TT|). Under this assumption, 
the nontrivial IGG in certain systems is a natural con¬ 
sequence of the global symmetry, even in the absence of 
specific Hamiltonians. 

Consider a spin-^ system on a square lattice; i.e., the 
physical leg on every site tensor is a 2 -dimensional spin- ^ 
Hilbert space. For this system, we will show a symmet¬ 
ric PEPS under the basic assumption must feature an 
IGG containing a Z 2 subgroup. Since SU{2) spin rota¬ 
tion group has no projective representations, the basic 
assumption ensures that every virtual leg must form a 
representation of SU{2), which generally is a direct sum 
of a number of half-integer spin representations and a 
number of integer spin representations. Eq.Q now has 
the following simple interpretation: the site tensors are 
spin singlets formed by the virtual spins and the physical 
spin-i, and the bond tensors are spin singlets formed the 
virtual spins only. 

Now we can consider the particular 27r SU{2) rota¬ 
tion, and denote the corresponding lU(s,f) matrix on a 
virtual leg (s, i) as J(s, i), which is simply a direct sum of 
the minus identity transformation in the half-integer spin 
subspace and the identity transformation in the integer 
spin subspace. Next, consider the combination of trans¬ 
formations {J(s, *)} acting on the virtual legs only — this 
is a particular gauge transformation. Since the physical 
spin-i only picks up an overall —1 in the 27r SU{2) ro¬ 
tation, and the bond tensors are spin singlets, we know 
that the gauge transformation {J(s,i)} is an element in 

JUG. 


To see this system featuring a nontrivial IGG, we only 
need to show IGG X~ group. We will demonstrate 
that the gauge transformation {J(s,f)} ^ x ~ group. To 
do this, we impose the G 4 rotational symmetry and the 
translation symmetry of the square lattice. Note that 
{J(s, f)} G X ~ group if and only if for every virtual leg, 
the dimension of either the half-integer spin subspace or 
the integer spin subspace vanishes. However, this cannot 
be true. The site tensor is a spin singlet, which requires 
the virtual legs to combine into a spin- ^ so that it can fur¬ 
ther combine with the physical spin-f to form a singlet. 
Therefore, if {J(s, i)} € x~group, on a single site tensor, 
we must have an odd number of virtual legs which con¬ 
tain purely half-integer spins while the remaining virtual 
legs contain purely integer spins. This explicitly breaks 
the Ci rotational symmetry. 

Consequently, there is at least one element J = 
{J(s, I)} in IGG but not in x — group, and = e. This 
tells us that IGG at least contains a Z 2 subgroup {I, J}. 

The above argument can be easily generalized to other 
symmetries, such as the time reversal symmetry. For 
the time reversal symmetry, consider a system with 
one Kramer doublet on every physical leg. To form a 
Kramer singlet PEPS, one must combine an odd number 
of Kramer doublets on virtual legs of every site tensor. 
However, for site tensors on a square lattice, there are 
even number (four) of virtual legs per site, and the G 4 
symmetry dictates that the transformation on virtual 
legs only gives a nontrivial element of the IGG which is 
at least Z 2 . 

We point out that translational symmetry itself is 
enough for the above argument and one does not nec¬ 
essarily consider G 4 . This is because translational sym¬ 
metry relates the left (down) virtual leg with the right 
(up) virtual leg connected to the same site tensor via the 
fact that the virtual legs connected by a bond need to 
form a spin singlet (or a Kramer singlet). What is really 
important for the above argument is the existence of a 
half-integer spin (or a Kramer doublet) per unit cell. One 
way to see this is to consider a honeycomb lattice with 
spin-f per site, i.e., two spin-Us per unit cell. In this 
case, every site has three virtual legs and it is possible 
to construct symmetric PEPS wavefunctions with purely 
half-integer spins on virtual legs, in which case the 27r 
spin rotation on the virtual legs only becomes an element 
in the x ~ group. 

Next, let us consider a system with fractional filled 
hard core bosons and see how a nontrivial IGG naturally 
emerges. As an exercise, we can simply translate the 
previous discussions on spin-f systems into ^-filled hard¬ 
core boson systems on the square lattice. The physical leg 
for the hard-core bosons is two dimensional Hilbert space 
with basis labeled as |0) and |1). When mapped to a spin- 
5 system, | 0 )(| 1 )) is identified as the down spin (up spin). 
The U(l) charge transformation for the hard-core boson 
system can be written as exp[i0(S'* -I- |)] using the spin 
operator on the leg-i. Note that spin-0 is identified as 
charge-^, a projective representation of the charge U(l). 
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Since a bond tensor is a spin singlet formed by two virtual 
spins in the spin language, the representation of U(l) 
group in the hard-core boson language on a bond tensor 
is 


[eie(s:+5) . = e"'® (16) 

where the complex conjugation comes from the fact that 
bond virtual legs transform as conjugate representation 
of site virtual legs, and we have used = 0 for the 

two virtual legs a and b. So, every bond tensor carries 
charge —1. 

Further, since the site tensor is also a spin singlet, we 
require = 0, where i = 0 labels the physical 

leg and other i ^ 0 label virtual legs. Therefore the 
representation of U(l) symmetry on a site tensor reads 

+ =e‘S^ (17) 

i=0 

Namely, every site tensor carries charge-1. Consequently 
each unit cell carries charge- 

Note that in this exercise, the bond tensor transform 
nontrivially under U(l), so the virtual leg transformation 
W = does not satisfy Eq.Q in our definition 

of symmetric PEPS. But one could easily redefine the 
virtual leg transformation VF’s so that the charge carried 
by the bond is absorbed to a neighboring site, and Eq. 0 
is satisfied using the redefined IP’s. 

The essential results from previous discussions on the 
spin-i systems can now be translated as following state¬ 
ment: the virtual leg hosts both integer charges and half 
integer charges of U(l), so 2tt rotation of U(l) symmetry 
on all virtual legs gives the nontrivial Z 2 IGG. 

In the following, on the square lattice, we provide 
a general argument that a nontrivial minimal required 
IGG emerges for a symmetric PEPS with fractional-filled 
bosons under our basic assumption. Further, this mini¬ 
mal required IGG is given by the 27r rotation of the U(l) 
symmetry on the virtual legs only. 

Firstly, we have the physical legs carrying integer 
charges. And if the tensor network is symmetric under 
the U(l) symmetry, for site tensors and bond tensors, we 
can rewrite Eq. 0 as 


WsS o T" = ©sT" 

WsS o Bb = ©s^b (18) 

where symmetry operation S can be any U(l) group ele¬ 
ment. Note that we put ©g operation on bond tensors as 
well to pick up the possible phase factors. As mentioned 
before, this phase factor on the bond can always be tuned 
away by redefining Wg. But for the moment, let us keep 
it since we want to include the previous exercise. 

We can view the left side as the U(l) action on a 
site/bond tensor. Under the basic assumption, the above 
equation indicates every site/bond tensor carries a fixed 
U(l) charge, which can be a fractional charge. In the 


presence of the lattice symmetry, we expect all virtual 
legs of site tensors share the same U(l) reducible projec¬ 
tive representation. (Virtual legs of bond tensors have 
the conjugate representation). Our plan is to assume the 
27r rotation of U(l) symmetry is trivial (only a phase 
factor) on the virtual leg, and then demonstrate a con¬ 
tradiction. This assumption dictates that the irreducible 
charges carried by a virtual leg can only differ by integer 
numbers. Namely, the basis for virtual legs of site tensors 
can be written as 


{\x),\x + ni),\x + n 2 ),...} (19) 

where x can be any fractional number and rii are integers. 
Under symmetry operation Ug, state \x-\-i) transform as 

[/e|a:-I-rij) = -I-rii) (20) 


So, 2tt rotation on any state of the above Hilbert space 
will give the same phase factor e'^®. Similarly, the basis 
for bond legs are 

{| - a:), I - a; - m), I - X - 712 ),...} (21) 


Recall that a single tensor should carry a fixed charge. 
Consequently a bond tensor should carry charge — 2x — 
ng, where ng is some integer. And a site tensor should 
carry charge Ax + Ug- Since the physical leg only carries 
integer charges, should also be an integer. We then 
conclude that, for a single unit cell, the charge should 
be nig — ng, which must be an integer. This contradicts 
with the fact that the system is at a fractional filling. 
Therefore to construct a symmetric PEPS at a fractional 
filling under our basic assumption, the 27r rotation of 
U(l) symmetry must be nontrivial on all virtual legs, 
and the nontrivial IGG naturally emerges. 

We discussed the naturally emerged IGG in certain 
quantum systems. It is possible for the ground state 
symmetric PEPS to have a larger IGG which contains 
the naturally emerged IGG as a subgroup. We call the 
naturally emerged IGG as the minimal required IGG. A 
larger IGG than the minimal required IGG has impor¬ 
tant implications in both conceptual understandings and 
numerical simulations. We will come back to this point 
in Sec.(Ill) and Sec.(VII). 

The minimal required /GG’s in systems at fractional 
fillings are consistent with the Hastings-Oshikawa-Lieb- 
Schultz-Mattis (HOLSM) theorem. Consider a 2-l-lD 
system with an odd number of spin-^ per unit cell, the 
HOLSM theorem states that it is impossible to have a 
featureless trivial insulator. In other words, the ground 
state must either be gapless, break the spin rotation or 
the lattice translation symmetry, or be topological or¬ 
dered with a ground state degeneracy. 

In our formalism, a half-integer spin per site on the 
square lattice (and similarly on the kagome lattice) en¬ 
forces a minimal Z 2 IGG, consistent with the HOLSM 
theorem. Eor instance, if IGG = Z 2 , the system could be 
in either a deconfined phase with a toric code topological 
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order, or a confined phase. But the confined phase cor¬ 
responds to either e or m condensation, which leads to 
spin rotation or lattice translation symmetry breaking. 

For a honeycomb lattice spin-^ system, there are two 
spin-i per unit cell and the HOLSM theorem does not 
apply. As mentioned above, we expect that symmet¬ 
ric PEPS on the honeycomb with a trivial IGG can be 
constructed, which is consistent with the possible trivial 
symmetric insulator phase in this system as pointed out 
in Ref. [71]. 


E. An example 

Here, we will give a simple PEPS with IGG = de¬ 
fined on the kagome lattice. In particular, we will write 
the PEPS description for a nearest neighboring (NN) 
resonating valence bond (RVB) state that preserves all 
lattice symmetry. The lattice symmetry generators for 
kagome lattice are shown in Fig.Q. 

As shown in Ref.l^, there are four different kinds of 
symmetric NN RVB states defined on kagome lattice with 
spin-i per site. Also, by solving projective symmetry 
group (PSG) equations for the Schwinger-boson mean 
field ansatz on the kagome lattice, one finds eight distinct 
PSG classes. And four of them can be realized by NN 
pairing term^^. One can check that the four NN RVB 
states are exactly representative states for these four PSG 
classes. Here, we will focus on on e part icular PSG class, 
named as Qi = Qi state in Ref.ESEZi, This pa, rticular 
PSG class is a promising candidate phase^^^^^^^ for the 
spin liquid reported in recent DMRG simulation^i^lli^l. 
Here, we will explicitly write down this NN RVB state in 
the PEPS language. 

In fac t, this state has already been studied extensively 
in PEP^MZH. Here, we will slightly modify the con¬ 
struction. Every physical leg is a spin-f and virtual 
leg accommodates spin representation 0 0 i, with ba¬ 
sis {|0), I t)) I i)}- Bond tensors are spin singlets, which 
can be written as a matrix in this basis, 

1 0 0 \ 

0 0 -i (22) 

0 i 0 / 

where the direction of bond tensor is shown in Fig.(|^). 
A bond tensor with the inverse direction is transpose of 
the above matrix. Tensors for different sites are equal to 
each other, and can be written as, 




FIG. 4. (a): kagome lattice and the elements of its symmetry 
group. ai ,2 are the translation unit vectors, C(, denotes tt/S 
rotation around honeycomb center and cr represents mirror 
reflection along the dashed red line, (b): Site tensor and 
bond tensor for kagome lattice in one unit cell. Virtual legs 
of site tensors are labeled as {x,y,s,i), where {x,y) denotes 
the position of unit cell, s = u,v,w is the sublattice index 
and i = a,b,c,d specifies one of four legs, (c): One possible 
orientation of kagome lattice. Particularly, for NN RVB state, 
the orientation of bonds denotes the direction of spin singlets. 


chosen to make PEPS symmetric under lattice symme¬ 
tries. One can verify the state defined above is consistent 
with the PEPS representation of NN RVB given in Ref.^^ 
up to a gauge transformation. 

As discussed before, the Z 2 IGG here is generated by 
the 27r spin rotation of all virtual legs. Since all tensors 
are spin-singlet, they are invariant under this operation 
up to —1 factors on the site tensors. This NN RVB PEPS 
belongs to one of the crude classes proposed in this paper. 
Roughly speaking, according to global symmetry, we can 
find the generic sub-Hilbert space that the building block 
tensors must live within for each given crude class, which 
vastly generalize the one-dimensional sub-Hilbert space 
defined as in Eq.(23). 


III. ALGORITHM FOR SYMMETRIC PEPS 


T" =1 t) G) (U 000) + |0 i 00) - i|00 j. 0) - i|000 j.))- 
I 4 ,) 0 (I t 000) + |0 t 00) - i|00 to) - i|000 f)) 

(23) 

where the order of site virtual legs is given in Fig.(|^). 
We can view site tensors as superposition of singlets 
formed by one physical leg and one of the four virtual 
legs, while the coefficient of singlets need to be carefully 


For a given quantum model with certain given symme¬ 
try groups, we propose a general simulation scheme to 
study its phase diagram as follows: 

1. One classifies symmetric PEPS according to their 
short-range physics. More precisely, crude classes 
are distinguished by ways of implementing symme¬ 
tries on virtual legs. 
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2. For each class, by enforcing symmetry transforma¬ 
tion rules, one finds constraint Hilbert spaces for 
the building block tensors in the PEPS representa¬ 
tion. 


3. One performs the energy density minimization for 
every class in the constrained Hilbert space, and 
determines the class which gives the lowest energy 
density. The quantum phase of the model will be a 
member phase of this crude class. This finishes the 
short-range part of the simulation task. 


4. At last, one could try to completely determine the 
quantum phase diagram by studying the long-range 
physics, e.g., by measuring correlation functions for 
the symmetric PEPS with the minimal energy den¬ 
sity. With a careful scaling analysis, together with 
the sharp information on the long-range physics ob¬ 
tained from the short-range physics (see Sec VI for 
details), possible long range symmetry breaking or¬ 
ders may be identified. 


As the main example, we will demonstrate this simula¬ 
tion scheme for a half-integer spin system on the kagome 
lattice. We will start with classifying and constructing 
generic symmetric PEPS with IGG = that preserve 
the full lattice symmetry as well as the spin rotation and 
the time reversal symmetries. As we will show shortly, 
the condition IGG = Zi actually dictates that the virtual 
legs form (projective) representations of on-site symme¬ 
tries. Therefore when we consider IGG = Z^ symmetric 
PEPS, we already made our basic assumption in an im¬ 
plicit way. In addition, although we focus on the minimal 
required IGG under our basic assumption, the discus¬ 
sions can also be easily generalized to symmetric PEPS 
with a larger IGG. 


A. General framework for classification 


From now on we assume IGG = IGG x y — group, 
which is always true if IGG is a simple finite abelian 
group Zn- 

Consider the gauge transformation associated with a 
symmetry R: Wn, and the corresponding phase on site 
tensors: 0^;. We have T® = QrWrRoT^ and Hb = Wro 
Hb, as shown in Sec |H C| However, since both site tensors 
and bond tensors are invariant under the IGG action (up 
to phases for site tensors), we conclude that tensors are 
also invariant under a new symmetry operation defined 
as W'^ = -nnWR and = prQr, 

T" = Q'rWrR o T" 

B^, = W’rRoB^, (24) 

where r]R € IGG and pR = {pr{s)} is a set of phase 
factors on site tensors associated with rjR, such that 
fJ-RVR o r® = T®. For instance, for a half-integer spin 
system described by PEPS with IGG = {I, J}, if 77 _r = J 


corresponds to the 27r SU (2) rotation on the virtual legs, 
then pr{s) = — 1 for all sites. 

Similarly one could modify Wr and Qr with any ele¬ 
ment in the x~ 9 foup, i.e., bond dependent phase factors 
{£fl(s,i)} as: 

WR{s/h,i) eR{s/h,i)WR{s/h,i) 

0i?(s) 

i 

where we have £R{s,i) = £ij(b,j)* if (s,i) and (b,j) are 
connected. Further, £fl(b, 1) = £_R(b, 2)* for the two legs 
of the same bond tensor, as required in the definition of 
the X ~ group. 

Basically, the symmetry transformation on the virtual 
legs Wr is ambiguous since it can be combined with 
any element in IGG. Mathematically, the representation 
of R on the Hilbert space of PEPS (including both the 
virtual and physical Hilbert spaces) form a new group, 
which is the original symmetry group SG extended by 
the IGG. This extension is related to the 2-cohomology 
H‘^{SG,IGG) and H‘^{SG,U{1)). (For details about 
projective representations and the 2-cohomology, see Ap¬ 
pendix]^) Particularly, we can view those IGG elements 
as “representations” of the identity element in the sym¬ 
metry group on virtual legs. 

Keeping these discussions in mind, let us consider a dis¬ 
crete symmetry group SG as an example. S'G is always 
defined by a collection of group identities. For instance, 
elements i?i, i? 2 , ■ • •, -Rn € SG satisfy the following rela¬ 
tion: 


R1R2 . . . Rn — C (^ 6 ) 

Then, acting R1R2 ... on a symmetric PEPS, one ob¬ 
tains a combined transformation sending every tensor 
back to the same tensor: 

= Qr.Wr.RiQr.Wr.R^ ... o r® 

B^ = Wr,RiWr, i?2... Wr„ i?„ o Hb (27) 

By definition, the transformation leaving all tensors in¬ 
variant (up to phases on site tensors) can only be an ele¬ 
ment in IGG. Explicitly writing down Eq.( |27[ ) on virtual 
legs of site tensors, we conclude that 

WR,is,i)WR,iR^\s,i))... 

IVi?„(i?;(li...i?C\s,i)) =r]{s,i)x{s,i) (28) 

where ri{s,i) is the action of g G IGG on the virtual 
leg (s,i). Further, {x(s,z)} is an element in the % — 
group. We point out that since WR{s,i) = [IE^^(b,j)]* 
if (s, i) and (b,j) are connected, Wr on virtual legs of 
bond tensor gives us no extra equation. However, phase 
factors on site tensors will give an extra condition, which 
reads 

0fli(s)0fl2(^rHs))... OrJR-^ . .. AC^(s)) 

= M(s)n^*(®’*) (29) 
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Here /r*(s) is the phase factor obtained after applying 77 
on the s-site tensor. 

Our goal is to solve Eq.(|28|) and Eq.p^ for all group 


identities and obtain the representations of symmetry op¬ 
eration on virtual legs {Wr) as well as phase factors on 
site tensors {Qr)- Recall that the same physical wave- 
function can be represented by many PEPS which dif¬ 
fer from each other by gauge transformations (note that 
these are general gauge transformations wh ich may not 
be in / GG.). One should really solve Eq. (28) and Eq. (29) 
up to gauge equivalence. 

Under a gauge transformation U = {U(s, i)} on virtual 
legs, (r®)' = V oT® and B^ = VoB^ satisfy the following 
conditions: 


(T^y = vQrWrRot^ 

= {VeRV-^)(yWRRV-^R-^)RV oT^ 

= eRwyRoiT^y, 

(30) 


and 


B'^ = VWrRoB^ 

= {VWrRV-'^R-^)RV o Hb 


= wyRB'^. 


(31) 


Here we use the fact that V commutes with Or in the last 
step of Eq.@. Here, = VWrRV-^R-\ Writing 
the above expression explicitly on virtual leg (s, i), we get 

WR{s,i) V(s,i) ■ WR{s,i)V~\R~^{s,i)) (32) 

while Or is invariant. Particularly, rj G IGG changes as 

V(.s,i) V{s,i) ■ r]{s,i)V~\sG) (33) 


And phase factors /i and y in Eq.(29) are invariant. 


Apart from the above gauge transformation, one can 
change site tensors by phase factors, which do not affect 
physical observables. Note that one could also change 
bond tensors by phase factors, but such a modification 
is always equivalent to a gauge transformation together 
with a changing of phase factors on site tensors. Unlike 
gauge transformations, a modification of phase factors 
on site tensors may change the physical wavefunction up 
to an overall phase. When site tensors change as T® —>• 
$ o T® = <i)(s) • r® = Wr associated with the 

symmetry R is invariant, but Or goes to ^OrR^~^R~^. 
Namely, the phase factor Or = will change as 


0fl(s)^0i?(s)$(s)$*(i?-'(s)) 


(34) 


Basically, we should solve for the Wr and Or in 
Eq.(28l and Eq.(p^ up to two kinds of equivalences. 


First, if two sets of Wr and Or are related by Eq.(32) 


and Eq.(34), they are equivalent and we denote this sit¬ 
uation as the gauge equivalence. The gauge equivalence 
contains the U-ambiguity in Eq.( 32) and the ^-ambiguity 
in Eq.(34). 


Second, if two sets of Wr and Or are different by an 


IGG element, they are also equivalent and we denote 
this situation as the group extension equivalence. Sum¬ 
marizing our discussion in Eq.( 24|2'5 l, it means that one 
could modify Wr and Or as Wr -G = PrSrWr 
and 0i^ —?► 0)^ = PrSrOr, where tjr € IGG and 
Er G X ~ group and 


WyisG) = ?7i^(s,^)ei^(s,^)Wi^(s,^) 

0fl(s) =Mit(s)n^fl(®u)®(®)' 


(35) 


Note that to save notation, we define erOr as multiply¬ 
ing n.4(su) on 0(s). The group extension equivalence 
contains an ? 7 -ambiguity and an £-ambiguity in Eq. (351. 


Note that different from the gauge equivalence, we have 
an rj-ambiguity and an E-ambiguity for each symmetry 
element R. 


We will solve Eq. (28 1 and Eq. (291 for the whole sym¬ 


metry group up to both the gauge equivalence and the 
group extension equivalence. Eventually we will obtain 
many classes of PEPS satisfying inequivalent Wr and Or 
transformation rules. Among all combinations of Wr and 
Or within the same equivalence class, we can choose a 
particular representative, and construct explicit forms of 
Wr and Or by fixing the ? 7 -ambiguity, the ^-ambiguity, 
the U-ambiguity and the ^-ambiguity. These Wr and 
Or specify the sub-Hilbert spaces for the building block 
tensors in each class. We sometimes call the whole pro¬ 
cedure of hxing the four ambiguities as gauge fixing. 

Practically, we often firstly use the group extension 
equivalence to simplify Eq.(28l and Eq.([2^. For in¬ 
stance, one can use the e-ambiguity to simplify {x(s,z)} 


in Eq.(28) and Eq.(29): under a transformation Wr^ 


£Ri Wr. , according to Eq.(28 1 , we find 


Xis,i) -G ER^{s,i). ..er^{R„\ .. .i?i \s,z))x(s,f). 


(36) 


Moreover, one can use the ry-ambiguity to simplify the 
{ 77 (s, i)} and {/r(s)} in Eq.(28) and Eq.(29). For example, 
if some symmetry operation R appears only once in the 
group identity = e, one could use the 77 - 

ambiguity for R to make sure { 77 ( 3 , i) = 1} and {^(s) = 1} 
for this group condition. 


After the group extension equivalence is used, we will 
use the gauge equivalence (the U-ambiguity and the d*- 
ambiguity) to solve for explicit forms of Wr and Or. 
Note that the group extension equivalence and the gauge 
equivalence are not completely independent. For exam¬ 
ple, after fixing the U-ambiguity and the ^-ambiguity, 
it is possible some part of the e-ambiguity and the 77 - 
ambiguity are also fixed. In the following we demonstrate 
this procedure in an example: the half-integer spin sys¬ 
tems on the kagome lattice. 
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B. Classification of kagome PEPS 

Here, we will classify symmetric kagome PEPS wave- 
function with a half-integer spin-S' per site, which pre¬ 
serves all lattice symmetries, the time reversal symmetry 
as well as the spin rotation symmetry. We will only as¬ 
sume IGG = Z 2 = {I, ■}} without specifying the physical 
meaning of J. Later we will prove that J can always be 
chosen to be the 2tt spin rotation on the virtual legs. Let 
us begin with setting up some useful facts. 

First, we can use the H-ambiguity to diagonalize 
J(x, y, s, i) for every virtual leg (x, y, s, i), where (x, y, s) 
labels a site on the lattice by the coordinates of the 
unit cell x,y and the sublattice index s = u,v,w, and 
i = a, 5, c, d labels one of the four virtual legs coming 
out of the site tensor, (see Fig|3for illustrations) In this 
gauge, V(x, y, s, i), the matrix J(x, y, s, i) is a direct sum 
of an identity matrix and a minus identity matrix. Let 
us denote J(xo, yo, sO) * 0 ) = Idi ® {~Ii32) for some given 
virtual leg (xq, yo, Sqj *o)j where Di + D 2 = D. We will 
consider the generic case in which Di ^ Z? 2 - 

Using the lattice symmetry, it is straightforward to 
prove that one can always redefine {J(x, y, s, i)} by mul¬ 
tiplying with an element e in the x~g'>^oup: e{x, y, s, i) = 
±1 so that J(x, y, s, i) = Idi ® (— 11 ) 2)5 y, s, i). (Such 
a modification is allowed in our definition of IGG.) For 
example, consider a particular lattice symmetry opera¬ 
tion R, which could be the 60° degree rotation Ce or the 
lattice translation Ti or T 2 of the kagome lattice (see Ap¬ 
pendix]^ for precise definitions), we always have a group 


relation R ^ ■ e ■ R = e. Using Eq.(28) for this group 


relation and choosing J to replace the e on the LHS: 


Wj, (i?(x, y, s, i))I{R{x, y, s, f))IUfl(i?(x, y, s, i)) 


=? 7 (x,y,s,j)x(x,y,s,i). 


(37) 


The t] on the RHS must be J, otherwise we would find J 
to be an element in the x ~ group, violating IGG = ^ 2 - 
Therefore we know that J(i?(x, y, s, *)) and J(x,y, s,*), 
which are generally on two different virtual legs, are 
related by a similarity transformation IFA(i?(x, y, s, i)) 
and an overall phase factor x(x,y,s, i). But we are al¬ 
ready in a gauge such that J(x,y, s,i) are all diagonal. 
We then conclude that J{R{x,y, s,i)) = ±J{x,y,s,i). 
Since all virtual legs are related by lattice symmetries, 
we know J(x,y, s,i) = £(x, y, s, i)J(xo, yo, sqj * 0 )) where 
£{x,y,s,i) = ±1 V(x,y, s,j). 

Next, we show {e(x,y, s,i)} Gy — group. Namely, if 
(x, y, s, i) and (x', y', s', i') are connected by a bond ten¬ 
sor Bi,, then e(x, y, s, i) = e{x', y', s', i'). This is because 
if £(x,y,s,i) = —£{x',y',s',i'), then the matrix {Bii)ai 3 
satisfying Eq. (14) for lU = J would not have a full rank. 


since D\ ^ D 2 . This means that some singular value of 
(Hb) vanishes, dictating an IGG larger than Z 2 . For in¬ 
stance, one can multiply an arbitrary U(l) phase on the 
zero singular value eigenstate on one of the two virtual 
legs, leaving the bond tensor Bi, invariant. 

Therefore {e{x,y,s,i)} € y — group and we can al¬ 
ways redefine J such that J(x,y, s,i) = loi ® (—Id 2 )j 


y{x,y, s,i). From now on we will work within this gauge 
and denote the matrix I_Di ® (—lua) simply as J. 

This allows us to denote the rj{x,y,s,i) transforma¬ 
tion in Eq.(28) simply as y since it is site and virtual 


leg independent. In addition, according to Eq.(33), the 


remaining U-ambiguity: V (x, y, s, i) must commute with 
J. In other words, V{x,y, s,i) are block diagonal with 
two blocks, and the sizes of blocks are Di and £>2 respec¬ 
tively. 

Now we can consider an arbitrary symmetry transfor¬ 
mation R, which could be either a lattice symmetry or 
an on-site symmetry. Eq. (37) still holds for R and the y 


on the RHS must be J. Consequently we have: 


Wj^^{R{x,y,s,i)) 

=x{x,y,s,i)I. 


J • WR{R{x,y,s,i)) 


(38) 


Squaring this equation leads to y(x,y, s,i) = ±1. How¬ 
ever only the -|- sign is possible since otherwise the ma¬ 
trix Wji{R{x, y, s, i)) will not have a full rank, again due 
to Di ^ D 2 . Thus we have proved that WR{x,y,s,i) 
commutes with J, V(x,y, s,i) andWR. Mathematically, 
this means that when we extend the symmetry group by 
IGG = IGG X y — group, IGG is in the center of the 
extended group. 

Let us consider the phase factors yj(x, y, s) on site 
tensors obtained when applying the nontrivial element 
J on the virtual legs. This determines whether the 
site tensor is Z 2 even or Z 2 odd. Now we are ready 
to show that fj,j{x,y,s) is site independent in the cur¬ 
rent gauge. Namely if one site tensor is Z 2 even (odd), 
the same is true for all site tensors. Consider a lat¬ 
tice symmetry R which send a site {x,y,s) to the site 
(x', y', s'), Eq.( 131 states that the two site tensors are re¬ 


lated by a possible permutation of virtual indices (e.g. 
induced by a lattice rotation) together with multiplica¬ 
tions of Wr matrices on the virtual legs as well as a 
overall phase factor 0ij(x,y, s). Because Wr matrices 
all commute with J, it is straightforward to see that the 
yj(x, y, s) = yj(x',y', s'). Because all sites are related to 
each other by lattice symmetries, pj{x, y, s) are identical 
for all sites. Thus in the discussion below we will simply 
denote the y S IGG associated phase factors p{x, y, s) in 
Eq.(29) as /i, since it does not depend on the site. 


By applying the condition IGG = Z 2 to the kagome 
lattice with the symmetry group described in Appendix 
[A] we are able to solve the equations for symmetry oper¬ 
ations, i.e. Eq.( 28|29 ), by gauge fixing. For the purpose 
of presentation, here we only demonstrate the calculation 
for the translation symmetry, and list the full results of 
the classification. The calculation for other symmetries 
is in Appendix (We encourage interested readers to 
follow a simpler and pedagogical example on the square 
lattice in Appendix [F| before reading this more technical 
Appendix for the kagome lattice.) 

Let us consider the translation symmetry group. This 
group is isomorphic to Z x Z: the group is defined by its 
generators Ti, T 2 as well as the relation between genera- 
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tors, 


T-^t-^T 2 Ti = e 


(39) 


As shown in Eq.(13), for PEPS symmetric under {i = 
1, 2), we have 


V{s,i) and ^{x,y,s) = ‘h(s); namely they are only de¬ 
pendent on the sublattice index and the virtual leg index 
from a site, but are independent of the unit cell coor¬ 
dinates. Eurther, in this gauge, site tensors are transla¬ 
tional invariant (but could be sublattice dependent), 


rp(x,V,s) ^ o T<^x,y,s) 

^(xysi\x'y's'i') ® ^(xysi\x'y's'i') (^^) 

From the group relation T2^T^^T2Ti = e, we have 

WV/(r2(a;, y, s, i))W^^^{TiT2ix, y, s, i))WT2 (7iT2(x, y, s, i)) 
WT^{Ti{x,y,s,i)) = 'ni 2 Xi 2 {x,y,s,i) (41) 

as well as 


0 T 2 {T2{x, y, s))0r^ (TiT 2 (x, y, s)) 0 t 2 {TiT2{x, y, s)) 
0Ti(ri(a;,y,s)) = y,i2Wxi2{x, y,s,i) (42) 


where 7712 G {I, J}, and {xi 2 (a;, y, s, i)} G y - group. 

Under transformations Wri —t £TiWTi and 07 ^ —>■ 
STi ^Ti, we have 


X 12 ^ £t 2 (x, y + 1, s, i)er^ (x + l,y + l, s, i)- 

£t2(x + l,y + 1 , s, i)£T, (x -f 1, y, s, z)xi 2 (a;, y, s, i) (43) 


Thus, we are able to set all Xi2ix, y, s,i) = 1 via the 
ffTi-ambiguity. 

According to Eq.(32) and Eq.([M|, by doing a gauge 
transformation V(x,y, s,i) and multiply phase factors 
$(x,y,s): 

WT2{x,y,s,i) ^ V{x,y,s,i)WT2{x,y,s,i)V~^{x,y - l,s,z) 
0T2(a;,y,s) 0T2(a;,y,s)$(a;,y,s)$*(a;,y- l,s) 

(44) 


We are able to set WT 2 {x,y, s,i) = I as well as 
QT2{x,y,s,i) = 1. Thus we obtain T^x,y,s) _ rp{Q,y,s)^ 
The remaining U-ambiguity preserving the form of Wt '2 
should satisfy V (x, y, s, z) = U(x, 0, s, z), and the remain¬ 
ing ^-ambiguity preserving the form of 0 t 2 should sat¬ 
isfy <i>(x,y, s) = $(x, 0, s). In addition, any nontrivial 
6 X 2 transformation will change the form of Wt 2 = I, 
so £t 2 is fixed to be 1. Together with the condition 
Xi 2 (a^) y, S) z) = Ij ths remaining -ambiguity satisfies 
£T 2 (x, y, s, z) = £Ti (x, 0, s, z). 

Similarly, for Ti transformation, using the remaining 
U-ambiguity and <I>-ambiguity, we have 

WTi(x,y, s,z) U(x, 0, s,z)WTi(x,y, s,z)U“^(x - 1,0, s,z) 
0Ti (x, y, s) ^ 0Ti (x, y, s)$(x, 0, s)$*(x - 1,0, s) 

(45) 


Thus we can set Wti(x, 0,s,z) = I and 0 Ti(x,O, s) = 1. 
To maintain this form of Wr^, we find that there is no re¬ 
maining -ambiguity: £ti is fixed to be 1. The remain¬ 
ing U-ambiguity and ^-ambiguity satisfy V{x,y, s,i) = 


rf{x,y,s) ^J,(x,0,s) ^^{0,0,3)^ S = M,U,l(; (46) 


Thus, in the gauge that we choose so far, we can solve 
Eq.(41), and get the implementation of translation sym¬ 
metry on PEPS as 


Wti (x, y, s, z) = 77^2 
Wt 2 (a:, y, s, z) = I 
0Ti {x, y, s) = p\2 

0T2 (a:, y, s) = 1 (47) 


So for systems with translational symmetries and 
IGG = ^ 2 , there are at least two distinct classes of wave- 
function. In the context of quantum spin liquids, these 
two classes are known as zero flux state and tt flux state, 
corresponding to 7712 = I and 7712 = J respectively. Con¬ 
densations of spinons in these two spin liquids lead to 
different types of magnetic order^^. In the above gauge, 
although all site tensors related by the translation sym¬ 
metry share the same form, bond states related by the 
translation symmetry are in general different if 7712 is 
nontrivial. 

The calculation for other symmetries is similar as the 
above procedure. The basic idea is to keep fixing gauge 
by the four ambiguities. And when we find certain alge¬ 
braic data, such as the 7712 introduced above, that cannot 
be removed by the ambiguities, they describe different 
symmetric PEPS classes. We only list the result here, 
and put details in Appendix [B| 

This classification scheme will always lead to three fi¬ 
nite sets of algebraic indices 77 ’s, x’s and 0 ’s and we will 
discuss their physical meanings in Sec |IV| Although in 
general systems every set of indices is nonempty, for a 
half-integer spin system on the kagome lattice described 
by PEPS with IGG = Z 2 , we have: 


• V 12 , VCe Vcn where 77 G {I, J}. The correspond¬ 
ing / 7 i 2 ,yceIMcr are determined by 77 ’s. 


• X-T and XT, where x = ±1- 


• There turns out to be no tunable 0 indices in this 
example. 


So the number of classes equals to 2^ =32. By choosing 
a gauge, the symmetry operations on PEPS can be solved 
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as 


(a;, ?/, s, i) = 

WT2{x,y,s,i) = I, 

Tir f -N xy+^x{x+l)+x+y , 

Wce{x,y,U,l} = ^^2 WCe{u,l), 

TTT f ■\ xy+^x(x+l)+x+y 

W^Ce (a;, 2 /,-y, i) = ? 7 i 2 

Wce.{x,y,w,i) = ^“ 1+2 “(“+!) ^ 

Wa{x,y,s,i) = ? 7 i 2 ^^+^"u><,(s,i), 
Wr{x,y,s,i) = wr{s,i), 

Wgfi{x,y,s,i) = 0(I«i 


(48) 


In this gauge all Wr matrices are unitary. The last 
equation is for the SU(2) spin rotation along ft direc¬ 
tion by an angle 9. In addition, in this gauge we choose 
J = W 2 Tv{x,y,s,i) = (8) namely J is the 

direct sum of for the integer spin subspace and — 1^2 
for the half-integer spin subspace and Di + D 2 = D. 

For the rotation transformation wc^{u,i), we have 


wcs{u,a) = wcs{u,c) = I, 

WCe {u, b) = WCg {u, d) = T]i2VCe , (49) 

For the reflection transformation Wa-{s,i), we have 


w„{u,a) = I, 

Wa{u,c) = XcrVl 2 VCeVa, 
Wa{v,a) = ?7i2, 

Wa{v,c) = yCeVa, 
Wa{w,a) = XcrVCe, 
Wa{w,c) = yi 2 Va, 


W^{u,b) = X^ll2VCni 
Wa{u,d) = rja] 

Wcr{v,b) = X<jV12, 

Wcr{v,d) = XoVC^iVa] 
Wcr{w,b) = rjCe^ 

W^{w,d) = XcTVl 2 Va', (50) 


And for the time reversal transformation w-j-, we have 


W'f{u, a) = wq-, 

Wt{u,c) = m2VCeV<yWr, 
wr{v,a) = rii2'nCe.WT, 
wr{v,c) = T]aWr, 

W']-{w, a) = Wq-, 
wq-{w, c) = m 2 VCeVaWr, 


Wq-{u,b) = rji2r]CsWT, 
wr{u,d) = rjawr; 
wt{v, b) = wr, 

Wq-{v, d) = r]l2VC6V<yWr-, 
wq-{w,b) = r]l2VCeWT, 

wriw,d) = rj^wr; (51) 


where 


wj- 


0,(I«i 8 6^’"'^?) ifx-7- = l 

0* i^ni 8 ) if XT- = -1 


(52) 


Here is dimension of the extra degeneracy associated 
with spin-S'i. Namely, the total degeneracy for spin-S^ 
living on one virtual leg equals rii x {2Si -1-1). We have 
the virtual bond dimension 


D = Y,n^{2S, + l) (53) 

i 


And, = icTy 8 lrai /2 is a rii dimensional antisymmetric 
matrix. 

For &rS, we have 


0Ti(a;,y,s) = 2^12, 

<dT 2 {x,y,s) = 1 , 

fc>C 6 (a;,y,u) =/ri 2 Bce(u), 

„ / X xy+\x(x+l)+x+y 

0 C 6 (a:, 2 /,y) = Mi 2 

r\ f \ xy+^x{x+l) 

ec^{x,y,w) = 


e^{x, y, s) = /Zi^^+“*' 0 <,(s), 
QTix,y,u/w) = 1, 
0 r(a;, 2 /,y) = f^i 2 yce^ 

06 (n = Ij 


(54) 


where 


0 C 6 (w) = {fJ- 12 tJ-CeV^-, 

0 CT(y) (/^(t)^5 

0^(y) = y,CeQceiu)Qaiu); 

Qaiw) = fJ.^HCe{QCeiu)&aiu))~^. (55) 


Note that in Eq.(55) 0Ce(u) and 0 ct(u) contain square 
roots so there appear to be two possible values of each of 
them differing by a minus sign, giving rise to 0 -indices. 
However, these minus signs can be tuned away using the 
T^-ambiguities in the definition of Wc^ and Wg- since ev¬ 
ery site tensor is Z 2 odd. So one could simply fix an 
arbitrary choice for the square roots here. This is the 
reason why there turns out to be no tunable 0 indices in 
this example. 

Even after all these transformation rules are deter¬ 
mined by gauge fixing, we still have some remaining V- 
ambiguity for each class. (Note that there is no remain¬ 
ing nontrivial ry,£ and 4> ambiguities.) To preserve the 
lattice symmetry, the remaining E-ambiguity is indepen¬ 
dent of sites and legs. To preserve the form of Won, the 
remaining V-ambiguity must have the following form: 


4^ = 0(V5. 8 l 2 Si-i-i), (56) 

i 


where Vs^ is a rii dimensional matrix. In addition, the 
time-reversal transformation JT 7 - further constrains the 
form of component matrices hs.. When xt = 1) on® can 
show that Vsi must be a real matrix. For the purpose of 
presentation we only consider XT = 1 classes here. The 
XT = ~ 1 cases involve quaternion matrices and we leave 
the general and detailed discussions in Appendix [B| 
Next, we are at the stage to construct the con¬ 
strained sub-Hilbert spaces for building block tensors for 
all classes, according to the Wqi transformation rules. 
The basic idea is to determine the generic form of a 
single site/bond tensor using the WA’s with R leav¬ 
ing the site/bond invariant, and then generate all other 
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site/bond tensors using all Wrs. The generic forms of 
site tensors are straightforwardly determined in this fash¬ 
ion, with a set of real continuous variational parame¬ 
ters whose number basically equals the dimension of the 
constrained site sub-Hilbert space. However, for bond 
tensors, we will use the remaining H-ambiguity to bring 
them into canonical forms which are maximal entangled 
bond states containing no continuous variational param¬ 
eters. 

To make sure a bond tensor Bi, to be invariant under 
the SU{2) spin rotation, it must have the following form: 


M 


Bh = 


®( 

i=l 


Bt 


•Ks, 


(57) 


where B^' is Hi dimensional matrix, and Ks^ is the fixed 
(25'i-l-l) dimensional matrix representing the spin singlet 
formed by two spin-S'i on the two virtual legs shared by 
Bf,. For example, we get Ks=o = 1, ^s=i = io'y 

As shown in Appendix]^ when XT = 1 3 'iid a given Si, 
depending on the four possible values of rja- and Xa-, the 
component matrix B^' must be a purely real/imaginary 
symmetric/antisymmetric matrix. Then we can use the 


remaining H-ambiguity in Eq.(56) to simplify B?\ be¬ 


cause under a Vs- transformation, B^' transforms as: 


Bt 


VsrB^^-Vl 


(58) 


Clearly we can use a real orthogonal Vg. to diagonalize 
(block diagonalize) B^' if B^' is a symmetric (antisym¬ 
metric) matrix. After this, the eigenvalues of B^' could 
have arbitrary norms. But then we can use another real 
diagonal V 5 . matrix to normalize the eigenvalues so that 
they are only ±1 (if B^' is purely real) or ±i (if B^' is 
purely imaginary). 

This procedure fixes Bf, to be maximal entangled states 
with no continuous variational parameters. However, the 
relative number of -l-l(-l-i) eigenvalues and — 1 (—i) eigen¬ 
values cannot be further tuned away by gauge fixing and 
will serve as discrete variational parameters on the bond 
tensors. 

The previous discussions in the subsection are general 
for any half-integer spin-B. Below we focus on the ease 
with S = ^. For simplicity, we demonstrate the results 
for with D = 3. The basis of virtual legs of site tensors 
are {| 0 ),| f)) I i)}- Namely, virtual legs are formed by 
one spin singlet and one spin doublet. Note that virtual 
legs of bond tensors are dual to those of site tensors, so 
the basis are (0|, (f |, (i |. Symmetric PEPS with larger 
D are also conceptually straightforward but technically 
involved to obtain, and we leave the general construction 
in Appendix [B] 

As discussed in Appendix only classes satisfying 
T^cr = J, Xo- = 1 and XT = 1 can be realized with D = 3. 
So the realizable classes reduce to 2^ = 4 with D — 3. 
At such a small D, it turns out that each class has only 
two continuous variational parameters. (Note that for 


B = 6 , i.e. two spin singlet and two spin doublet on the 
virtual leg, we find that all the 32 classes can be realized. 
And each class has 47 continuous variational parame¬ 
ters.) Following the above procedure we can bring the 
bond tensor on a given bond bo into the canonical form: 

/±1 0 0 \ 

Bf,„ = 0 0 -i (59) 

VO i oy 

All other bond tensors are generated by combination of 
translation and rotation symmetries as: 

BrQj) = R ^WrRoBi,^ (60) 

where R = with ni,n 2 , ncg G Z. 

One can view a bond tensor as a quantum state living 
in the Hilbert space formed by the tensor product of two 
virtual legs. Namely, we have 

Bf,„=±(0,0|-i(t,;|+i(i,tl (61) 


Here we use notation B),„ as the quantum state represen¬ 
tation while Bhg as the matrix (tensor) representation. 

At a given site Sg, the generic form of the site tensor 
for all classes can be summarized as: 


={Ko + ki2{pi,P2)} + QceikW ^ b,C ^ d} + Qaiu)- 
{a -(rT d,b -(tT c} + Aii 2 MC 6 0C6(w)0<T('u){a c, 6 o d} 

(62) 

with real continuous parameters pi,P 2 - Here a, b, c, d de¬ 
note virtual leg of sites, as shown in Fig.(|^. Kq and ki 2 
denote linear independent spin singlet states, which can 
be expressed as 

iFo =1 t) ® U 000) - U) 0 I t 000) 
ki 2 =Pi - (It)® |o m) + u) ® |o nt))+ 
P2-(it)®iout) + u)® iom))- 
{pi +P 2 )- (It)®|om) + U)®|o;tt)), (63) 


where the first spin lives on the physical leg, while the 
following four spins live on virtual legs a, b, c, d respec¬ 
tively. Note that we have chosen a particular gauge such 
that all site tensors share the same form. 

By direct comparison, the NN RVB state (Qi = Q 2 
state) given in Sec.(HE) is represented as the PEPS de¬ 
fined in Eq.(59) and Eq.(62), with pi = p 2 = 0 and: 


m 2 = VCe = I. Va = J; 

X<y=XT = 1; 

(64) 


C. Algorithm for minimization 

As demonstrated above, after the symmetry trans¬ 
formation rules are determined for each class, for a 












19 


fixed bond dimension D together with the specified on¬ 
site symmetry (projective) representation on this virtual 
Hilbert space, we can construct the generic symmetric 
PEPS. The strategy is to firstly construct one site ten¬ 
sor and one bond tensor, and to use spatial symmetries 
to generate all site tensors and bond tensors. For each 
class, in general, the maximal entangled bond tensor will 
be completely fixed up to a finite number of ±1 signs. 
These signs are physical (they could modify the energy 
density) and should be treated as discrete variational pa¬ 
rameters. The site tensor, however, will be fixed in a 
sub-Hilbert space for each give class, whose dimension 
(minus one if one considers normalized PEPS) represents 
the number of continuous variational parameters. Differ¬ 
ent classes give different sub-Hilbert spaces for the site 
tensors. 

For efficient PEPS simulations, one generally chooses 
an open-boundary sample and introduces certain bound¬ 
ary conditions. In order to determine which symmetric 
PEPS class that the ground state of a model Hamiltonian 
belongs to, one needs to minimize the energy density near 
the center of the sample by tuning all the variational pa¬ 
rameters of each given class. And the ground state crude 
class is identified as the class which gives the lowest en¬ 
ergy density. The effect of boundary conditions should 
decay quickly into the bulk of sample and the optimal en¬ 
ergy density near sample center is expected to converge 
to a boundary condition independent value for interme¬ 
diate sample sizes. 

By construction, PEPS belonging to different classes 
has distinct short-range physics (the site tensors and 
bond tensors transform according to inequivalent alge¬ 
braic equations), therefore we expect that quite generi- 
cally these optimal energy densities for different classes 
are significantly different. But, there are three scenar¬ 
ios that the optimal energy densities for distinct classes 
can be identical, at least in the thermodynamic limit. 
The first scenario is that the quantum Hamiltonian un¬ 
der investigation is right at a phase transition point, so 
that the optimal energy densities for the two symmetric 
PEPS on the two sides of the phase transition are the 
same. This scenario, however, is expect to occur only at 
phase transition points. 

The second scenario has a deeper physical origin — it 
is possible that the invariant gauge group of the PEPS 
representation of a quantum phase is larger than the as¬ 
sumed IGG for the symmetric PEPS classification. We 
will comment more on this scenario at the end of the 
paper. For instance, in the example of kagome spin- 
1/2 systems, we classify all the symmetric PEPS whose 
IGG = Z 2 in the previous section. But it is possible that 
a given quantum phase actually require IGG larger than 
Z 2 . A PEPS with a larger IGG indicates more invari¬ 
ance of the tensor network, and therefore it could be well 
approximated (with arbitrarily small wavefunction dif¬ 
ferences) by PEPS with the minimal required Z 2 IGG. 
Physically this is related to the Higgs mechanism break¬ 
ing the large IGG down to ^ 2 - In particular, it is possible 
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FIG. 5. Illustration of truncating an infinite PEPS to a finite 
size PEPS with a choice of boundary condition. 


that the same PEPS with a larger IGG can be well ap¬ 
proximated by two or more distinct IGG = Z 2 PEPS 
classes. This picture can actually be used as a feature of 
the algorithm proposed here. Namely, if numerically one 
finds that two distinct classes of IGG = Z 2 PEPS give 
the same optimal energy densities over a finite region in 
the phase diagram, it could be due to the fact that the 
true IGG of the quantum phase is larger than Z 2 . 

The third scenario is the one that we currently do not 
fully understand. It seems possible that certain symme¬ 
try breaking conventional phases such as a valence bond 
solid phase that can be represented by distinct classes 
even with the minimal required IGG. It would be very 
interesting to perform numerically simulations on quan¬ 
tum models to understand whether this really occurs or 
not, and if this scenario occurs, how does it occur. We 
will comment further on this issue later in Sec IVIl 

Note that although the bond tensors are position de¬ 
pendent in our construction, it is straightforward to ab¬ 
sorb each bond tensor to a neighboring site tensor so 
that all the bond tensors become identity matrices. After 
this transformation the PEPS is in the conventional form 
used in numerical simulations. As a result, for a PEPS 
with given variational parameters, the measurement of 
the energy density near the sample center, which is one 
or few terms in the translational symmetric Hamiltonian, 
can be carried o ut in the same way as in the original 
PEPS algorithmd^l^. The only new ingredient of the 
current minimization scheme is an algorithm to perform 
the energy density minimization by tuning the continu¬ 
ous variational parameters within each class. Below we 
propose such an algorithm based on the conjugate gradi¬ 
ent method. Similar minimization algorithm was investi¬ 
gated in the context of one-dimensional matrix product 
states with translational symmetrjP^. 

Before moving towards the minimization algorithm, 
we introduce the well-established numerical me thod s in 
PEPS to perform energy density measurementd^*^. As 
a demonstration of principles, let us consider a spin- 
1/2 symmetric PEPS on the square lattice; namely, for 
any site (bond) on the infinite square lattice, we have a 
well-defined site (bond) tensor. We firstly cut the infi¬ 
nite square lattice to a La, x Ly finite sample. And the 
open boundary condition around the edge of the sample 
is specified by a certain choice of the virtual quantum 
states living on the dangling bonds. For instance, one 
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FIG. 6. Illustration of measuring wavefunction norm and en¬ 
ergy density of a PEPS. 


can choose a simple boundary condition by requiring all 
the virtual states on the dangling bonds to be the spin 
singlet |0) in virtual Hilbert space (see Figj^for illustra¬ 
tion). 

With a well-defined finite size PEPS, the computa¬ 
tion of wavefunction norm and the energy density near 
the center of the sample is illustrated in Figj^ Note 
that the total Hamiltonian H is a summation of many 
terms which are related to each other by lattice symme¬ 
try. For instance one can consider H = ‘ 

Ss> — J2<xx'> ^xx', where < xx' > represents the near¬ 
est neighbors and each energy density term is given by 
= JSg ■ Sg'. Because by construction the PEPS is 
lattice space group symmetric, indicating that as system 
size increases, expectation values of all the energy den¬ 
sity terms near the center of the sample would quickly 
converge to the same value. Therefore it is legitimate to 
optimize only a single energy density term as shown in 
FiglU 

These measurements become the standard problem in 
PEPS — contracting a finite-size two-dimensional tensor 
network with the bond dimension D^. Here the bond di¬ 
mension is squared because we view a pair of a site tensor 
with its conjugate as a single tensor as shown in Figj^ In 



T*{x) 


FIG. 7. Each site tensor in a symmetric PEPS is a linear 
superposition of the states in a sub-Hilbert space. 


general there is no way to exactly contract such a tensor 
network as long as the system size is not very small. How¬ 
ever PEP^ and other tensor network method^SSHUl pro¬ 
vide efficient algorithms to approximately contract such 
tensor networks to high accuracies. The basic idea of the 
PEPS contraction algorithm is to view the first row of the 
2d tensor network as a matrix product state (MPS), and 
view the other rows as matrix product operators acting 
on the MPS. Consequently the existing algorithm of com¬ 
pressing a MPS allows one to approximately contract the 
whole 2d tensor network in a row-by-row fashion. The 
details of the PEPS contraction algorithm can be found 
in review articles such as Ref|75| and |H31 

Now we present an algorithm to perform the energy 
density minimization, based on the well-established con¬ 
jugate gradient minimization algorithm. Namely, the 
quantity we need to minimize is: 




{i^\hxx'\'fp) 

(V’lv’) 


(65) 


where {pi} represents the finite collection of continuous 
variational parameters in a given symmetric PEPS class 
with a bond dimension D. When {pi} are tuned, every 
site tensor in the PEPS wavefunction \ij}) is modified. In 
order to apply the conjugate gradient minimization algo¬ 
rithm, one needs to compute the following derivatives: 

dE _ 1 d{tjj\hgs'\ip) {ip\hgff\'ip) dj-tpl-tp) 

dp^ {ip\ip) dpi dpi 


The only new quantities that we need to compute are 
^ and these quantities can also be 

efficiently computed using PEPS algorithms. 

Due to the symmetric PEPS construction, we know 
that each site tensor is a linear superposition of the states 
in the symmetry constrained sub-Hilbert space (see FiglT] 
for illustration): 


T(x)=^p,T,(x), (67) 

i 

where x labels the real space position of the site tensor. 
(Because the overall factor of all piS does not change 
the normalized wavefunction, it may be convenient to 
set one of the pi’s, say po to be unity and only study 
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structing symmetric PEPS classes still works in higher 
dimensions, the PEPS algorithm to contract the tensor 
no longer applies. However, it is still possible to approx¬ 
imately contract tensor networks in higher dimensions 
while giving up certain accuracy, efficiency and larger sys¬ 
tem sizes. For exa mple, the tensor renormalization group 
techniques (TRGp2Hll could be used to contract ten¬ 
sors in higher dimensions. It would be interesting to see 
whether the TRG algorithms, when combined with the 
symmetric PEPS construction studied here, can be used 
to efficiently study quantum phase diagrams in higher di¬ 
mensional systems. We leave this as a subject of future 
investigations. 



FIG. 8. The linear order variations of d{'ip\%l}) and {'ip\hggi\il)). 


the variation of other pi’s.) Note that for different site 
X, the form of Ti{x) are generally different due to the 
nontrivial symmetry transformation rules of tensors and 
the fact that we absorb the bond tensors into neighboring 
site tensors. 

It is then straightforward to see that to the linear or¬ 
der, the variation of (V’l^) (('01^x5'I'*/’)) is a summation 
of Lx X Ly X 2 terms, as illustrated in Fig|^ where the 
factor of two is due to the fact that both the variation of 
\tjj) and the variation of ('i/'j contribute. Each term can be 
computed by the standard PEPS contraction algorithm. 

Basically one needs to compute the derivatives 

each variational parameter pi. For 
each derivative one needs to compute ^ x Ly terms. 
(This number can be reduced by a factor if one chooses an 
open boundary sample respecting the point group sym¬ 
metry.) However, one notes that every derivative term 
only involves one modified site tensor. So one could pro¬ 
gram the contraction sequence so that all the other ten¬ 
sors are contracted first, leaving an “environment tensor” 
E{x) for the only remaining site tensor at x, which is ex¬ 
actly the modified site tensor. Then, different pi deriva¬ 
tives on the given site x can be efficiently computed by 
contracting the different Ti{x) with the same E{x). We 
then end up with a factor of ~ L^xLy in the computation 
complexity due to the choices of modified sites. However 
this factor ^ L^ x Ly in the computation cost can be 
straightforwardly parallelized on ^ x Ly computer 
nodes, simply because no communication is needed be¬ 
tween different contractions. Gonsequently with a com¬ 
puting cluster, the minimization algorithm proposed here 
can be efficiently implemented. 

Finally we comment on possible algorithms in spatial 
dimensions higher than two. Although our scheme of con- 


IV. PHYSICAL INTERPRETATION OE 
CLASSES 

We will discuss the physical meanings of different 
classes, which are labeled by 0 _r, xr well as pr. In 
this section, we will focus on the non-symmetry-breaking 
liquid member phase in each crude class. We will com¬ 
ment on the meaning of these indices in the long-range 
ordered member phases in Sec|VI[ 


A. Interpretation of Qr and xr 


Although it happens to be true that the kagome half¬ 
integer spin example has no tunable Qr indices. Or in¬ 
dices do appear in general quantum systems. In Ap¬ 
pendix]^ we perform the crude classification for the half¬ 
integer spin systems on the square lattice with a space 
group generated by translation symmetries and the C 4 
rotation symmetry only. Assuming IGG being the min¬ 
imal required Z 2 , we find a tunable 0 C 4 index. 

In fact, the Or indices and the XR indices generally 
appear even when the IGG is trivial. For instance, we 
could consider a system on the kagome lattice with no 
on-site symmetry (i.e., remove the spin SU(2) rotation 
and the time-reversal symmetry in our main example), 
and consequently the minimal required IGG is trivial. 
Assuming IGG being trivial in this system, we will not 
have the p indices but still have the x indices. The calcu¬ 
lation procedure of transformation rules almost remains 
the same as before if we simply limit all the 77 ’s to be 
identity. Eventually we will arrive at Eq.(55l replacing 


all the pr by -|-1. Note that there is no Ty-ambiguities 
to tune away the signs for the square roots as in the 
half-integer spin case. In this system, apart from the x 
indices, we do have two tunable 0 indices in the PEPS 
classification: 0c6(w) = ±I and 0 ^{u) = ±1. 

Different Or indices can be viewed as different sym¬ 
metry quantum numbers (for either on-site symmetries 
or space group symmetries) carried by each site tensor. 
These quantum numbers of the site tensors, generally 
speaking, directly contribute to the quantum numbers 
of a finite size sample. The physics of Or indices is 
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similar to the physics of the so-called “fragile Mott in¬ 
sulator” discussed by Yao and KivelsorP^l. And similar 
indices in one-dimensional matrix product states have 
been investigated recently^. For instance, in the fragile 
Mott insulator example^, a Mott insulator wavefunction 
is constructed on the checkerboard lattice which carries 
nontrivial point group quantum numbers on the odd- 
by-odd unit cell lattices. This distinguishes the fragile 
Mott insulator from trivial insulators which carries triv¬ 
ial quantum numbers on the same lattices. And such 
nontrivial quantum numbers can be traced back to the 
quantum numbers carried by the wavefunction on every 
square cluster on the checkerboard lattice. If one tries to 
use a site tensor in PEPS to represent the square clus¬ 
ter wavefunction, it is clear that this site tensor forms a 
nontrivial representation of the point group symmetry. 

The physical meaning of XR may be more well-known. 
These are generalizations of the symmetry fractionaliza- 
tions in the 2d AKLT modeP^. Let’s hrstly briefly de¬ 
scribe the PEPS construction of the SO{3) symmetric 
spin-2 AKLT state on the square lattice. In this con¬ 
struction, each virtual leg forms a spin- 1/2 projective 
representation of the SO{3) symmetry group of the spin- 
2 system. Each site tensor is given by the only singlet 
state formed by the physical spin -2 and the four vir¬ 
tual spin-l/ 2 ’s, and each bond tensor is formed by the 
only spin singlet formed by the two spin-l/ 2 ’s on the 
two ends of the bond. Such an AKLT wavefunction can 
be shown to be the unique gapped ground state of the 
AKLT Hamiltonian on the square lattice with periodic 
boundary condition^^. 

However, when the system has an open boundary, one 
needs to specify a symmetric boundary condition. But 
one encounters the following problem: each site tensor on 
the boundary has only three virtual spin- 1 / 2 ’s and it is 
impossible for form a spin-singlet with the physical spin- 
2. Basically each site on the boundary can be viewed as a 
half-integer spin — which is a projective representation 
of the original group. One sometimes calls this 

phenomena as the symmetry fractionalization in 2 d in 
the absence of topological orders. When coupled together 
along a translational symmetric edge, the low energy dy¬ 
namics of the edge states can be effectively described by 
a translational symmetric half-integer spin chain, which 
would give a gapless excitation spectrum assuming no 
spontaneous translational symmetry breaking. Clearly, 
in the PEPS construction, the origin of such symmetry 
fractionalization behavior is due to the fact that projec¬ 
tive representations appear in the virtual legs. 

For an on-site symmetry R, this is exactly the physics 
that xr captures. For instance, the XT index appear¬ 
ing in the kagome example is really about the projective 
representations of the symmetry group SU(2) x T on 
the virtual legs. As mentioned before, when XT = Ij 
the half-integer (integer) spins on the virtual legs form 
Kramer doublet (singlet) under the time-reversal trans¬ 
formation. This is the usual representation of SU ( 2 ) x T. 
However when XT = ~lj the half-integer (integer) spins 


on the virtual legs form Kramer singlet (doublet) un¬ 
der the time-reversal transformation. This is a nontriv¬ 
ial projective representation of SU{2) x T. We expect 
that XT = would give rise to nontrivial signatures in 
entanglement spectra and physical edge states. 

For a spatial symmetry i?, the physical meaning of XR 
is less obvious. But it’s one-dimensional analog has been 
investigated in the context of matrix product state^^^J^^. 
In our example, the x<t is capturing similar physics in 2d 
kagome lattice, which basically describes how the tensor 
network forms possible projective representations of the 
spatial reflection. We speculate that nontrivial X(t would 
give rise to signatures in entanglement spectra when the 
partition of the system respects the a reflection. 

In summary, Qn is capturing local contributions to 
symmetry group quantum numbers, and XR is capturing 
the symmetry fractionalizations not due to topological 
orders. 


B. ga and symmetry fractionalization 

Here, we will show that rj^s are directly related to the 
symmetry fractionalization of spinon excitations (char- 
gons). To see this, let us firstly introduce the concept of 
symmetry fractionalization in the presence of topological 
orders. We will use the unitary on-site symmetry as an 
example. Related discussions can be found in RefHOl and 
ReflSIl 

Starting from a topologically ordered ground state 
with a global symmetry group SG, consider an ex¬ 
cited state, having n—quasiparticles (which do not have 
to be of the same type) spatially located at position 
ri,r 2 ,... ,r„, far apart from one another. Let’s denote 
this state by |'0(ri, r 2 ,..., r„)). For any symmetry trans¬ 
formation U{g) by a group element g G SG, U{g) will 
generally transform this state to another state: 

C/(g)o|' 0 (ri,r 2 ,...,r„)) |• 0 (rl, r 2 ,..., r„)) ( 68 ) 

One way to describe the symmetry fractionalization on 
quasiparticles is the following condition: there exist lo¬ 
cal operators Ui{g), 1 / 2 ( 5 ), ■ • ■: Un{g), such that Ui{g) is 
a local operator acting only in a finite region around the 
spatial position r^, and does not touch the other quasi¬ 
particles; in addition, Ui{g), 1 / 2 ( 5 ),..., Un{g) satisfy: 


Ui{g) ■ 1/2(5) ■ • • Un{g)\tp{ri,r2, ... ,r„)) 
=[/(5)|V’(ri,r2,...,r„)) = |^(ri, rz,..., r„)) ( 69 ) 


Pictorially, this condition is shown in Fig.(|^. 

Note that technically Eq.(69l is not a general condi¬ 
tion for symmetry fractionalization phenomena. For ex¬ 
ample, let us consider SG to be an on-site U{1) symme¬ 
try, and assume that Eq.(69l holds for a wavefunction 
|'!/'(i'i, r 2 , ■ • ■, r„)). We can then just add one extra 1/(1) 
charge outside the regions that Ui[g) {i = l,..,n) act 
and obtains a new wavefunction |'!/'(ri, r 2 ,..., fn))- It is 
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FIG. 9. Illustration of symmetry fractionalization phenom¬ 
ena: Under a symmetry transformation U{g) with 'ig £ SG, 
an excited state is transformed by the product of local trans¬ 
formation operators Ui{g), with each operator only acting on 
one quasiparticle locally. 


perfectly fine to imagine the extra charge as if it already 
exists in the ground state. Physically the local operators 
that transform quasiparticles: Ui{g) for \ijj) should be 
exactly the same as before, since |V’) and |V’) are locally 
identical around ri, r 2 ,..., r„. However, clearly Eq.(69) 
is no longer true for \ijj), because the global symmetry 
U{g) picks up an extra U{1) phase from the added 17(1) 
charge. 

In fact, Eq.(69) implicitly assumes that, under a global 
symmetry transformation, there is no phase “locally ac¬ 
cumulated” in the ground state wavefunction. But, 
as demonstrated above, generally there could be such 
“locally accumulated” phases in the ground state, and 
Eq.(69l should be modified up to the “locally accumu¬ 
lated” phases outside the blue regions in Eig[^ 

How to sharply define such “locally accumulated” 
phases in general ? The answer to this question is impor¬ 
tant to provide a general sharp definition of Ui{g). But 
to answer this question, one needs a tool capable to di¬ 
agnose wavefunctions locally, which is exactly the power 
of PEPS. For the moment, let us postpone answering 
this question in the framework of PEPS, and have some 
further discussion on symmetry fractionalizations. 


First, fractionalized symmetry transformations are 
local operators and cannot change the quasiparticle’s 
species (or more precisely, the superselection sector of 
a quasiparticle). Thus, we can investigate the transfor¬ 
mation rules of each anyon species individually. However, 
anyons do not need to form a representation of 5'G' due 
to the nontrivial fusion rule. For example, in a Z 2 topo¬ 
logical ordered phase, two chargons fuse to one trivial 
particle. We can multiply each chargon in the system by 
a fixed element in an IGG' = Z 2 = {1,-1}. Clearly, the 
total phase becomes unity, and physical wavefunction is 
invariant. Here IGG' is the subgroup of U(I) describ¬ 
ing the fusion rule of chargons. Quite generally for a 
topological order, IGG' = Z^- 

A PEPS with IGG = can describe a deconfined 


phase with a Z^ topological order. We will only con¬ 
sider this case and we do have IGG' = IGG. So IGG 
tells us that when we implement the global symmetry 
transformation on chargons, it is perfect fine to have a 
phase ambiguity, if this phase ambiguity is an element in 
IGG. Consequently, a single quasiparticle could form a 
projective representation of SG with coefficient in IGG, 
which is classified by second cohomology H‘^{SG,IGG) 
(see Appendixfor details). 

Now, let us translate the above discussion into the 
PEPS language. The main task is to construct the lo¬ 
cal symmetry transformation operators for a small patch 
of PEPS with a nontrivial IGG. Here we focus on 
IGG = Z 2 case. Without loss of generality, we assume 
that tensors of the PEPS are all Z 2 even. Then we cut 
a small patch A from the PEPS. We can view the tensor 
associated with patch as a linear map from boundary 
virtual legs to physical legs living in the bulk of the patch, 
which is labeled as T^. Here 0 denotes that there is no 
quasiparticles inside A. Namely, 

f^ = Y.(TA)iv\I){V\ (70) 

I,v 

where |/) labels ket states of all physical legs inside A, 
while {V\ labels bra states of all boundary virtual legs. 

Before studying excitations inside A, we firstly discuss 
properties of T^. As a tensor, is Z 2 even. Namely, 
action of the nontrivial Z 2 element g on the boundary legs 
of leaves the tensor invariant. This property implies 
that r^, as a linear map, can never be injective. To see 
this, consider an arbitrary boundary state \V), we have 

Tl\V)=f'^,\goV) (71) 

So, the inverse map of is not well defined. To have 
a reasonable definition of the inverse map, one observes 
that an arbitrary boundary state \ V) can be rewritten as 

\v) = \m + \ 9 oV)) + \{\V) - \goV)) 

= nu\V) + [l-nu)\V) (72) 

where U is the Z 2 even sector of boundary legs. Namely, 
V|t7) G U, we have \g oV) = \V). Au is a projection 
operator which projects a boundary state into U. Under 
T^, the second term in the above equation is mapped to 
zero. For a generic PEPS with IGG = ^ 2 , we can further 
assume that is injective on the subspace lA when the 
patch A is not too small. This is because the dimension of 
the physical Hilbert space increases parametrically faster 
than the dimension of the boundary virtual Hilbert space 
as the patch size increases. Such a PEPS is named as a 
Z 2 injective PEPS in ReflSSl Namely, generically one 
can find a linear map (T^)“^ from bulk physical legs to 
boundary virtual legs, such that 

(fO)-i ■f\ = Au 


(73) 
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Next, let us study the case with topological excitations 
inside patch A. One could create odd number of chargons 
near the center of the patch A by modifying to some 
Z 2 odd tensor T^. Opposite to the previous case, we 
have 


niv} = o, y\v)&u 


(74) 


Generically we can further assume is injective on the 
Z 2 odd sector of boundary legs. Namely, one can con¬ 
struct as linear map from bulk legs to Z 2 odd 

sector of boundary legs, such that 


{TX)-^ • = % 


(75) 


where Ajj = 1 — 11 ;^. 

Similarly, one can construct patch tensors with even 
number chargons inside the patch by modifying to 
any other Z 2 even and Z 2 injective tensors. For example, 
let us assume T\ to be such a tensor. Then, one can find 
it inverse on the subspace U, such that (T^)“^ • 

f\=Au. 

In the following, we will study the local physical oper¬ 
ator acting on small patches for a symmetry R. Starting 
with a PEPS wavefunction with topological excita¬ 
tions inside small patches A,B^..., while the region out¬ 
side these patches share the same tensors as the ground 
state wavefunction Idto). The action of the symmetry R 
on 14') is ob tained b y acting R on all tensors, which is de¬ 
fined in Eq. ( | 6 | 10 | 12 D . Since we try to construct local sym¬ 
metry operators only on patches A,B,..., we can apply 
gauge transformations Wr on all virtual legs in the re¬ 
gion outside all small patches as well as on the boundaries 
of all small patches, but leave virtual legs inside small 
patches untouched. Note that this gauge transformation 
does not modify the i?-transformed physical wavefunc¬ 
tion at all. Because tensors outside small patches are the 
same as tensors of ground state, the following relations 
still hold for them: 


patch. And Wr in Eqj^is defined to only act on bound¬ 
ary virtual legs of T 4 . Note that is either Z 2 even 
or Z 2 odd, which corresponds to even number chargons 
or odd number chargons inside A. Note that we should 
always choose the patch that is large enough so that all 
quasiparticles exist in the patch before the transforma¬ 
tion keep staying in the patch after the transformation. 
The above equation can be viewed as the definition of 

In fact, Eq.((7^ is a very general result which is ap¬ 
plicable even when the condition of symmetry fraction- 
alizations breaks down. For example, it is possible that 
certain symmetry transformation interchanges quasipar¬ 
ticle superselection sectors. In the PEPS formulation this 
happens when and T 4 describes distinct quasiparti¬ 
cle species, and consequently there is no way to use a 
local physical operator in A to send TA to T^. For the 
kagome example this would never happen. For example, 
we showed that Wr matrices all commute with the non¬ 
trivial IGG element g = 3, and therefore the parity of 
the number of chargons would be the same in and 
T 4 . But in a symmetric PEPS with a larger IGG (e.g. 
IGG = Z 2 X Z 2 ), we expect that it is possible that Wr 
does not commute with a g G IGG. In this case the R 
may interchange quasiparticle species. 

Below we only consider the situation that and T 4 
support the same superselection sector and consequently 
share the same Z 2 parity. This allows us to construct the 
fractionalized local physical operator for the symme¬ 
try R acting on patch A that realizes Eq. ([77]); namely: 

Liofj, = WRRofA, (78) 

at least for those T 4 describing the relevant low energy 
states. One should keep in mind that only acts on 
physical legs, without touching boundary legs; i.e., 

Li = Y.^Li)RR\I){I'\. (79) 

ij' 


T^ = &rWrRoT^ 

B^ = WrRoB^ (76) 

Thus, under the symmetry R together with the gauge 
transformation Wr defined above, tensors outside 
patches will be invariant up to an “locally accumulated” 
phase risGOutside emphasize that this actu¬ 

ally provides the sharp definition of the “locally accumu¬ 
lated” phases mentioned earlier in this section. As dis¬ 
cussed in the previous subsection, 0i?(s)’s exactly cap¬ 
ture the local phases picked up after applying a global 
symmetry transformation. Without the tool of PEPS, it 
is actually difficult to sharply define this object. 

For tensors inside patches, we have 

f^ = WRRofA (77) 

Here, T 4 is the linear map associated with patch A, 
which is obtained by contraction of all tensors inside A 


To obtain the explicit form of this local operator, let us 
consider a particular tensor T^, which supports an odd 
number of chargons in A. We have 

fe^R ^ , ye (gQ) 


where = WrR o T^, and (T^) ^ is defined in 

Eq.(75). In the above equation we assume that both 
and is Z 2 odd as well as injective in the Z 2 odd 
subspace of boundary legs, which is expected to be gener¬ 
ically true. Note that ■ (T^)”^] can be viewed as 
an operator acting only on physical legs. 

To study the transformation rules for a number of char- 
gon excitations, let us consider a finite set A of tensors: 
A = {T^\i = 0,1,...} in the patch A. These ten¬ 
sors may describe states with chargon number equal to 
zero, one, two, etc, and are injective in the correspond¬ 
ing boundary Z 2 sectors respectively. But tensors in A 
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FIG. 10. (a): Tensor Ta and its “generalized inverse” 
associated with patch A. Ta is obtained by contracting all 
bond tensors and site tensors inside patch A. As a linear 
map from boundary legs to bulk legs, Ta is either Z 2 even 
or Z 2 odd. (b): The local i?-symmetry operator on patch 
A. is an orthonormal basis, where every state in the 

basis is either Z 2 even or Z 2 odd. (c): The local symmetry 
operator for a series symmetry operations Ri ... R„ , where 
Ri ... Rn = I. If rju is nontrivial, action of this operator on 
Z 2 even or Z 2 odd tensor gives different phase factor. This 
indicates symmetry fractionalization of chargons. 


contain no fluxon excitations in A. (we will study the 
symmetry fractionalization of fluxons later in this paper.) 
We assume that any symmetry transformation as shown 
in Eq. ([77]) transform within the linear space spanned by 

A. 

In addition, we assume the tensors in A to satisfy 
= 0, Vf 7 ^ j. Physically, this can be 
achieved by choosing A so that all tensor states in it can 
be sharply distinguished from each other by a set of mu¬ 
tually commuting local physical measurements. Mathe¬ 
matically these local physical measurements are Hermi- 
tian operators acting near the center of the patch where 
quasiparticles live. For instance, these measurements 
could include a measurement of the locations of chargons 
by inserting small fluxon loops. Then are chosen 

to be the eigenstates of these measurements with dis¬ 
tinct eigenvalues. Since these measurements are locally 
near the center of the patch, the boundary condition (i.e., 
the virtual boundary state) will not affect the measure¬ 
ment when the patch is large enough, and the condition 
= 0, Vi 7 ^ j is expected to hold. 

We then can construct a local operator to transform 
states in A under a symmetry R: 


Li = T.^T 






(81) 


as shown in Fig.llOr). One can easily verify, de¬ 


fined above indeed satisfies Eq.(78) for all states in A. 


Moreover, such local operators in patches satisfy 

the symmetry fractionalization condition Eq. (69) up to 


the “locally accumulated” phase outside these patches 
HsGOutside 

After the local symmetry operator is defined, we are 
able to study the symmetry fractionalization of chargons. 


Consider a relation between symmetry group elements 
R 1 R 2 ■ ■ ■ Rn = e, we can construct a local symmetry op¬ 
erators as 

^Ri...R„ — ^Ri ' " ( 82 ) 

By inserting Eq.([ 8 T]) into the above equation, we get 

Lt...R„ = ^[(f(83) 

i 

where 

fA,Ri...R„ = ^ ^ ^ Wn^Rn o 

= XR.VR o (84) 


Here, the Z 2 element 77 ^ and the phase factor xr act on 
boundary virtual legs, as shown in Fig.(10:). The second 
line of the above equation is obtained by the following 
fact: 


VR{s,i)XR{sA) = WV?i(s,i).. .fEfi„(i?„li.. .i?i ^s,!)) 

(85) 


When rjfi = I, the action of ^ on an arbitrary 
tensor T 4 G A gives the same phase. When 77 ^ is the non¬ 
trivial Z 2 element, a Z 2 odd tensor picks up an extra 
— 1 comparing to a Z 2 even tensor Tj under the action 


of . This is exactly the phenomena for symme¬ 

try fractionalization of chargons: for nontrivial 777 J, under 
symmetry Ri... Rn, a single chargon picks up an extra 
— 1 comparing to a topologically trivial excitations. 

Note that XR only serves as a global phase, thus does 
not contribute to the symmetry fractionalization of char¬ 
gons. It appears in Eq.(84) even for the ground state 
tensor patch. In fact, this result is expected and is con¬ 
sistent with the physical interpretation of x discussed in 
the previous subsection. One way to see this is to repeat 
the above analysis only for the ground states of the Id 
spin-1 AKLT model on an open chain, with the patch A 
covering one end of the chain. Here one should instead 
consider an injective matrix project state since the IGG 
here is trivial. The appearance of x in this example can 
be simply interpreted as the projective representation of 
the edge states in the AKLT model. 


V. FLUXONS AND THE DECORATED PEPS 

In this section, we will construct the decorated PEPS 
from the original symmetric PEPS with IGG = ^ 2 - The 
decorated PEPS explicitly captures all Z 2 gauge fluctu¬ 
ations, and is a good tool to study properties of fluxons. 
In particular, we can extract fractional lattice quantum 
numbers of fluxons by constructing local symmetry op¬ 
erators on a small patch of the decorated PEPS. The 
result shows that lattice quantum numbers of fluxons is 
completely determined by the chargon distribution. 
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A. The decorated PEPS 

Given a symmetric PEPS with IGG = Z 2 , we can cre¬ 
ate topological excitations such as chargons and fluxons. 
As shown in Sec |IID[ creation of chargons is by locally 
changing site tensors from Z 2 even(odd) to Z 2 odd(even) 
and the wavefunction would vanishes if we modify odd 
number of tensors on a closed manifold. Fluxons are al¬ 
ways created in pairs as ends of strings of the nontrivial 
Z 2 action on virtual legs. 

Note that there is a major difference between chargons 
and fluxons in this symmetric PEPS language: chargon 
strings are always “hidden” while fluxon strings are ex¬ 
plicitly present. Fluxon strings, which are extended ob¬ 
jects, will cause inconvenience as one tries to define local 
lattice symmetry operations on small patches with flux¬ 
ons. 

To overcome this inconvenience, we define a decorated 
PEPS, on which we can create fluxons by changing ten¬ 
sors “locally” without creating strings. The decorated 
PEPS can be viewed as a dual description of the original 
symmetric PEPS, and they represent the same physical 
state if one does not consider the boundary effects. In 
particular, we will show that fluxon strings are always 
“hidden” while the chargon strings are explicitly present 
in the decorated PEPS. For simplicity, we consider the 
decorated PEPS on the square lattice first, and develop 
the decoration method. Then we apply the method to 
the kagome PEPS with Z 2 odd site tensors. 

Now, let us discuss the method to obtain the decorated 
PEPS defined on the square lattice, as shown in Fig. 0- 
We will first discuss the case where all tensors are Z 2 
even. The procedure to construct the decorated PEPS is 
as follows: 
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FIG. 11. (a): We decorate the PEPS by changing bond ten¬ 
sors, adding plaquette tensors, while leaving site tensors in¬ 
variant. New legs connecting plaquette tensors and bond ten¬ 
sors are depicted by green dashed dotted line. These tensors 
is invariant under action of g on old virtual legs and g on new 
legs, (b): Any patch of PEPS are invariant under action of g 
and g on boundary of this patch. The whole PEPS is invariant 
under IGG x IGG. (c): By changing a plaquette tensor from 
Z 2 even to Z 2 odd (blue tensor), one creates a fluxon at the 
plaquette center. Odd number of fluxons inside a patch can 
be detected by acting p loop around the boundary new legs, 
where one gets an extra minus sign. The physical meaning of 
7j loop is the Wilson loop of chargons, and the extra minus 
sign encodes the braiding statistics of chargons and fluxons. 


states with even numbers of — I’s: 

p, = |T, T, T, T) + |T, T, -T, -T) +... (87) 

where Pc labels the tensor at a plaquette center c. 


1. One adds two new virtual legs pointing to plaque¬ 
tte centers for every bond tensor. A new leg has 
dimension D = 2. And we change the bond tensor 
Bh to Bh as 

Bb = Bb (Z \l,l) I- g ■ Bb ^ \ -T, -T), 

( 86 ) 

where | ± 1) labels the basis for the new virtual 
leg, and g is the nontrivial Z 2 element. Here in 
Eq.( ]86[ ), g is defined to act only on one end of the 
bond tensor Bb. This is shown in the second figure 
of Figpr^a). 


Let us visualize the decorated PEPS. Since only configu¬ 
rations with even number of | — 1) contribute to plaquette 
tensors and bond tensors, we conclude that | — 1) always 
form loops. Further, every loop configuration contributes 
equally to the wavefunction. To see this, consider an ar¬ 
bitrary loop configuration C. Then, we can transform 
the configuration to a PEPS wavefunction defined on the 
undecorated lattice. The PEPS wavefunction asso¬ 
ciated with this loop configuration is obtained by mod¬ 
ifying the undecorated PEPS I'll). For bond tensors in¬ 
tersecting with loops, according to Eq.(86), we have 


Bb g • Bb 


( 88 ) 


2. All site tensors are the same as those in the 
undecorated PEPS, as shown in the first figure of 
Figl^a). 

3. One adds plaquette tensors at all plaquette cen¬ 
ters as shown in the third figure of Fig{T^a). Pla¬ 
quette tensors connect the modified bond tensors 
nearby with four new virtual legs. Plaquette ten¬ 
sors are simply superpositions of all (new) virtual 


While all other tensors are unchanged. Namely, the ac¬ 
tion of g forms the same loop configuration C. In fact, C 
can be viewed as a fluxon loop in the original PEPS be¬ 
fore decoration. Since the tensor obtained by contracting 
site tensors and bond tensors inside any region is also Z 2 
even, the action of g on any loops gives the same quan¬ 
tum state. We then conclude that every loop configura¬ 
tion |1E'£) contributes the same wavefunction as I'll). The 
decorated PEPS, which is the linear superposition of all 
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|^'£), describes the same quantum state as the undeco- (a) 
rated one up to normalization (for the moment, let us 
postpone the discussion on boundary conditions). 

The decorated PEPS still has a Z 2 = {!,<?} invariance 
with g acting on the old virtual legs only. In addition, 
the single-sided g-action along a loop of old virtual l^s is 
equivalent to flipping all the new Ising variables | ± 1) —>■ 

=pl) along the corresponding loop C of the new virtual 
legs. 

In the following, we will study fluxon excitations in 
the decorated PEPS. First, let us point out that for the 
decorated PEPS, there is an additional Z 2 gauge trans¬ 
formation IGG leaving all tensors invariant. We define 
the action 'g as 

go|±T)=±|±T), (89) 




FIG. 12. (a),(b):The decorated PEPS on square lattice with 
every site tensor to be Z 2 odd. Bond tensors with blue dot 
and blue line should be modified as B'l, = g ■ Bb to ensure 
that every site is the end of some chargon string. One can 
verify (a) and (b) are gauge equivalent, thus represent the 
same state. 


while g acts trivially on the original virtual legs. The 
nontrivial element in IGG = Z 2 is the action of g on all 
new virtual legs only. Then, if the action of g on all new 
virtual legs of a tensor leaves the tensor invariant, we call 
this tensor a Z 2 even tensor. Similarly, we can define the 
Z 2 odd tensor. To create a fluxon living on plaquette 
center c, one can simply change the plaquette tensor 
from Z 2 even to Z 2 odd. For instance, a plaquette ten¬ 
sor PJ" which supports one fluxon projects out quantum 
states with even numbers of | — 1): 

p™ = I - T,T,T,T) + I - T, -T,-T,T) +... (90) 

In order to see that P™ indeed supports a fluxon, we 
translate back to the undecorated PEPS. Pc only has con¬ 
figurations with odd numbers of | — 1), Thus, in the un¬ 
decorated lattice, there are always odd numbers of mod¬ 
ified bond tensors g • Bb around the plaquette c. In other 
words, a single fluxon lives at the plaquette c. 

Fluxons inside a small patch can be detected by acting 
g on boundary of that patch, see Fig. ©)■ An odd num¬ 
ber of fluxons contributes an additional (—1) due to the 
Z 2 oddness of the patch tensor. It is natural to interpret 
the loop of g (acting only on one end of an involved bond) 
as the chargon loop which detects fluxons. One can eas¬ 
ily show that this is exactly true and the end points of 
an open g-string are chargons. In the decorated PEPS, 
the chargon strings are explicit while the fluxon strings 
are “hidden”. 

Before we discuss the symmetry fractionalization of 
fluxons, let us study the case with site tensors being Z 2 
odd. In this case, if one exactly follows the above decora¬ 
tion procedure, one finds that the similarly constructed 

wavefunction satisfies: 

I'kc) = (-1)”1^) (91) 

where Us is the number of sites enclosed by the loop £. 
To construct a decorated PEPS with in-phase contribu¬ 
tions from all Ifflc), we require an additional —1 for loops 
enclosing odd numbers of sites, and need to modify the 


decoration procedure. As shown in Fig.(12), we simply 


modify some bond tensors by the action of g only on one 
of the two new virtual legs. More precisely, we can choose 
bond tensors in x direction to be modified for every other 
column in the following way: 

Bb = g-Bb 

= |T,T)-g-Pb(g)|-T,-T) (92) 


Note that the g action here is defined to only act on one 
of the two new virtual legs, and picks up the (—1) for 
the second term in the second line. One can easily verify 
that for the modified decorated PEPS, loop configura¬ 
tions enclosing odd numbers of sites contribute the same 
wavefunction as the undecorated PEPS. 

As we pointed out before, g strings can be interpreted 
as chargon strings. Further, the Z 2 oddness of site ten¬ 
sors can be interpreted as one chargon per site. The 
action of g on modified bond tensors in fact creates short 
“chargon strings” connecting background chargons. A 
loop of I — 1) enclosing an odd number of sites intersects 
with chargon strings for an odd number of times, and 
contributes an extra —1. Keeping this in mind, one can 
easily construct other possible decorated PEPS, once en¬ 
suring that every background chargon (Z 2 odd tensor) 
is an end of a g string, as shown in Fig. We point 
out that all different decorated PEPS obtained for the 
same quantum state are gauge equivalent, in the sense 
that they can be transformed to each other by acting g 
on both ends of a collection of virtual legs. 

The above point of view is very useful for us to con¬ 
struct the decorated PEPS in more complicated lattices 
such as the kagome lattice. In the kagome lattice, pla¬ 
quette centers form a dice lattice, as shown in Fig.(131. 


Following the similar procedure as in the square lattice 
case, we decorate the PEPS by adding plaquette tensors 
and changing bond tensors. The new virtual legs are 2- 
dimensional Hilbert space, with basis | ± 1) . Note that 
unlike the previous case, there are three kinds of plaque¬ 
tte tensors. Two lie at centers of triangles while one lies 
at the honeycomb center. These tensors project out con- 
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FIG. 13. The decorated PEPS for the kagome lattice with 
Z 2 odd site tensors. Fluxons live on the dual lattice, which 
is the dice lattice in this case. There are three plaquette 
centers in one unit cell, labeled as p, q and r. Green/blue 
legs are two dimensional Hilbert space with basis as | ± 1). 
Plaqurtte tensors project out configurations with odd number 
of I — 1). Bond tensors need to be modified in two different 
ways according to their color. 


be identity. While for the case with site tensors being 
odd, Wxi is nontrivial. If we decorate the PEPS as in 
Fig.(lT2}t), we get 

WTi(a:, j/,i) = P, {x, y, *) = I (95) 

where i labels the four new virtual legs coming out of 
the plaquette tensor at {x,y), and J = (o-?i) is the 
representation of g on new virtual legs. If the decoration 
has the form as in Fig. ([II 3 ), we get 

Wti (x, y, i) = I, Wt 2 [x, y, i) = (96) 

As we can see, values of Wt^ depend on the way to deco¬ 
rate the PEPS. However, similar to the undecorated case, 
we have 


figurations with odd numbers of | — 1). Further, due to 
the Z 2 oddness of site tensors, we should ensure that ev¬ 
ery site tensor in the decorated PEPS is connected with 
a chargon string. Thus, we can modify bond tensors as 

Bb = Bb® |T,T) + g ■ Bb(^ {\-T, -T)) 

B'b = BbZ) |T,T) + g-Bb®g-{\ -T,-T)) 

= H6®|T,T)-5-56 0(1-1,-1)), (93) 

where the pattern of bond tensors is shown in Fig. Q. 


B. Symmetry fractionalization of fluxons 

In this part, we will develop a general method to ex¬ 
tract (fractional) quantum numbers carried by fluxons. 

Comparing with the undecorated PEPS, there are 
more gauge freedoms associated with new virtual legs for 
the decorated PEPS. We call these new gauge freedoms 
as V. For the purpose of the discussion in this section, 
we only need to consider V gauge transformations such 
as acting g on both ends of a number of new virtual legs, 
which leave the physical wavefunction intact. If the quan¬ 
tum state is invariant under a symmetry i?, tensors of the 
decorated PEPS should satisfy the following conditions: 

T® = QrWrRoT^ 

Bb = WrWrR o Bb 

Pc = WrRoP^ (94) 


W^^\T 2 {x, y, i))W^^\T^T 2 {x, y, i))Wt, [T^T^i^x, y, i)) 
WTAT2^T^T2{x,y,i)) = m2 (97) 


as a gauge invariant quantity. Inserting Wp. into the 
above equation, we conclude that rji 2 = I for the case 
with site tensors of undecorated PEPS being Z 2 even, 
while rji 2 = J for the case with site tensors being Z 2 odd. 

Next, we will show that rji 2 is directly related to the 
translation symmetry fractionalization of fluxons. To see 
this, let us consider the decorated PEPS with only fluxon 
excitations inside some small patches. (Note that the 
decorated PEPS is inconvenient to study chargon exci¬ 
tations since chargon-strings are explicit.) Following the 
similar procedure in Sec.(IVB), one can construct local 
translation operators for a single patch A of the deco¬ 
rated PEPS, labeled as , * = 1,2. Labeling the state 
associated with A as T^, we get 


M WTJ^of^ (98) 

By series connecting local symmetry operators, we can 
define 

^ • (M)-^ • Li ■ Li (99) 

as the local operator associated with T 2 ^T^'^T 2 Ti. Act¬ 
ing this operator on T^, we have 


fA 

T-^T-^T^Ti 


■ = X 12 V 12 V 12 o fj 


( 100 ) 


where Wr labels the gauge transformation associated 
with symmetry R on old virtual legs while Wr labels 
that on new legs. Here Wr takes the same value as in 
the undecorated PEPS. And we will solve Wr in the fol¬ 
lowing. 

Let us firstly consider a simple example: the transla¬ 
tion symmetry group generated by Ti , T 2 for the deco¬ 
rated PEPS defined on the square lattice. For the case 
where site tensors are Z 2 even, we have Wt^ and TErb to 


For the case with nontrivial yi 2 , if there are odd num¬ 
bers of fluxons inside the patch A, ^ will pick 

J -2 -^1 L2L1 

up an extra —1. This indicates the nontrivial transla¬ 
tional symmetry fractionalization of fluxons. Since rji 2 
only depends on the Z 2 parity of site tensors, we con¬ 
clude that the translation symmetry fractionalization of 
fluxons is fully determined by the background chargon 
distribution. The above argument can be easily general¬ 
ized to arbitrary lattice symmetry operations. 
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In the following, we will figure out Wr of the decorated 
PEPS for the kagome lattice example. As we discussed 
above, this directly implies the symmetry fractionaliza- 
tion pattern for fluxons. One can easily work out sym¬ 
metry transformation rules directly from the decorated 
PEPS shown in Pi^ The result is listed as follows. 

Wri(a;,?/,s,i) = J^, 

WT2{x,y,s,i) = I, 

WcAx,y:P/q,i) = 

Wc,{x,y,r,i) = ry+"^^^^+^'>+^+y, 

wAx,y,p,i) = ry+\ 

WAx,y,q/r,i) = ry, (101) 


where s' can be any one of the plaquette sublattices la¬ 
beled by pjqjr as shown in Eigp^ and i labels the new 
virtual legs coming out of the plaquette tensor. These 
symmetry transformation rules are not gauge invariant. 
The symmetry fractionalization of fluxons is determined 
by gauge invar ia nt qu a ntitie s p’s. By re placi ng Wr in 
Eq.pil, IbTII , IbItI , liji] , ( |b 35| and ( |b4i| with Wr 
obtained above, we find 


Vl2 = VTicr = r]T2a- = J, 

VTiCe = lT 2 Ce = VCe = Vcr = PaCe = I- (102) 

For on-site symmetries such as the spin rotation and 
the time reversal symmetry, symmetry transformation 
rules on new virtual legs are trivial. In other words, we 
expect fluxons constructed here to be spin 0 as well as 
Kramer singlet. Our result for symmetry fractionaliza¬ 
tion of fluxons is consistent with eariler results in ReflTZl 
andiTbl 


VI. SYMMETRIC PEPS ON TORUS AND 
LONG-RANGE ORDER 


In the above discussion, we mainly focus on symmetric 
PEPS on infinite lattices. In this section, we will consider 
the symmetric PEPS on a finite torus. Namely, site ten¬ 
sors on the left (up) boundary are connected to sites on 
the right (down) boundary by bond tensors. Further, we 
will provide some of our partial understandings on how 
long-range ordered phases fit into the current symmetric 
PEPS formulation. 

To make the discussion concrete, we will focus on spin- 
1 systems to demonstrate the principle. We firstly con¬ 
sider the PEPS description of a Z 2 spin liquid phase on 
finite tori in subsection [VIA and construct the topo¬ 
logical degenerate ground state sector. Then we study 
the finite size effects due to the spinon (vison) condensa¬ 
tion in subsection VIB which gives rise to the MO (VBS 
order) in the long range. One result here is that the un¬ 
decorated PEPS provides the natural basis to represent 
a MO state, while the decorated PEPS is the natural 
language to represent a VBS state. 


A. Topological degeneracy in the PEPS 
formulation 


It is well known that for the toric code topological or¬ 
dered phase, the ground state degeneracy (GSD) equals 
four on torus in the thermodynamic limit. As mentioned 
in subsection IID1 in the Z 2 invariant PEPS, one can 
construct these four states on a finite torus by acting the 
nontrivial Z 2 element q on the non-contractible loops of 
torus, as shown in Fig.((T^).l^ We label these four states 
as I'ko.o)) 1^0 ,tt) and |'k,r, 7 i-)- Recall that g strings 

can be interpreted as flux loop. So these four ground 
state basis can be visualized as different ways of insert¬ 
ing non-contractible flux loops. We call this set of basis 
for the ground state manifold as the m-basis. Note that 
in general, these four states have different energies on 
any finite torus. However, if the system is in the de- 
confined phase, the energy difference between these four 
states goes to zero in the thermodynamic limit. 


Now, let us consider decorated PEPS defined in Sec|V| 
The decorated PEPS describes the same wavefunction as 
the undecorated one on the infinite plane. However, this 
is not true if we consider finite samples on tori. Using the 
method developed in the last section, we can easily con¬ 
struct the decorated PEPS on a torus based on an undec¬ 
orated PEPS state |'ko,o)- Similar to the infinite PEPS 
case, any configurations with all | — 1) forming loops will 
contribute to the wavefunction of the decorated PEPS. 
For a torus sample, we should consider both contractible 
loops and non-contractible loops of | — 1). First, any con¬ 
figurations with only contractible loops contribute l^'o.o) 
to the decorated wavefunction. Note that for configura¬ 
tions with an even number of non-contractible loops, one 
can always decompose them to only contractible loops. 
For configurations containing an odd number of non- 
contractible loops in the x/y/hoih direction, one gets 
±|^ 7 r,o)/±|^o, 7 i-)/±|^ 7 r, 7 r)- Here, the ± signs depend on 
the way of decoration: as pointed out in Sec. to cap¬ 
ture the background chargons of original PEPS, some 
bonds are modified to B'f^ = g ■ Bh. When loops of | — 1) 
intersect with bond we will get an extra —1. Thus, 
it is always possible to choose the distribution of such 
that the ± signs are -|-1. The obtained decorated wave- 
function |'I'o,o) is then 


1 ^ 0 . 0 ) = I^'O.O) + |^,r,o) + (103) 


up to a normalization factor. 

Other three states in the decorated language can be 
generated from |'I'o,o) by threading chargon strings on 
non-contractible loops of torus, or in other words, by act¬ 
ing nontrivial Z 2 element g on new virtual legs along non- 
contractible loops. We label them as and 

respectively, as shown in Fig.(14D). It is straight- 
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FIG. 14. Topological degenerate ground states on finite torus. 
Site tensors at left (up) boundary are connected with those at 
right (down) boundary . (a) labels m-basis while (b) labels 
e-basis. 


forward to see that 

I^TT.o) = l^'o.o) + 

1^0,^) = I^'o.o) - 

(104) 

Here, the basis formed by four jd!) states is named as the 
e-basis. 

As concrete examples, in Appendix we explicitly 
constructed the four-fold ground state sectors of the 
Qi = Q2 QSL and Qi = —Q2 QSL on an even by odd 
(4n -I- 2)-unit-cell torus sample. Interestingly, on such 
samples, we find that the ground state sectors of the two 
QSL share no identical lattice symmetry irreducible rep- 
resentations(irreps); namely, any irrep in the Qi — Q 2 
QSL ground state sector is different from any irrep in the 
Qi — —Q2 QSL ground state sector. This result is consis¬ 
tent with a recent study based on parton construction^^. 


B. Long-range ordered phases represented by 
symmetric PEPS 

We claimed that each crude class contains many pos¬ 
sible member phases, and these phases are distinguished 
by long-range physics. Here we mention some concrete 
examples. In the 32 crude classes on the kagome lattice, 
we know that each crude class contains one Z 2 QSL mem¬ 
ber phase which has no spontaneous symmetry breaking. 
Let’s consider the crude class which contains the Qi = Q 2 
Z 2 QSL, firstly constructed using the Schwinger-boson 
approactPSl. There are a few known neighboring phases 
of this QSL in which symmetry is spontaneously broken 
in different fashions. These symmetry breaking phases 
and the Qi = Q 2 Z 2 QSL are in the same crude class 
proposed in this work. 

For example, the 120° Q = 0 long-range magnetic or¬ 
dered (MO) state can be obtained if the sp inons in the 
Qi = Q2 QSL condense in the long rang^^MZI And 
a valence bond solid(VBS) state with 12-site per unit 
cell, which breaks the translational symmetry, can be 


obtained if the visons in the Q\ = Q 2 QSL condense in 
the long rang^^. To capture the long-range physics of 
these symmetry breaking phases, we expect that scaling 
of bond dimension D and scaling with system sizes need 
to be performed in our symmetric PEPS methods. For 
instance, it was shown that generically the entanglement 
entropy in symmetric states which breaks a continuous 
symmetry in the long-range contains additive logarith¬ 
mic corrections^. To our knowledge, there is no known 
PEPS construction with a finite D that can be proven to 
host this behavior. 

Next, let us try to understand the effects of the 
spinon/boson condensations in the PEPS formulation on 
finite tori. Note that practical PEPS simulations on a 
torus sample can be computationally very expensive, so 
the discussion here is mainly for conceptual purposes. 
Before the condensation, the system is in a Z 2 spin liq¬ 
uid phase with a four-fold ground state sector "Hgs- Af¬ 
ter the vison-m condensation, the system is expected to 
have a multi-fold ground state sector due to the lattice 
symmetry breaking, and the number of ground states de¬ 
pend on the geometry of the sample and the VBS order 
pattern. On the other hand, after the spinon-e condensa¬ 
tion, there will be gapless Goldstone excitations assuming 
spin-rotational symmetry. Exactly how the ground state 
sector on a torus evolves as the spinon/vison condenses? 

Below we provide our partial answer to this interest¬ 
ing question. Note that, in the ordered phases, for the 
system to “know” the spinon/vison condensation, we al¬ 
ways consider torus samples whose sizes are larger than 
the vison/spinon confinement length scale. Eor simplic¬ 
ity, we will focus on samples whose geometries are com¬ 
mensurate with the spatial patterns of the long-range 
orders. We can consider a Ginzburg-Landau theory de¬ 
scribing the boson condensation phase transition. In 
the Ginzburg-Landau language, this commensuration is 
achieved by choosing the finite lattices together with 
a proper boundary condition such that the momentum 
space minima of the condensed bosons are available. 

Here one immediately sees that the undeeo- 
rated(decorated) m-basis (e-basis) is exactly suitable 
to describe the spinon(vison) condensation. Eor in¬ 
stance, the existence of m-loops trapped in the torus 
holes is explicit in the undecorated m-basis. Con¬ 
sequently the boundary conditions for the spinon 
condensation Ginzburg-Landau theory are sharply 
defined. We conclude that only one of the four states 
in the m-basis is representing the true ground state in 
the MO phase. The spinon condensation minima in the 
momentum space are clearly not available for the other 
three states due to the J momentum shift. Physically, 
we know that the other three states trap MO vortices 
in the torus holes. Therefore they have finite excitation 
energies (but zero excitation energy density) in the MO 
phase based on simple nonlinear sigma model analysis. 

Similarly, the decorated e-basis is suitable to describe 
the VBS order, since the boundary condition for the con¬ 
densed visons is explicit. Only one of the four states 
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in the e-basis corresponds to a true ground state in the 
VBS phase, while the other three states host VBS do¬ 
main walls wrapping around the torus holes. Note that 
the confined spinon-e corresponds to a VBS vortejiP^ and 
the e-string corresponds to a VBS domain wall in the 
VBS phase. These three states are separated from the 
true ground state by excitation energies proportional to 
the linear sample size due to the energy cost of the do¬ 
main wall. 

Note that the decorated PEPS only captures one state 
in the true ground state sector in the VBS phase. The 
other “cat states” in this symmetry breaking phase are 
still missing in the current construction. In fact on a 
torus sample with even by even unit cells, our construc¬ 
tion leads to physical states with the center of mass mo¬ 
mentum at P = (0,0) only. We expect that, in order to 
capture the other cat states at different center of mass 
momentum, one needs to perform the symmetric PEPS 
classification on finite tori instead of on the infinite plane. 

Next, we discuss some general guiding principles about 
these symmetry breaking member phases. It is crucial 
to note that these symmetry breaking patterns in the 
vicinity of the Q\ = Q 2 QSL are not arbitrary. In fact, 
the sharp way to understand how the symmetry breaking 
patterns arise is exactly to use the Ginzburg-Landau the¬ 
ory — the golden tool to investigate symmetry breaking. 
In order to write down a Ginzburg-Landau theory for the 
spinon-e (vison-m) condensation, the only information 
that we need to know is how these particles transform 
under the global symmetry group, which is nothing hut 
the projective symmetry group transformation rules for 
spinons and visons investigated in Sec \IV\ and S'ec]^ In 
the past, this is exactly how the 120° Q = 0 MO state 
and the VBS state with a 12-site unit cell were identified 
to be the neighboring phases of the Qi = Q 2 ^2 QSL. 

Therefore, even though we emphasize that our algo¬ 
rithm can be used to efficiently determine the crude 
class of the quantum ground states based on short range 
physics, we still learn sharp constraints on candidate 
symmetry breaking patterns in the long range physics: 
these orders must be consistent with the Ginzburg- 
Landau theory of the given crude class. With a careful 
scaling with larger system sizes and bond dimension D, 
our algorithm could be practical useful to pin down the 
possible long range orders in quantum phase diagrams. 

Note that the spinon condensation and vison conden¬ 
sation in the 32 Z 2 QSL studied here are quite differ¬ 
ent: The rj indices are really the symmetry transforma¬ 
tion rules for spinons, while the symmetry transformation 
rules for visons are completely fixed. Gonsequently, two 
QSL with different p indices are expected to connect to 
different MO orders after spinon condensations, since the 
Ginzburg-Landau theories are different. For instance, the 
Qi = —Q 2 Z 2 QSL is connected to the -s/S x -s/S MCP^, 
fundamentally different from the 120° Q = 0 order in the 
vicinity of the Qi = Q 2 QSL. 

However, the same Ginzburg-Landau analysis indi¬ 
cates that all the 32 Z 2 QSL could give rise to the same 


long-range VBS order pattern after vison condensation. 
There are two possible explanations for this phenomenon: 
(I) it is possible that the long range VBS orders emerg¬ 
ing from different QSL, although sharing the same real 
space pattern, are still in different quantum phases. (2) 
it is possible that certain VBS phase could appear in dis¬ 
tinct crude classes .1^ 

We tend to believe that either scenario could be correct 
under certain conditions, although we do not have rigor¬ 
ous understandings. It would be very helpful to perform 
numerical simulations based on the algorithms proposed 
here in models supporting relevant symmetry breaking 
phases and see exactly what happens. But we have to 
leave this as a subject of future investigations. At this 
moment we only can provide some physical speculations. 

First, the Ginzburg-Landau theory for vison or spinon 
condensations completely misses x indices in the classifi¬ 
cation, and the long-range real space VBS pattern does 
not capture the physics described by these indices. Gon¬ 
sequently we expect that scenario-(I) is correct if differ¬ 
ent classes have different x indices. In fact, we could 
compare two classes with XT = +I and XT = respec¬ 
tively. The Z 2 QSL member phase in the second class is 
expected to host nontrivial projective representations of 
SU{2) X T on the physical boundary. These boundary 
degrees of freedom are expected to lead to measurable 
effects even if the bulk VBS order is established, which 
is absent in the first class. 

Second, we speculate that scenario-(2) could be correct 
when different classes share the same x indices but have 
different p indices, e.g., the Qi = Q 2 QSL and the Qi = 
—Q 2 QSL. Note that this speculation is not as naive as 
it appears. In particular, in Appendix[D| we demonstrate 
that the 4-fold ground state sectors of these two QSL have 
completely different lattice symmetry irreps on (4n -f 2)- 
unit-cell samples. Therefore the symmetric PEPS of the 
two classes in the VBS ordered phase must be describing 
distinct ground state wavefunctions on such samples. In 
Appendix we trace the origin of the lattice quantum 
number discrepancy of the two PEPS classes in the 12- 
site VBS phase, and show that it is related to the phase 
factor due to the quantum fluctuation of valence bonds 
along a VBS domain wall. However, we expect that this 
particular phase factor is not a universal feature, and the 
VBS orders in the two classes correspond to the same 
quantum phase. 


VII. DISCUSSION AND CONCLUSIONS 

In this paper we attempt to construct generic sym¬ 
metric ground state wavefunctions for integer or fraction¬ 
ally filled correlated systems using PEPS, under certain 
assumptions. Here we review the assumptions that we 
made and discuss the limitations and generalizations of 
our results. 

We firstly highlight an assumption that is, to our 
knowledge, due to a more fundamental difficulty. And 
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we currently do not know how to solve this difficulty gen¬ 
erally. This assumption is that the on-site symmetry is 
implemented as the simple tensor product of local rep¬ 
resentations or projective representations on the virtual 
legs in PEPS. For instance, this is the origin of the mini¬ 
mal required IGG in the half-integer spin systems on 
the kagome lattice. 

This assumption is known to have problems at least in 
the long range physics when attempting to describe SPT 
phases. For instance, let’s attempt to construct a U(l) 
charge-conserving Chern insu lator using the fermionic 
version of PEPS (fPEPSjP^^^^^^. Here the exact con¬ 
structions of free fermion sta tes w ith a nonzero Chern 
number using Gaussian fPEPS^i^, in which the virtual 
legs transform as U(l) representations, are shown to host 
power law correlation functions in the real space. It has 
been pointed by Hasting^ that for a general U(l) sym¬ 
metric PEPS with a bounded bond dimension D which 
is a fully gapped ground state of a local Hamiltonian, the 
assumption that the virtual legs transform as U(l) rep¬ 
resentations and the assumption that the PEPS carries 
nonzero Chern number generically lead to contradictions. 

Another important piece of information can be ob¬ 
tained by understanding the exact PEPS constructions of 
available short-range correlated SPT ground state wave- 
functions of exact solvable model^^. In particular, for a 
finite on-site symmetry group, it is shown that the virtual 
legs do not form representation (or projective represen¬ 
tations) of the symmetry. Instead the virtual degrees of 
freedom transform in a “non-on-site” fas hion, which can 
be described by matrix product operator^SlM!, 

We currently do not know how to generically repre¬ 
sent an SPT state in two and higher dimensions with 
correct long-range physics using PEPS. However, it is 
possible that our assumption about symmetry represen¬ 
tations on virtual legs does not cause problems in captur¬ 
ing the short-range physics of SPT phases under certain 
conditions. For instance, given an finite size sample, it is 
possible that the SPT ground state can be accurately ap¬ 
proximated with a PEPS after a scaling with respect to 
bond dimension D is performed. As demonstrated in an 
example using Gaussian fPEPS in ReflS^ the required 
bond dimension D in practical simulations on intermedi¬ 
ate sized sample may not be very large. 

Moreover, we speculate that even the short-range 
physics of a SPT phase may not be captured using the 
current symmetric PEPS construction. For example, 
it is known that for inversion symmetric system, the 
ground state wavefunction of a Chern insulator with an 
odd Chern number is inversion odd.^ And the inver¬ 
sion quantum number should be completely short-range 
physics. 

We made a second assumption: we study only those 
symmetric quantum ground states that can be repre¬ 
sented by a single tensor network on the infinite lat¬ 
tice. This assumption is made here mainly for technical 
simplicity rather than fundamental difficulty. Note that 
this assumption is weaker than the assumption that the 


ground state sector is composed of one-dimensional rep¬ 
resentations of the symmetry group on any finite size 
samples. For instance consider a Z 2 QSL studied in 
this paper with a four-fold ground state sector on tori. 
When considering a finite size torus, some of them could 
form multi-dimensional irreducible representations of the 
space group. 

This assumption could be violated in general model 
simulations. As a trivial example we could consider a 
ferromagnetic state in an SU{2) symmetric model. In 
this case the number of degenerate ground states scale 
linearly as the number of sites, which certainly cannot 
be represented by one or few PEPS. 

As a slightly nontrivial example, we refer to the chiral- 
spin-charge-Chern liquid (SCCL) in Refl^ The spin 
dynamics in SCCL is described by a chiral Z 2 QSL, 
which is a Z 2 QSL breaking the time reversal symme¬ 
try and has nonzero spin-chirality order parameter (e.g., 
< Si ■ Sj X Sk > 7 ^ 0 for three nearby spins i,j,k.). 
This state breaks both time-reversal and mirror reflec¬ 
tion symmetries, but leaves the combination of the two 
respected. In this situation, we found 8 = 4x2 ground 
states on symmetric torus samples (compatible with the 
PSG transformations). The factor of 4 is related to the 
topological degeneracy of Z 2 gauge theory. And the extra 
factor of 2 is due to the fact that the time reversal, the 
mirror reflection and the lattice rotation form nontrivial 
2-dimensional irreducible representations. The latter fact 
dictates that it is impossible to represent such chiral liq¬ 
uids by a single symmetric PEPS, in which case the extra 
factor of 2 degeneracy cannot be captured. The simple 
way to proceed is to instead only consider the combina¬ 
tion of the time reversal and the mirror reflection as a 
symmetry, which allows a description of one of the two 
time-reversal images using PEPS. The PEPS description 
of the other state can be obtained by the time-reversal 
transformation. 

We now comment on another fact in our construction. 
In the half-integer spin systems on the kagome lattice, we 
show that a spin-singlet symmetric PEPS has an IGG 
that at least contains a Z 2 subgroup. If IGG = Z 2 for 
a PEPS, and if the PEPS is describing a fully gapped 
QSL, we showed that the topological order is toric-code- 
like in Sec IHD 11 This remains to be true if we con¬ 
struct some Z 2 QSL in the absence of the time-reversal 
symmetry, using ou r formu lation. However, there are 
known construction^i^OHioaj gapped Z 2 QSL on the 
kagome lattice in the absence of the time-reversal sym¬ 
metry whose topological order is the same as the one in 
the double-semion model, fundamentally different from 
toric-code. 

Interestin gly, in a PEPS construction of the double 
semion QSLP^, in which spin rotation is still imple¬ 
mented as representations on the virtual legs, the con¬ 
structed tensors are actually Z 4 invariant. Naively, such 
a state should have a 16-fold degenerate ground state sec¬ 
tor on torus, but it was shown that only 4 of them are 
linearly independent. 
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Next we comment on the connection between our work 
with previous works. For readers that are familiar with 
the parton constructions and pro jective s ymmetry group 
analysis of parton wavefunctions^^Mlll^ clearly part of 
our results can be viewed as generalizations of these 
analyses into PEPS wavefunctions. In particular, in the 
kagome half-integer spin S example presented here, ev¬ 
ery crude class contains a distinct Z 2 QSL as a member 
phase. Part of our results can be viewed as a classification 
of Z 2 QSL on the kagome lattice. Comparing with previ¬ 
ous investigations on this top ic spe cifically for S = 1/2, 
based on parton construction^^ZEH^ -vve find that our re¬ 
sult captures every phase present in the Schwinger-boson 
const ructiorP^, and finer than that. Basically the previ¬ 
ous PSG analysis of the Schwinger boson construction is 
related to the ry-indices and ©-indices in our formulation, 
while in this work y-indices are revealed. 

However, comparing with the classihcation based on 
the Abrikosov-fermion construction of Z 2 QSL on the 
kagome lattic^^, we find that some of them cannot be 
described in our result. Similar observation was made 
by Hef l75l when directly comparing Schwinger-boson and 
Abrikosov-fermion constructions. We currently do not 
have a full understanding of the physics behind this phe¬ 
nomenon. But it is worth pointing out that the missing 
Abrikosov-fermion Z 2 QSL are all found to be g apless (at 
least perturbatively) on the mean-field levepS. 

Finally we comment on the hierarchical structure of 
the crude classes. Sometimes there are physical reasons 
to believe that the IGG needs to be larger than the min¬ 
imal required one in order to correctly capture certain 
quantum phases. The double semion PEPS mentioned 
above may be viewed as such an example. 

A more conventional example in which this is ex¬ 
pected to happen is the collinear Neel ordered phase on 
the square lattice, for which we expect IGG = 17(1) 
in our PEPS constructiom In fact, the non-compact 
CP^ (NCCP^) descriptioip22l for the Neel state signals 
that the natural gauge dynamics in this state, although 
Higgs’ed out in the long-range, should be 17(1). 

One can ask the following question: for instance, sup¬ 
pose we have one parent crude class of symmetric PEPS 
with a large IGGi, what are the possible descendant 
crude classes with a smaller IGG2 C IGGG Simi¬ 
lar questions we re inv estigated in the context of par- 
ton construction^^HUnil^ Generally one expects that there 
could exist multiple descendant crude classes with IGG 2 , 
which eventually gives a hierarchical structure of the 
crude classes with different IGG’s. This hierarchical 
structure may be useful to understand certain exotic 
quantum criticalities. For example, two member phases 
belong to distinct crude classes could be smoothly con- 
nect ed v ia a critical point described by their parent crude 
clas^^^. 

As one can see from the above discussions, the cur¬ 
rent work, which is based on the point of view of di¬ 
agnosing ground state wavefunctions using symmetric 
PEPS, brings up many open questions and needs fu¬ 


ture investigations to clarify. In addition, the algorithms 
proposed here for simulating strongly interacting mod¬ 
els need benchmark tests to have a understanding of its 
practical performance. Nevertheless we believe that sep¬ 
arating the short-range part of the physics from the long- 
range part is a useful idea in investigating quantum phase 
diagrams of strongly correlated systems. While generally 
the long-range part is still a difficult task, we expect that 
the method introduced here can be used to provide sharp 
information for the short-range physics efficiently. 

We thank for help discussions with Fa Wang, Jung- 
Hoon Han, Hyunyong Lee, Panjin Kim, Xie Ghen and 
Michael Hermele. This work is supported by the Alfred P. 
Sloan fellowship and National Science Foundation under 
Grant No. DMR-1I51440. 


Appendix A: Symmetry group of the kagome lattice 


As shown in Fig.Q, we label the three lattice sites in 
each unit cell with sublattice index {s = u,v,w}. Fur¬ 
ther, we specify the virtual index {i = a, &, c, d} of a 
given site. We choose Bravais unit vector as di = i and 
«2 = ^{x + VSy). Thus, we are able to specify the virtual 
degrees of freedom of site tensors as (x, y, s, i). The sym¬ 
metry group of such a two-dimensional kagome lattice is 
generated by the following operations 


Ti : {x,y,s,i) ^ {x + l,y,s,i), 

T2 : ix,y,s,i) {x,y+ 

(7 : {x,y,u,i) -t {y,x,u,iai), 

{x,y,v,i) -t {y,x,w,icr2), 

{x,y,w,i) {y,x,v,ia2), 

Ce : {x, y, u, i) {-y + l,x + y - 
{x,y,v,i) {-y,x + y,w,i)- 

(x, y, w, i) [-y + l,x + y, u, icj- 

together with time reversal T. Here, 

{acri,6cri,Co.i,do-i} = {d, c,b,a} 
{adr2, 6 <t2, Ca2,dcr2} = {c, d, O, b} 
{aCe > bce : Cce = {b, a, d, c} 


(Al) 


(A2) 


The symmetry group of a kagome lattice is defined by 
the following algebraic relations between its generators: 


T^^T^^T 2 Ti = e, 
tT“^Tf^CTT 2 = e, 
a~^T 2 ^crTi = e, 
G^^T^^CeTi = e, 

— Irr-i—1 




TiCeTa = e, 
(7 ^CqctCq = e, 
Cq — a — / — e, 


(A3) 


g T 5 T= e, V 5 = Ti, 2 ,cr,C '6 
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where e stands for the identity element in the symmetry 
group. 

Further, consider system with spin rotation symmetry 
operator Ro-n, which means spin rotation about axis ft 
through angle 9. We mainly consider half-integer spins 
{SU{2) symmetry) in this paper. The spin rotation sym¬ 
metry commutes with all lattice symmetries as well as 
time reversal symmetry: 

9 ~^R9n9Ren = e, Vg = Ti^ 2 ,cr,Ce,T 

(A4) 


Appendix B: Classification of PEPS wavefunction 
with Z 2 IGG on kagome lattice 

In this Appendix, we will classify symmetric PEPS 
wavefunctions defined on kagome lattice with a half¬ 
integer spin per site. We first obtain the symmetry trans¬ 
formation rules for all different classes. Then using the 
symmetry transformation rules, we can solve the con¬ 
straint Hilbert space for all classes. 


as well as 


0 T 2 y, s))0Ti {TiT 2 {x, y, s)) 0 t 2 {TiT 2 {x, y, s)) 
0Ti(F2"^TiT2(a;,y,s)) = yi2W_Xi2{x,y,s,i) (B3) 


where 7712 S {bJ}) and Xi 2 {x,y,s,i) is a bond de¬ 
pendent phase. Under ambiguity Wt^ —ETiHri, 
0Ti —et; 0Ti, we get 

X 12 2 / + 1, s, i)eT^ {x + l,y + l, s, i)- 

ETs (a; -b 1, 2 / -b 1, s, i)eTi {x + l,y,s, i)xi 2 ix, y, s, i) (B4) 


Thus, we are able to set Xi 2 (a:, y, s, i) = 1, V(x, y, s, i). 

According to Eq.(32) and Eq.(34), using gauge 
transformation V(x, y, s, i) and phase factor ambiguity 
$(x, 2 /,s), we get 


WT^{x,y,s,i) V{x,y,s,i)WT 2 {x,y,s,i)V ^{x,y- l,s,i) 

'dT 2 {x,y,s) eT 2 {x,y,s)<^{x,y,s)<^*{x,y - l,s) 

(B5) 


1. Solving symmetry operation on PEPS 


In this part, we will work out symmetry transformation 
rules for PEPS defined on a infinite kagome lattice with 
a half-integer spin per site. We will focus on the case 
with minimal required IGG, which equals Z 2 , as shown 
in the main text. Eurther, every site tensor is a Z 2 odd 
tensor as we will see later. 


As shown in Sec.(BIB), the representation of Z 2 IGG 


on virtual legs can be set as the same form {I, J}, where 
J = Idj 0 (—Idj) with a ±1 ambiguity. The remain¬ 
ing U-ambiguity V{x,y,s,i) commute with J. Namely, 
V(x, y, s, i) is block diagonal with two blocks with blocks’ 
size to be Di and D 2 respectively. Further, for any sym¬ 
metry R, we have proved the associated Wr transforma¬ 
tion should also commute with J. 


a. Implementation of lattice symmetry on PEPS 


For completeness, we copy the calculation for transla¬ 


tion symmetry transformation rules done in Sec.(BIB|. 


According to the definition of symmetric PEPS, for R 
{i = 1,2) transformation, site tensors and bond tensors 
satisfy the following condition: 

rpG,y,s) ^ o y{x,y,s) 


R{xysi\x'y's'i') ^ R(xysi\x'y's'i') 


(Bl) 


Prom group relation T 2 ^T 2 ri = e, we have 
W:f:^{T 2 {x, y, s, i))W:f:^{TiT 2 {x, y, s, i))WT 2 (7’iT2(x, y, s, i)) 
WT^{T2^TiT2{x,y,s,i)) = 7712X12(2:, 2/, s, z) 

(B2) 


Then, we are able to set Wra ( 2 ;, 22 , s, i) = I as well as 
0T2 ( 2 :, 2 /, s, z) = 1. And we get The 

remaining ambiguity which leaves Wt 2 and invariant 
should satisfy the condition: V(x,y, s,i) = V{x,0,s,i) 
and ^{x,y,s) = $(a;, 0, s). Any ejb transformation will 
change Wra, so Ets is fixed to 1. 

Similarly, for Ti transformation, using remaining 
gauge transformation, we have 

Wti ( 2 :, y, s,i) {x, 0, s, i)WTi (x, y,s,i)V~^{x - 1,0, s, i) 
0Ti{2:,22, s) ^ 0 Ti{ 2 :, 22 , s) 4’(2:, 0, s)$*(a; - 1,0, s) 

(B6) 

Thus we can set Wt^ {x, 0, s, z) = I and 0 Ti (x, 0, s) = 0. 
Then Sp, ( x, y, s, i) = eTi{x,0, s,i) = 1. Further, accord¬ 
ing to Eq.( |Bl[ ), site tensors are translational invariant 
under this gauge: 

_ ^(x,0,s) _ rps ^(0,0,s) fB7') 


To keep this property, the allowed transformations are 
only sublattice dependent: V{x,y,s,i) = V{s,i) and 
'^ix,y,s) = $(s). 

In the gauge we choose above, we can solve Eq.(B2) as 


Wti ( 2 :, y, s, i) = 77^2 
Wt 2 ( 2 :, 2 /, s, z) = I 
0Ti {x, y, s) = 

0T2 ( 2 :, y, s) = 1 (B8) 


Now, let us add Gq rotation symmetry. Under Cg sym¬ 
metry, tensors will transform as 


^{xysi\x'y's'i') ^ ^[xysi\x'y's'i') 
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where 


as well as 




(BIO) 


f \ xy+},x(x+l)+x+y ^ ^ 

©Ce (a;> y> w) = Mi 2 ©CeW 

r\ f \ xy+\x(x+l)+x+y^ , ■. 

©Ce (a;, y, u) =/ii 2 '^c^v) 

r\ f \ xy+^x(x+l)^ , > 

©Ce (a;, y, u) = /ii 2 ^Cb[W) 


(B16) 


From group 

t—Irr-i—Ir 


relation 




C^^T^^T^CeT 2 = e, we get 


= where we define Wfi{s,i) = W/j(0,0, s,i) and Qr {s) = 
0fl(O,O, s). I nserting the above result back to Eq.(Bll) 
, , and Eq.(B12|, we conclude that all yT,c« = 1- 

Wc^iCe{x,y,s,i))Wj.^^{T 2 CG{x,y,s,i))WcB{T 2 Ce{x,y,s,i)) purtheF^lince = e, we get 

WTACe^T 2 Ce{x,y,s,i)) = T^T,CBXT,CB{x,y,s,i) 

-wui-i/rm f irrn t ■.■,Wce[x,y,s,t)Wce{Ce ix,y, s,i))Wcb{Cq {x,y,s,i)) 

^Ce i<^6{x,y,s,i))Wj,^ {T2CG{x,y,s,i))WTAT2Ce{x,y,s,i)) 


Wcb {CgT 2 {x, y, s, i))WT 2 ('72(x, y, s, i)) = Vt^CbXt^Cb {x, y, s, i 

(Bll) 

as well as 

©Cs (<^6(2;, y, s))©r^ {T 2 Cg{x, y, s))Qcb {T 2 Cg{x, y, s)) 
©Ti(C'6'^r2C'6{x,y,s)) = MTiCs 


,^06(^6 ^(a;,y,s,i))tFc6(Cg ^(a;,y,s,i))VFc-s(C'g ^{x,y,s,i)) 
= VCBXCB{x,y,s,i) (B17) 


Using Eq.(B151 and Eq.(B17), we can simplify the above 
equation as 

WCb {w, i)wcB (I'l (u, i)wcB (w, iCe ) ’ 
WCe(.X,iCe)wCe{u,ic<i) = m 2 Vc<iXCe{x,y, s,i) (B18) 


©C6(C'6(a;,y,s))©T2(72C'6(a;,y,s))0Ti(T2C6(a;,y,s)) So, we conclude that xce{i) = XCB{0,0,u,i) = 

Irr. ( trr t TT ♦ . .^Ce(a^, y, s,*) and xce(0 = XCe(^Ce)- Under remaining 

^Cg\C^T 2 {x^ y, 5 ))©t 2 (^2(^5 y^'V) — ^t^Cq _[ J_ Xt 2 Cq\^^ '^€Cq{sj i), xCq (i) changes as 


(B12) 

Due to the ry-ambiguity, we can always redefine Wr —> 
yWR and Qr —>■ y0ij, which has no physics consequence. 
Thus, by redefining 

Wti yTiCeWri, Wt2 ^iTiCeyTaCeW^Ta, 

0Ti yT2C6©Ti, 0T2 yTiCsMTaCe) (B13) 


we set the right side of Eq.(Bll) and Eq.(B12) to be I 
and 1. Then, by performing transformation 

U(a:, y, s, i) = 


{x+y) 


®(a^,y.a) =MTiC6/^T2Ce 


(B14) 


Wt; and are changed back to its original value in 

Eq.‘@. 

Using scb ambiguity, we are able to set 
XTiCBix,y,s,i) = 1 and xtjCsCO, y, s, 0/6) = 1. The 
remaining ecg is independent of unit cell coordinate, 
namely ecg(a;, y, s^_7)= eceis,i). Then, by solving 
Eq.(Bll I and Eq.(B12), we get 

WcB{x,y,u,i) = 

Wce (a;, y, u, i) = 

VFc,(a;,y,zc,i) = 

(B15) 


XCsW XC6(*)n^C’6(s,*)ec’6(suC6) (B19) 

S 

By choosing proper ecg, we can set xCg = 1- 

Then, by performing unit cell independent gauge trans¬ 
formation V{s,i), Wce transforms as 

Wce{x,y,u,i) V{u,i)Wce{x,y,u,i)V~'^{w,icB) 
WcB{x,y,v,i) V{v,i)WcB{x,y,v,i)V~'^{u,i) 
WcBix,y,w,i) -)■ V{w,i)WcBix,y,w,i)V~'^{v,i) 

(B20) 

Thus, we can set wce(v,i) = wce{w,i) = wce{u,a) = 
wce (at, c) = I by choosing proper gauge. Now, ecg is fixed 
to 1. And the remaining U-a mbig uity satisfies V{s,i) = 
V(i) = ViicB.)- We solve Eq.(B18|) under this gauge as 


wcsiu^b) = wcB{u,d) = m 2 VCe 


(B21) 


Similarly, for phase factor, the corresponding equation 
reads 

©Ce (a;, y, s)Qcb {Cg^{x, y, s))Qcb {Cg^{x, y, s)) 

©Ce {Ce^ix, y, s))0c6 [CQ^ix, y, s))©C6 (Cg"® (x, y, s)) 

= yce (B22) 


From Eq.(B16), we conclude 

^Ceiu)^Ceix)QCe{w) = ±(/ii2/iC6)^ 


(B23) 
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Due to 77 -ambiguity for Cq, we can always redefine sym¬ 
metry transformation rules as 


Wcg J • VDCe? ©Ce -J' -Qce 


(B24) 


to absorb the minus sign in Eq.(B23). Further, Wcg can 
be transformed back to the original form by performing 
gauge transformation 


V{u/w,a/c) = V{v,b/d) = I, 

V{u/w,b/d) = V{v,a/c) = r]. (B25) 


Inserting the solution back to Eq.(B31), we get all 

XcrT2{x,y,S,i) = 1. 

Similarly, the equations for phase factor reads 


y, s))0^2 (T2cr(a;, y, s))0^(T2cr(a:, y, s)) 

0Ti(o'“^T2Cr(a;,7/, S)) = flaTi 

y, s))Qti y, s)) 0 a(Ticr(a:, y, s)) 

0T2 {( 7 ~^Tia{x, y, s)) = 

(B33) 


The remaining E-ambiguity satisfies V(x, y, s, i) 

V{t) = ViicJ. 

Under phase transformation $(s), we get 


By solving the above equations, we get 

0„ (x, y, s) = 0,, (s) (B34) 


0(s) ^ 0 (s)$(s)4>(C'6-1(s)) 


(B26) The equation corresponding to cr^ = e reads 


So, we can set 0(u) = 0(u>) = 1. Now, ^-ambiguity is 
only le ft wit h an overall phase factor. Further, according 
to Eq.(B23l (without minus sign), we have 


Wa{x, y, s, i)W„{a{x, y, s, i)) = yaXaix, y, s, i) (B35) 


Combine Eq. (B32) and Eq. (B35), we have 


^Ceiu) = {yi 2 yCe)^ 


(B27) 


Notice that in the above gauge, according to Eq.(B9), 
all site tensor are equal, namely, 


rj^U _ rj-yV _ rpU 


(B28) 


Now, let us add reflection. For reflection symmetry tr, 
we have 

rj^ix.y.s) ^ 

^{xysi\x'y's'i') ~ ^ ^(xysi\x'y's'i') (B29) 

where 

a o 
a o 

(B30) 

From group relation a~^T^^aT 2 = e and 

a~^T 2 ^crTi = e, we can list the corresponding equations 
as 


y, s, *))VEj:^(T 2 cr(a;, y, s, i))IT<^(T 2 cr(a;, y, s, t)) 
WTi(cr"^r20'(a;,y,s,'j)) = y^TtX<yTt{x,y,s,i) 

W~'^{a{x, y, s, 7))lUf/(Ticr(a;, y, s, i))IT<^(Ticr(a;, y, s, t)) 
Wt 2 {(y~^Tia{x, y, s, i)) = y„T 2 XcrT 2 {x, y, s, i) 

(B31) 

Using Ea ambiguity, we are able to set Xo-Ti ix,y,s,i) = 1 
and XcrT 2 (0, y, s, c/d) = 1, with remaining £a{x,y,s,i) = 
£a{s,i)- 

Then, we can solve the above equation as 

IU<,(x,y,s,i) = vlTi'naT 2 'ni 2 '>Xa{s,i) (B32) 


{yaTiVaT 2 T^^Wa{s,i)wa{cr{s,i)) = y^Xa{x,y, s,i) 

(B36) 


So, we conclude that rjaTi 
Xa{s,i). By applying cr on 


= d<TT 2 , and Xcr{x ,y,s i) 
both sides of Eq.( B36| ), 


get 


we 


Wcr{cr{s, i))'Wa{s, l) = yaXa{(j{s, i)) 


(B37) 


Since the left side of Eq.(B36) and Eq.( B37| ) are equal, 
we have X(y{s,i) = X(t(o'(s, i)). In particular, we have 


Xa{w,a/b) = Xa{v,c/d) = Xa{w,a/b) 


(B38) 


which means 


Xa{w,a/b)=±l (B39) 

Using the remaining £a-{s,i), we can set Xa-iu,a/b) = 
Xcr(u, a/ 6 ) = 1. Then, we are left with £cr{s,i) satisfying 
the following relations: 

£ry{u,a) = £a{v,b), £^{u,b) = £„{v, a) (B40) 

For group relation cr^^CecrCe = e, the corresponding 
equation is 

W~^{(jix, y, s, i))Wce {(^{x, y, s, z))IU<^(C'e'V(a:, y, s, i)) 
Wcf,{(r{x,y,s,i)) =VaCeXaCeix,y,s,i) (B4I) 

By simplifying the above equation, we have rja-Ti = yi 2 , 
and X<yCe{x,y,s,i) = XaCsisJ), with 

Wa(v,i)Wa(v,ia 2 ) = VaVaCeXa{v,i)XaCe(.V,i) 
Wc,{u,i)wCs{u,ial)Wa{w,ia2) = dcrdcrCe dl2XfTC6 (u, *) 

= VaVaCeVl2X'y{w,icr2)X'yC6{w,i^2) (B42) 
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Using remaining e^j transformation, we are able to set 

XcrCe{u/v,a/b) = 1 

X^Ce {w, a) = Xaiw, a) = ±1 

X<yCe{w,b) = Xcr{w,b) = ±1 (B43) 


Then, we get solutions for 0cr(s) as 

0(T {u) = (/io-) ^ 1 

Qaiv) = flCgQCfi(.u)Qa{u); 

Ocriw) = tJ-crfJ-Ce{^Ceiu)(^a{u))~^. (B51) 


Further, by performing the remaining gauge transfor¬ 
mation V{x,y,s,i) = V{i) = V{icg), transforms as 

Wcr{x,y,u,i) V{i)Wa{x,y,u,i)V~^(id) 
Wa(x,y,v/w,i) V(t)Wa(x,y,v/w,i)V~^(ia 2 ) 

(B44) 


Then, we can set w{u, a) = I. The only U-ambiguity left 
is a block diagonal independent of sites and legs matrix 
V(x,y,s,i) = V. 

According to Eq.(B36) and Eq.(B42), we can solve the 
transformation rules for reflection as 


(u, a) = I, (u,b) = XaVaCg , 

Wa{u,c) = XcrVaVcrCei Wa(u,d) = rj^', 

W^{v,a) = yCeVaCg, Wa(v,b) = Xam 2 , 

Wa{v,c) = IJaVCg, Wa(v,d) = Xa77l2?7<T»?<TC6; 

Wa(w,a) = XaVCg, Wa(w,b) = T]i2r]aCe, 

'Wcr{w,c) = rj^yCgVaCe, Wa(w, d) = Xa'ni 2 Va] (B45) 

where = Xa(w,a) = x^(w,b). 

For the phase factor, the equations read 

[Qa(u)f = Qa(v)Q„(w) = (B46) 

as well as 

[0(j(u)] = ya^aCg 

&cr(u)Qa(w)Qcg{'u) = yi2yayaCg (B47) 


where we have used the fa ct Xa(s , i) = Ili Xo-Cg = 1- 
According to Eq.(B27), Eq.(B46) and Eq.( |B47 ), we 
have 


[0ct(u)0,x(w)0C6(m)]^ = (iJ-12tJ-ay<TCe)^ = 1 

= ■ yayaCe ’ M 12 MC 8 = i^l2yCgi^^Cg (B48) 

So, we get ycrCg = Mi 2 MC 6 - In our case, site tensors are Z 2 
odd, then the relation for y causes additional constraint 
for rj: 


VaCe = 'ni 2 llCg 


From Eq.(B46), we get 

0<t(m) = 


(B49) 


(B50) 


However, using the 77 -ambiguity for cr, we can absorb the 
minus sign in above equations by redefining Wa —>■ 
and 0CT —?► /r0CT- Further, by performing gauge transfor¬ 
mation V(a/b) = I, V(c/d) = J, we can transform to 
their original forms. 


Let us summarize the result for lattice symmetry: 

WTi(x,y,s,i) = 77 ^ 2 : 

WT 2 (x,y,s,i) = I 

TTr f ■\ xy+\x(x+l)+x+y , .. 

Wcg (x, y, u, 1 } = 7712 ■'nce (u, *), 

iy+ ^x(x+l)+x+y 


Wceix,y,v,i) = 7712 " 
fUce(U2/,n;,7)=77^2"^2"^"+'^ 

W^(x,y,s,i) = 77i2^^+“*'77;<^(s,i) 


(B52) 


where 


WCe(u,a) = WCeiu,c) = I, 
wcg(u,b) = wcg(u,d) = yi 2 ycg, 


(B53) 


and w„(s^i) are given in Eq.(B45l with additional con¬ 
dition r]„Cg = rii 2 r]Ce- 

For phase factor 0^^, we get 

0 Ti(a;,y,s) = y^^, 

QT 2 {x,y,s) = 1 

r\ / \ xy+^x{x+l)+x+y^ , . 

0 C6(a;,y,u) =/7 i2 &Ce(u), 

„ , . xy+\x(x+l)+x+y 

0C6(a;,y,u) = 77 i2 
Qcg{x,y.w) = yT'^-^^^"'^ 

&aix, 77, S) = /7i+^+“*'0<,(s) 

(B54) 


where 0Cg(u) and 0 ct(s) are given in Eq.(B27) and 
Eq.jB^. 


b. Adding time reversal symmetry 


Now, let us consider time reversal symmetry T■ The 
transformation rule for time reversal symmetry is defined 
on Eq.(lO). We should keep in mind that time reversal 
is antiunitary, so for symmetry operation W^, we have 
TWrT-^ = VF*. 

From group relation T(~^T~^TiT = T 2 ^T~^T 2 T = e, 
we get 


(x, y, s, i)[W^^{x, y, s, i)\*W^^ (x, y, s, i) 
W^T['^(x,y,s,i)) = yTiTXTiT(T^^{x,y,s,i)) 

Wt 2 (a:, y, s, i)[W^^(x, y, s, (x, y, s, i) 

W^T 2 '^(x,y,s,i)) = yT 2 rXT 2 r(T 2 ^ix,y,s,i)) (B55) 
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Similar to previous case, by using ej- transformation, we 
are able to set XTiT xt^T to be identity. The solution 
for the above equation is 

Wt{x, y, s, i) = VTiTVT 2 r‘^r{s, i) (B56) 
The remaining ej- is independent of unit cell coordinate 

{x,y)- 

For group relation <j~^T~^cfT = e, we have 

ivT,TVT 2 rT^^w-^is, *)]* 

w*{s, i)w^{a{s, i)) = riaTXaTi^^ix, y, s, i)) (B57) 

So, we conclude x<yTix,y, s,i) = XcrT('S,*) and t]TiT = 
riT 2 T- Inserting the solution for Wa{s, i) into above equa¬ 
tion, we get 

i)wr{cr{s, i)) = *)) (B58) 

Acting a on both side of above equation, we get 

wr^(o-(s,i))wr(s,i) = y^rxW(s,i) (B59) 

Since the left side of above two equations are hermitian 
conjugate to each other, we conclude that Xcrr(s, *) = 

Let us consider Cq^T~^CqT = e. The corresponding 
equation is 

I^Ce {x, y, s, i) [VFf ^(a;, y, s, *)]* (x, y, s, i) 

W^CQ'^{x,y,s,i)) = VCeTXCeTiCQ^{x,y,s,i)) (B60) 

Then, we get rjT^r = VT 2 r = I, and xCeT{x,y, s,i) = 
XCaT^s^i). Inserting values of Wcg, we have 


So once we can determine the value of wj- = W'j-{u, a), 
we get the complete the solution of time reversal symme¬ 
try Wj- with W-j-{x,y, s,i) = W 7 -(s,i). And W'f{s,i) can 
be expressed by wj- as follows: 


Wq-fyU, a) = Wq-, 
wriu, c) = VaTVCeTWT, 
WT{v,a) = rjCeTWT, 

Wr{v, c) = VaTWT, 

Wq-iw, a) = Wq-, 
Wq-(W, C) = VaTVCeTWq-, 


wt{u, b) = r]Cerwq-, 
wriu, d) = ri^TWT', 

WT{v,h) = wr, 

WT{v,d) = llaTVCeTWr; 
wr{w, b) = rjCeTWT, 
wr{w, d) = r]crTWT\ (B67) 


Now, let US determine © 7 -. The equations for © 7 - read 
as 


'd*g{x,y,s)er{x,y,s)e*g{x,y,s)e^{g ^{x,y,s)) = figr 

(B68) 

where g labels lattice symmetry generators Ti, T 2 , Cg, cr. 
Further, under action of global phase factor $, 07- 
changes as 

^rix, y, s) 0r(a:, y, (B69) 


Therefore, by choosing proper phase, we can always set 
Qq-{u) = 1. Combine with Eq.( |B54 ), we are able to solve 
Eq.(|B68). The solution is Qq-{x,y,s) = <dq-{s), where 


Qt{u) = 1, 

07-(w) = /il2/rC6, 

Qq-{w) = 1. 

(B70) 


Wj-^{s,i)wr{CQ ^{s,i)) = VCeTXCeTi^G ^(^>0) (B6I) And the constraint on y,gr is 

Under transformation Wq —^ eq-Wq-, xCeT changes as Mo-r = y<y, yCsT = M 12 MC 6 


(B71) 


XCer{s,i) xc^ris,i)£r{Cfi{s,i))£r{s,i) (B62) 

So, we can set all XCeT{s,'i-) = I; with remaining £ 7 - = 
£ 7 -(s, a/b) = £^{s, c/d). The above equation is simplified 
as 

Wr^{s,i)wriCQ^{s,i)) = riCgT (B63) 


Since site tensors are Z 2 odd, we also have constraint on 
VgT as 

VaT = Va, VCeT = Vl2VC<i (B72) 

Finally, let us consider group relation = e. For Wq-, 
we get 


Let’s try to fix x<tt{s, *) by remaining £ 7 -. We observe 
that XcrT{s,a/b) — XcrT{s,a/b){£^)‘^. So, we can set 
Xa-Tiu,a) = I. Further, we get 

Xariv, b) = xlq-{u, d) = x^riu, a) (B64) 

due to the relation x<tt{s,i) = XCTr('^('®’*))• have 

w^^(u,a)wq-{u,d) = TjuT (B65) 

From Eq.(B63) and Eq.(B 6 ^, we conclude 


wt{s, i)w^{s, i) = yrxrix, y, s, i) 


(B73) 


where we use th e fact that Wq-(x,y, s,i) = wq-{s,i). In¬ 
serting Eq. (B67) back to the above equation, we conclude 
that XT = Xt{x, y, s, i) = ±1. So, we have 


WTis,i)w^{s,i) = rirXT 
Similarly, for phase factor Qq-, we have 
07 -(s,z) 07 (s,z) = gr 


(B74) 

(B75) 


Wq-^{s,i)wq-{a{s,i)) = T]aT 
Namely, we have Xo-Tis,i) = 1. 


(B 66 ) 


We conclude that yq- = 1. 

In our case, a physical leg lives Kramer doublets. 
Namely, we have L 7 - = —I, where T = Uq-K is the 
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action of time reversal operator on a physical leg, and 
K is the complex conjugation. In the following, we will 
prove that 777 - cannot be trivial, and site tensor must be 
Z 2 odd tensor. 

To see this, we can act time reversal twice on site tensor 
T®. Then, we get 


T® = QrWrTQrWrT o 

= n ° T 

i 

= -77r O (B76) 


where we use Eq.(B73) and Eq.(B75l to get the second 
line. And in the third line, we use the fact that = 1 , 
Xr[s,i) = XT = ill while o T® = —T®. So, we 
conclude 777 -oT® = —T®, namely, rtn- must be nontrivial, 
and r® is Z 2 odd. 

Let us try to understand the physical meaning of 
the above statement. We want to construct a PEPS 
wavefunction on the kagome lattice with every site as a 
Kramer doublet, and preserving all lattice symmetries as 
well as time reversal symmetry. However, the above proof 
tells us that, we are forced to introduce Z 2 gauge struc¬ 
ture due to the nontrivialness of 777 -. The Z 2 gauge struc¬ 
ture leads to either a spin liquid phase or a symmetry 
breaking phase in the long range physics. In other words, 
we can never be able to write a trivial symmetric wave- 
function with Kramer doublets on physical sites with no 
ground state degeneracy in this formulation! The above 
argument can be viewed as manifestation of Hastings- 
Oshikawa-Lieb-Schultz-Mattis theorem on PEPS. 


c. Adding spin rotation symmetry 

At last, let us consider the spin rotation symmetry. 
The action of a group element of the spin rotation sym¬ 
metry on site tensors is defined as 

Uen o r® = (B77) 

where S labels physical spins. In our case, the system is 
formed by half-integer spins. For PEPS invariant under 
spin rotation symmetry, we have 

= QeuWenUeu o T® 

Bb = Wen o Be (B78) 

Here, Wen and 0en are projective representations of spin 
rotation symmetry. To see this, let us consider the group 
multiplication relation: 

Ue ifJi ■ Ue^n-i — Ue^n^ (B79) 

The corresponding equation on virtual legs reads 

(a;, y, S, i) ■ We^n 2 {x, y, s, i) = 

XOirii ,02^2 VOirii ,02 ^2 ■We3n3ix,y,s,i) (B80) 


The above equation implies Wen{x^ y, s, i) form a projec¬ 
tive representation of SU( 2 ) symmetry, with coefficient 
U(I) X Z2. However, it is well known that SU{ 2 ) group 
has no nontrivial projective symmetry. Thus, we can al¬ 
ways set x’s and 77’s to be trivial ones by group extension 
ambiguities. Further, we have Qen always equal to I, 
since SU ( 2 ) has no nontrivial ID representation. 

The on-site spin rotation symmetry are commute with 
all lattice symmetry. Namely, we have g~^Ug}gUen = 
e, where g = Ti,T2,Cq,(j. So, for Wg and Wen- The 
corresponding equations for Wn are: 

77, s, i)Wg} {x, y, s, i)Wg{x, 77, s, i) 
Wen{g~^{x, 77, s, i)) = 'qg^enXgfin{x, 77, s, i) (B8I) 

According to the above solution, we have Wg(x^ 77, s, i) = 
I/J, which always commute with Wen- So, we get 

V’ «> i)Wen{g~^{x, 77, s, i)) = Vg,enXg,en{x, y, s, i) 

(B 82 ) 

One can prove that Wg^{x,y,s,i)Wenig~^{x,y,s,i)) 
form a ID representation of SU{ 2 ) symmetry. Thus, 
r]g^en and XgfiR can be set to trivial. We have 

Wen{x, 77, s, i) = wen (B 83 ) 

Namely, the representation of spin rotation symmetry 
shares the same form on virtual legs. 

Now, we consider the relation: 

We= 27 vix, 77, s, i) = rie=2rXe=2r{x, 77, s, i) (B 84 ) 

Using the fact that W6i=2ir(a:, 27 , s, *) = we=2-K, which is 
independent of sites and virtual legs, we conclude that 
Xe=2r{x, 77, s, i) = Xe= 2 r = ± 1 - 

At last, we have the relation Ug}T~^UenT = e, which 
is the equivalent way to say that a spin reverses its di¬ 
rection under time reversal symmetry. Then, for trans¬ 
formation rules on virtual legs, we get 

Wg}{x, 77, s, ^(x, 77, s, i)]*Wenix, y, s, i) 

w^{x, 77, s, i) = T]T, 9 nXT, 9 n{x, 27 , s, i) (B 85 ) 

We can easily conclude that ?7r,en = I, and 
XT,9n{x,y,s,i) = 1 . The result is reasonable, since we 
also expect spins living on virtual legs reverses direction 
under time reversal action. Then, we have 

Wen\'^r^(sB)]*w*efiWr{s,i) = I (B86) 

In our case, physical legs are spin doublets, so Ug=2-K = 
—I. Using similar proof as in Kramer doublet case, we are 
able to show that r]e=2Tv must be nontrivial and site ten¬ 
sors must be Z2 odd. So, we conclude rie=2-K = ??r = J- 
Further, it is straightforward to see that one can always 
redefine J such that X9=2-n = 1 - WLOG, we assume 
J = Idi © (-lua)) where Di + D2 = D is dimension of a 
virtual leg. 
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Now, let us fix the form of wgfi and wj- using the re¬ 
maining y-ambiguity, which is leg-independent block di¬ 
agonal matrix. Under gauge transformation U, we have 

Weft VwefiV~^ 

wr Vwr[V-'^]* (B87) 

So, first, we are able to set the (reducible) representation 
of spin rotation symmetry on virtual legs as 


matrix under consimilarity has already been studied in 
the mathematical literature Ref llOBl In the following, we 
will give results for the two cases XT = ± 1 - 

First, consider xt = 1- Then the rii dimensional extra 
degeneracy space should accommodate representation of 
Kramer singlet. By choosing proper basis, we are able to 
obtain 

Wr" = I". (B93) 


M 


Wei 


= 0 ( 1 ’! 


AOn-Si 


(B 88 ) 


i=l 


where Si labels spin quantum number while rii is the 
extra degeneracy associated with spin Si. In other words, 
a virtual leg is formed by Ui number of spin Si, where 
i = 1,2,, M. The dimension for spin Si is ni{2Si -I- 
1), and we get the total dimension of a virtual leg D = 
+ !)■ Further, we can arrange the order 
of Si, such that Si is integer (half-integer) for i <= mi 
{i > mi), and we have 5'i < • • • < Smi as well as Sm^+i < 
■■■ < Sm- Apparently, Di = YhA -I- 1) and £>2 = 

Yl,i=mi + 1 ’^i(2>5'i -I- 1). 

After fixing the form of wga, we still left with over¬ 
all gauge transformation V = 0j(Us, (g> 125 ^+ 1 ), where 
Vsi is arbitrary rii dimensional invertible matrix. We 
will use the remaining gauge degree of freedom to fix the 
representation of time reversal symmetry on virtual legs 
wr{s,i). Particularly, let us focus on wr = wr{u,a), 
since representations on other legs can be generated us¬ 
ing Eq.(|B^. 


According to Eq.(B 86 ) and Eq.(B 88 ), time reversal re¬ 


verses spin direction. So the most general form of wr 
reads: 


M 




(B89) 


Z=1 


Then, the remaining gauge transformation must satisfy 
Us. = Us* (B94) 

Namely, Us^ is real matrix. 

Then, we consider the case where XT = Then the 
Ui dimensional space are Kramer doublets, so rii should 
be even number in this case. By choosing proper basis, 
the canonical form for is 

= rim (B95) 

where flm = io’j, 0 Ini/ 2 - After choosing the basis for 
canonical Wr , we are left with gauge transformation U = 
©i(rsi <H) ISi(Si-i-i)) that satisfies 

Us. • fin. = fin. • Us*. (B96) 

Now, let us summarize the result. For a fully symmet¬ 
ric wavefunction on kagome lattice with a half-integer 
spin per site, there are at most 2^ x 2^ = 32 classes in the 
above framework. The different classes are distinguished 
by parameter ry’s and x’s- Here, J7i2, bCej^o- € 
while other 77 ’s are fixed. And X(n XT = ±1- 


2. Construction of PEPS state for different classes 


Here, Wr is Ui-dimensional invertible matrix. Further, 
according to Eq.(B74|, we have 


M 


©O' 

2=1 


[w. 


r ■ 


^ -i7r(Sf 


1 ) = 'nrXT (B90) 


where = I(—I) for integer (half-integer) 

spin. Note that we focus on the case where r]r = J, so 
we have 


wr'{wr'T = xr (B91) 

Then, under gauge transformation U = 0j(Us, ® 

^2Si + l)^ 

w^' -r (B92) 

If two matrices are related by above transformation, then 
they are consimilar to each other. The canonical form for 


In this part, we will use the symmetry transforma¬ 
tion rules obtained above as constraint to determine the 
Hilbert space for a symmetric PEPS state for all classes. 
We will present the general framework first, and then 
work out possible forms of bond tensors and site tensors 
separately. 


a. Genenal framework 

In the following, we will set up the framework to get 
the constraint Hilbert space for symmetric PEPS. 

Let us start by reviewing Hilbert space of PEPS with¬ 
out any symmetry. We require every virtual leg of site 
tensors is isomorphism to a H dimensional Hilbert space 
V. Virtual legs of bond tensors are all isomorphism to V, 
which is the dual space of V. Further, every physical leg 
is isomorphism a d-dimensional Hilbert space U. 
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The Hilbert space of a single bond tensor V b and a site 
tensor Yt have the following tensor product structure: 


Yb=Yi^Y 

(B97) 

Then, let us add the spin rotation symmetry. For a 
physical leg, there lives a half-integer spin Sq , so we have 


U = Vso 


(B98) 


where Yso accommodates a irreducible representation of 
SU(2) with dimension 2So + 1. 

Let us consider virtual legs of site tensors. Notice that 
all virtual legs of site tensors are related by lattice group 
transformations. So, in the presence of lattice symme¬ 
tries, all virtual legs of site tensors share the same repre¬ 
sentation (can be reducible) for spin rotation up to iso¬ 
morphism. In particular, as we show in the previous sub¬ 
section, we can always choose a proper basis, such that 
the spin operators have the same form on all virtual legs 
of site tensors, as in Eq.(B83). Further, we can decom¬ 
pose V = Vi © V 2 , where Vi with dimension Di denotes 
integer spin representations, and V 2 (V 2 ) with dimension 
D 2 accommodates half integer spin representations. 

We can further decompose V to SU (2) irreducible rep¬ 
resentations as 


M 

V^0(Ds. ©VsJ (B99) 

i=l 

Here D 5 . is an rii dimensional space that labels the extra 
degeneracy for spin-5'i. According to the decomposition, 
the orthonormal basis of V can be chosen as 

l^i, ta, nifs) = |S'i, ta) ® \Si, mp) (BlOO) 

where \Si,mp) S YSi labels an eigenstate of S'^ and S'z, 
while \Si,ta) G Di labels basis in the extra degenerate 
space. Under this basis , the spin rotatio operator shares 
the form as in Eq.( |B 88 ). 

Similarly, virtual legs of bond tensors can be decom¬ 
posed as 


M 


— 0(Dsi 


i=l 


with basis as 


'^p \ — ^ 


(BlOl) 


(B102) 


We point out that the nontrivial Z 2 IGG element is in 
fact 27r spin rotation on all virtual legs. In other words, 
we get Z 2 group {I, g} with a trivial representation on 
Vi (Vi) and a nontrivial representation on V 2 (V 2 ). 

After e stabli shing the structure of a single virtual leg 
as in Eq.(B99), we are able to get the structure of Vg 


and Yt according to Eq.(B97|. Yb and Yt are formed 


by tensor product of SU{2) representation, which can be 
decomposed to direct sums of irreps of SU (2) by Clebsch- 
Gordan coefficients. Further, we require bond states and 
site states to be spin singlets, which gives extra constraint 
to the possible Hilbert space. 

Now, let us add lattice symmetries. There are two 
kinds of constraint caused by lattice symmetries. First, 
a lattice symmetry may act as a linear mapping between 
the Hilbert space of different bonds and sites. In other 
words, if a single bond/site tensor is fixed, one can use 
transformation rules of lattice symmetries to generate 
other symmetry related tensors. Second, a lattice sym¬ 
metry may also be a self-mapping (automorphism) on 
the Hilbert space of a single bond/site. In this case, the 
possible Hilbert space of a single bond/site tensor will 
be further constraint by lattice symmetry transformation 
rules. 

At last, due to time reversal symmetry, we require all 
tensors to be Kramer singlets. 

In the following, we will apply the method developed 
above to solve the possible Hilbert space for the symmet¬ 
ric PEPS wavefunction of all classes. 


b. Constraint on bond tensors 

Let us consider bond tensors first. A bond tensor can 
be viewed as a matrix with dimension D x D. Let us 
define Bi, = B(^^ysi\x'y's'i') as the bond tensor connecting 
two virtual legs {x,y,s,i) and {x',y', s',i')- It is obvious 

that B(^xysi\x'y's'i') = ^{x'y's'i'\xysi)- 

Under the action of spin rotation symmetry vogn = 
0^1 (I^i G) e^®"''®’), the bond tensor is a spin singlet in 
the sense 

Bb = WQfi ■ Bb ■ Wgfi (BIOS) 

Then, we can explicitly write bond tensor By as a block 
diagonal matrix according to the spin quantum number 
of virtual legs as 

M 

= (B104) 

where B^' is an rii dimensional matrix, and Ks^ 
is a (25^ + 1) dimensional matrix labeling singlet 

state. More precisely, the quantum state Ks^ = 

{Si,ma, Si,mp\(Ks-)ap is a singlet state under Stot = 
S'®I + I©S'. Here, TOq = —Si + a—1 labels the quantum 
number of Sf. Namely, we have 

{S,m^,S,mp\{Ks)c.pSl, = 0 (BIOS) 

Using Clebsch-Gordan(CG) coefficients, we get 

{Ks)a.p = (B106) 

where (ps is an indefinite phase. We can absorb the phase 
factor to Bg , thus Ks is always real. For example, we 
have Ks=o = 1, = icfy. 
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From another point of view, the Hilbert space of a 
bond tensor Vs is the tensor produ ct of t wo virtual legs 
V, where V is decomposed as Eq.(BlOl). So, we can 
decompose Vs as 


^ ((Ds, < 8 ) Ds^.) 0 (Vsi ( 8 ) '^Sj)) 
i,3 

= 0 (Ds. 8 Ds, 8 8 Vs,) 

i,j,k 

(B107) 

where V|'‘g. is the “fusion space”, which means different 
ways to fuse spin Si and Sj to spin Sk- According to 
representation theory of SU(2), V|'‘g. is isomorphic to C 
if \Si — S'jl < Sk < Si + Sj. Otherwise, vanishes. 

Since we only focus on spin singlet bond states Sk = 0, so 
we conclude the possible Hilbert space of a bond tensor 
should be 


can easily verify that the above result is not limited to 
bond In fact, it is true for all bond tensors. 

Note that for integer (half-integer) spin S, Kg is sym¬ 
metric (antisymmetric). So, we conclude matrix B^' 
must be either symmetric or antisymmetric depending 
on values of Xa- and r]„. In particular, we can write 

— 0 • In,(2Si-|-l) (B113) 

where jiR = 1(—1) for rju = I(J). Since we also have 
Kg. = we conclude 

(B114) 

Finally, let us consider time reversal symmetry. Bond 
tensor should be a Kramer singlet in the sense 

^Isih\s2i2) = BV(siVl) • B(^sih\s2l2) ■ B4)-(S2,Z2) 

(B115) 


v|=° = 0(ID)s. 8 Ds, 8 V|=sO 8 Vs=o) (BIOS) 

i 


By inserting Eq.(B67), we conclude, for any bond Bt, we 
have 


where we use the fact Vf'^^ vanishes if Si yf Sj. Then 

' ' i..Si _ _ 

Bb e V|=o can be decomposed to S Dg, 8 Ds^ and 
Ksi G ^sTsi 8 Vs=o- Namely, we have 

i\Oci,CX2 ;/3l ,^2 

(B109) 

The above equation is just another way to express 
Eq.(B104). Notice that we use Bb to denote the quantum 


state associated with matrix (tensor) Bb. 

Let us add the lattice symmetry. Given a single bond 
tensor Bb^, we can generate all other bond tensors by 
using the relation R~^Bbo = R~^WrR o Bbg, where R is 
some lattice symmetry here. The explicitly expression is 

B (^H(^xysi)\R{x'y's'i')) 

^ V)) ' B(^xysi\x'y's'i') ‘ {B{x , y , S , Z )) 

(BllO) 

It is obvious that we can generate all bond tensors if we 
consider the group generated by Ti,T 2 and Cq. 

Eurther, reflection a will provide extra constraint on 
the Hilbert space of a single bond tensor. Let us consider 
B{vd\wb) = B(^00vd\00wb)- It is straightforward to see that 

B{yd\wb) = wl{v, d) ■ • w-^{w, b) (Bill) 


According to Eq.( |B45[ ), we have w,y{v, d) = Xa'ni 2 V<yVcyCfi 
and Wa{w,b) = rii 2 'ncrCe- Then, we get 


B, 


{vd\wb) XcrVo- ' 


(B112) 


Namely, for any block of it is either symmetric 

or antisymmetric, depending on values of Xa- and V<t- One 


B; = r],wr -Bb-w^r (B116) 

Eurther, using Eq.( |B89 l and Eq.(B104), we get 

{B^r = (B117) 

’^here we use the fact that Ks is invariant under time 
reversal symmetry: 


Oti-S" 


Kg ■ (e'"^")‘ = AT. 


(B118) 


Here we have w^' = for XT = 1 while w^' = flm for 
XT = -1- 

So, to summarize, the constraint on the Hilbert space 
of a single bond tensor is determined by parameter Xm 
rja- and XT- We will list the constraint case by case. 

1- XT = 1 

In this case, time reversal symmetry reads as 


M 


Wj- 


= 0 (I'> 




‘) 


2=1 


So, according to Eq.(B117), 


{B^r = 


(B119) 


(B120) 


Eurther, we are left with remaining global gauge 
transformation 


V = 0(Vs; 8 l2Si+i) 


where Vs- is matrix defined on 
transformation V, we get 


Bt 


Vn, ■ 


(B121) 
Under gauge 
(B122) 
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(a) 77 ^ = I, Xa = 1 

— S'" 

In this case, we get is real symmetric for 
integer Si, while real antisymmetric for half¬ 
integer Si- 

(b) rja = I, Xa = -1 

In this case, we get is real antisymmetric 
for integer Si, while real symmetric for half¬ 
integer Si- 

(c) 77 ^ = J, = 1 

In this case, we get B^' is real symmetric 
for integer Si, while imaginary symmetric for 
half-integer Si. 

(d) 77 ^ = 3, Xa = -1 

— S’’ 

In this case, we get B^' is real antisymmetric 
for integer Si, while imaginary antisymmetric 
for half-integer Si. 

By using the remaining gauge transformation V, 
we can set the bond tensor to a maximal entan¬ 
gled states. Namely, if B^' is real symmetric, the 
canonical form is 

=Diag(l,...,l,-I,...,-I) (BI23) 

where the number of ± sign ni± is not fixed. After 
doing this, we still left with gauge transformation 
Vm S 0{ni+) 0 0{ni-). Here, we point out that 
different ni± does not lead to new classes. If we 
do not require bond tensors to be maximal entan¬ 
gled, different ni± can be connected adiabatically 
by continuously tuning the entries of bond tensors. 

If B^' is real antisymmetric, the canonical form is 


remaining gauge transformation is block diagonal 
matrix V, which reads 

^ = ®l2S.+i) (BI27) 

i 

where Vs- satisfies 

Vs^ ■ • Fs* (BI28) 

Under the gauge transformation V, we have 

Bh' ^ K. • b!' ■ Vl (BI29) 

(a) 77a = I, x-T = 1 

In this case, we get Bj, ‘ is quaternion sym¬ 
metric for integer Si, while quaternion anti¬ 
symmetric for half-integer Si. 

(b) ?7a = I, Xa = -1 

In this case, we get B^‘ is quaternion antisym¬ 
metric for integer Si, while quaternion sym¬ 
metric for half-integer Si. 

(c) 77a = J, X<T = 1 

In this case, we get B^‘ is quaternion sym¬ 
metric for integer Si, while iB^’ is quaternion 
symmetric for half-integer Si. 

(d) 77a = J, X<T = -1 

In this case, we get B^‘ is quaternion antisym¬ 
metric for integer Si, while iB^’ is quaternion 
antisymmetric for half-integer Si. 

c. Constraint on site tensors 


B^ = U„./2 (BI24) 

where U = iay <S>lni/ 2 - The remaining gauge trans¬ 
formation satisfies U „./2 = • U „./2 • V^.. 

If By is imaginary, the canonical form is similar 
as the real case, except that all entries are replaced 
by ±i. 

2. XT = -1 

In this case, the time reversal symmetry reads as 

M 

Tur = 0(U„, (BI25) 

So, we get the constraint on the bond tensor to be 

(Bf-)* = ■ Bf • (BI26) 

Then, depending on values of /J,a and Si, either 
By' or iBy' is quaternion matrix. In this case, the 


T he H ilbert space of a site tensor Vt is defined in 
Eq.(B97). In the presence of SU{2) symmetry, Vt can 
be decomposed as 

VT=U(g)V(g)V(g)V(g)V 

Vs, (g) Vs,,, (g Vs, (g Vs,,,) 


- 0 


(BI30) 


where 

Dsi„Si,,Si„Si,, — (gDsi,, <gDsi„ <gDsi,^ (BI3I) 

labels the extra degenerate space associated with spins 
on four virtual legs. The basis of 
Ds-„Si,,Si„Si,, is labeled as 

\Si^,ta) g IB4V/3) g \Si,,t.y) g \Si^,ts) (BI32) 
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^So% Si Si Si fusion space, which denotes differ¬ 

ent ways to fuse spin Sq, Si^, 5'^, Si^, Si^ to spin Sk- The 
complicated fusion rules with six spins can be obtained 
by the fusion rules with only three spins as 




0 ^: 

q;,/S,7 


ios., 




(g)V|'‘ 


S^Si, 


(B133) 


Since site tensors are SU{2) singlet, we should focus on 
the Hilbert space with S = 0: 


vf ° ^ i^Si^Si^Si^Sij ^"^SoSi^Si^Si^Sij ^"^S=o) 

(B134) 

The basis for space g. g. g. 0 Vs=o can be ex¬ 

pressed as 

KsoSi,Si,S.,S,^ = iKsoSi,s,,s,,s,Jipj5>^ 

\So,mj) 0 \Si^,ma) 0 ISi^.rUfs) 0 0 \Si^,ms) 

(B135) 


where Kg^g, g. g. g. labels orthogonal singlet states for 
different 1 . 

Then, in terms of the tensor representation, we decom¬ 
pose the site tensor T‘^ as 

® ^'^kaSi,,Si^Si^® ^SoSi^s,^Si^SiJ (B136) 

'^a 

where the state 


We figure out lattice symmetries that maps site tensor 
T“ to itself as follows 

T 1 T 2 CI O = Ce O 

(B138) 

Besides, combining reflection a and rotation Ce, we get 
uCe o 

(B139) 

In the following, we will solve the constraint from above 
symmetry operations. Further, we can prove that the 
whole site tensor can be generated by lattice symmetries 
once we fix quantum states in the Hilbert space satisfying 
the two situations below: 

1 . Si^ is a half integer spin, while other three are in¬ 
teger spins; 

2 . Si^ is an integer spin, while other three are half 
integer spins. 

To see this, let us first consider reflection symmetry a. 
Under the action of tr, for the decomposed parts of the 
site tensor, we have 

^'^k^S,^Si^S,Jal3'yS = {Tgi^Si^Si^SiJS'YPa 
i^SoSi^Si^Si^SiJaPjS = (^SoSi^Si^Si^SiJs'ylSa 

(B140) 

It is obvious that we can choose KsgSi^Si^Si^Si^ to be ei¬ 
ther symmetric or antisymmetric under the permutation 
of 


'y i'^Si^Sii^Si^Sij)al3~/s\Si^,ta) (El \Si^,tp) 

CX,( 3,^,6 




(B137) 


is an arbitrary state lives in the extra degenerate space 

Due to the representation theory of SU{2), 
KsgSi^Si^Si^Sij^ does not vanish only if there are 
even number of half integer spins for Sq, , St^. Since 
^0 is a half-integer spin, we conclude that there should 
always be odd number of half integer spins living on 
virtual legs. For a site tensor, there are four virtual legs, 
so we get two different cases: 

1. Only one virtual leg are a half integer spin, while 
other three are integer spins; 


2. Three virtual legs are half integer spins, while the 
remaining one is an integer spin. 

We now consider the constraint from lattice symme¬ 
try. Remember that in the presence of translation and 
rotation, we can always choose a gauge such that all site 
tensors share the same form, as shown in Eq.(B7l and 


Eq.(B28). Then, in the following, we only need to focus 


on a single site tensor. 


- MKsgS.^S,,Si^SiJii3^s 
(BUI) 

where P is any permutation. The ± sign depends on the 
definition of K. Particularly, we have 

0-0 (R:SoS.„S.,Si,SiJa/375 = MKSgS.^Si^S.^SiJiiSjS 

(B142) 

Eor the two cases we consider here, Si^ is always 
different from spins of other three virtual legs. So, 
KsgSi^Si^Si^Si^ and KsgSi^Si^Si^Si^ are always indepen¬ 
dent tensors. Thus, we can absorb minus sign in the 
above equation by redefining KsgSi^Si^Si^Si^- Then, we 
get 

(J 0 T^= ^ {a 

)^c 5^d 

(B143) 

Remember that the site tensor is symmetric under cr, so 
we have 

r“ = 0 ,,iu<,o- o r“ 

= ©o-IFcr ^ {a O Tg.^g.^g.^g.^ 0 Kg^g^^g.^g^^g.^ 

?^d 5^ 

(B144) 
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As shown in the last subsection, we always choose the 
basis such that WR{x,y,s) S {I, J} for any lattice sym¬ 
metry R. So, we can always define the action of W^- 
trivially on g. g. g. ■ Namely, we can decompose 

WR{x,y,s,i) as 

Wr{x, y, s, z) = 0 {x, y, s, i) (g) l2Si-Hi) (B145) 


Then, from the above analysis, we conclude 


r: 


Si. Si. Si, Si, 


= O T‘ 


Si„Si,Si Si 


(B146) 


Writing the above equation explicitly, we get 

[^Si^Si^Sii^Si^]al3'y5 =Q<Tiu)[WcT ]aa'[Wa \pp'[Wu 
[Wa '^]<55'[TSi^SijSi.Si^]5'7'/3'a' 

(B147) 


According to Eq.(B45), we have 


w^{u,a)=l, wf (u, 6 ) = x„(^i 2 /rce)^'®, 
wf{u,c) = Xcrifj-uyc^ycrf^, W^{u,d) = 

(B148) 


Then, we can simplify the constraint as 

{"^Si^Si^Si^SiJaP'yS =Q<t{u) ' 

^ (B149) 

Since S'^, S'i. and Si^ are all integer spins or all half¬ 
integer spins, the above equation reads 


integer/half-integer spins, by the above lattice symme¬ 
tries, we are able to generate tensors Tg_ g, g, g, which 
satisfy one virtual leg to be a half-integer/integer spin 
and other three virtual legs to be integer/half-integer 
spins. 

At last, we add time reversal symmetry. The con¬ 
straint of time reversal symmetry reads T® = Qg-Wg-T o 
T®. Since Kg^g. g. g. g. is a Kramer singlet state and 
real, we have 


[K: 




inSy 


]/3/3'X 


]77' [e*’" '•i]ss' l^slsi^ Si Si. Si la'/Sg'i'; 


(B153) 


Then, according to Eq.(B89), the constraint on 


n 


reads 


0/375 = [u;f“(M,a)]a„/[u;./®(u,6)];3/3'X 

(«> c)hY{'^r‘‘ '^)]55'[^S*„S.,S.,S.Jo'/3'7'5' 

(B154) 


where according to Eq.(B67|, we obtain w^{u,i) as 


wf (u, a) = wf, (u, b) = {yi 2 fj.cef^w^, 

wriu, c) = (^i 2 /iCB/r„)^‘^wf, 

(B155) 


And depends on XT- Remember that we will focus 
on the case where Si^, Si^ and S^, must be all integer or 
half integer spins. Then, by inserting Eq.(B155l back to 
Eq.( p^ , we get 


i'^Si Si Si S, ]a/ 3 ’yS — Qa{u)[Tg g g g ]5-iPa (B150) 

a c ‘‘t o. ‘a ■'b ‘•c ‘■d 

where Qa{u) = (na)^ ■ 

We consider the constraint by the rotation symmetry 
now. Similarly, Tg, g, g. g. and Tg. g. g. g_ can be 

obtained from T^. g_ g, g, and Tg. g. g. g. by rotaion 
symmetry: 


"^Si^Si^Si^Si^ — ©Ceff^CeC'e o '^k^Si^Si^Si^ 

Tk.g.,g..s., = Qc,Wc,c, o (bi51) 


By inserting wCf, {u, i) 


defined in Eq.(B53), we get 


_ g_ _ 

Si. Si^] 0/375 = [rC7-“]aQ'['lC7-‘’]/3/3'X 
[Wg-'"]ry^' [Wg-'^\ssi [T'Si.Si,^Si.Si^]o'/3'7'5' 

(B156) 


When XT' = 1, we have then Eq.(B156l is 

simplified as 


rpl _ rpl *■ 

^Si.Si,Si.Si^ - ^Si.Si.Si.Si^ 


(B157) 


So, Tg. g. g. g. is real tensor. 

When XT = ~1) ^ve have Wj-' 
becomes 


= flrii, then Eq.(B156) 


[^Sij Si. Si_j Si Jo/375 — 0 C 6 (M)[ 7 Si.Si,^Si.Si^]/ 3 o 57 

[^Si.S,_j Si. Si Jo/375 = fJ‘ 12 yCfi^C 6 {u)Q<j{u)[Tg.^g.^g.^g.J^SaP 

(B152) 


[^Si.Si,^Si.Si^]o/375 — [f^ni.]aQ'[f^nij^]/3/3'X 
[f^ni.]77' [f^rai^]55' ["Jsi.SijSi.Si^]o',a'7'5' 


(B158) 


d. Examples 

where 0 Cg(u) = (^i 2 MCe)^- 

Thus, once we know tensors Tg.^g.^g.^g.^ with 5^. Now, let us focus on a special case where the physical 
to be a half-integer/integer spin and Si^, Si^, Si^ to be spin S'o = 5 and there are only spin -0 and spin-^ living 
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on virtual legs. Then, to obtain site tensors, according 
the above analysis, we restrict ourselves to two subspace 
of Vf=°: 

1. Si^ = i and Si^ = Si^ = = 0. Then the 

corresponding Hilbert space Ho is 

Ho = Di (g) (Do)^ (g) V|ippp (g) Vo (B159) 

where fusion space Vi i = V^ i with dimension one. 
2 2 2 2 

The singlet state = Kiiqqq is 

(B160) 

since /? = 7 = (5 = lin this case. One can easily verify 
that Kq is invariant under time reversal operator: 

{KDip,, = (ia^),y • {K*)i,p,s (B 161 ) 

Then, quantum state in Ho can be expressed by the 
tensor form as 


To (g Kq 


(B162) 


where Tq denotes an arbitrary quantum state in Di (g) 

(Do)3. 

2. Si^ = 0 and Si^ = Si^ = = ^. The the corre¬ 

sponding Hilbert space Hi is 

Hi = Do (g (Di)^ (g Vi„i 11 (g Vo (B163) 

2 2^222 

Using representation theory of SU{2), we get 



2 


)1 ii 
^222 


=(V°11 (gvf,,(gv°ii)0 
2 2 2 ^ 2 2 

(V °11 (gvf^ (gv\i) 

2 2 2 ^ 2 2 


(B164) 


where Ti, T 2 are tensor representation of arbitrary states 
in Do (g (D 1 )^. 

We can explicitly write down Ki as a quantum state by 
introducing basis | 0 ) for spin -0 and | t); I i) for spin-f: 

^0 =11) o U 000) - U) (g 11 000) 
^i=lt)®(|om)-|out))- 
U)®(iom)-iont)) 
k2 = 11 ) ® ( 2|0 tu) -10 m) -10 ut))+ 
U)g( 2 |o;tt)-|om)-|ont)) 

(B168) 

where we define Ki = (KiY^i^^glmj) (g \mampm~^ms)■ 
Now, let us consider case D = 3 virtual legs V = 00 f. 
Then, There is no extra degeneracy of spins. So, accord¬ 
ing to constraint from bond tensors, only classes with 
= J, Xcr = 1 and XT = 1 can be realized. Other 
classes require even dimensional extra degenerate spaces, 
since bond tensors of those classes are either antisymmet¬ 
ric or symplectic in the extra degenerate spaces. Thus, 
when D = 3, we only left with 7712 , rjc^, and the number 
of classes can be realized is 2^ = 4. 

Given a bond tensor we can fix it as a maximal 
entangled state with the following form 

/±1 0 0 \ 

= 0 0 -i (B169) 

Vo i oy 

Other bonds are all related to by translation and 
rotation symmetry, and can be generated as 

Bji^,)=R-^WRRoBhg (B170) 


So, Vi 1 1 1 has dimension 2 . And we choose the two 
2^222 

basis in Vi 1 11 0 Vq as 
2^222 


i^'^Yap-fS = 


Cl’"" 1 


(B165) 


where a = 1 in this case. Here, ^ = 

(S'iTOiS' 2 m 2 |Jmj) is the CG coefficient. And rrii = 
— S' — 1 0 J is the Sz quantum number. Simialr as the 
previous case, Ki and K 2 are also chosen to be invariant 
under time reversal operator: 

(^1(2) )a/37i = (i0-^)/3/3' (if^^)77' {KY2))il3'X 

(B166) 

Then quantum state in Hi can be expressed by the 
tensor form as 


Ti 0 All 0 r2 g K 2 (B167) 


where R = T"^T^^Cg‘^® with ni,n 2 ,nce € Z- 

For site tensors, they all share the same forr n. Th e 
spin singlet state KsgSi^Si^Si^Si^^ is fixed as Eq.( 


So we can express the site tensor using the quantum state 
representation as 


k ={ko + ki2{pi,p2)} + 0 c 6 (w){a ^ b,c^ d} + Qa{u)- 
{a yy d, & o c} 0 pi2PCe'd>Ce{k^<y{k{^ yy c, 5 yy d} 

(B171) 


where we have 

Ki 2 =a\Kx 0 027^2 

=pi '(It)® |o m) 0 u) 0 |o nt))+ 
P2-(|t)G|0Ut) + U)G|0m))- 
(Pi + P 2 ) • (I t) G |0 tU) 0 U) 0 |0 itt) (B172) 

where we define pi = oi — 02 and p 2 = —ai — a 2 as the 
two tunable parameters. 

Now, let us consider virtual legs V0O0O0|0| 
with D = 6. In this case, there are extra two dimensional 
degeneracy spaces for both spin-0 and spin-f. We believe 
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all 32 classes can be realized in this case. However, here 
we will focus on the 4 classes realized in Z? = 3 case. 

Fixing ?7cr = J, Xo- = 1 and XT = !> the bond tensor 
Hbjj now reads 



Other bonds can be generated by translation and reflec¬ 
tion symmetry as discussed above. 

For the site tensor, we have 

T® ={T'q 0 Kq Ti 0 Ki T 2 O K 2 } 

+ Ocg (u){a -H- b, c <->■ d} + Ocr(u){a <->■ d, b <->■ c} 
+ Mi2MCe(^Ce(u)0^(u){a <->■ c, b <->■ d} (B174) 


where Ti labels a quantum state in extra degenerate 
space, which has dimension 2^ = 16. F urther , the 
tra nsform ation rules of T’s are given in Eq.(B150) and 

Eq.(B152). So, there are three T’s (Tq, Ti and T 2 ) serv¬ 
ing as tunable parameters. Then the tunable parameters 
in T = 6 case should be 16 x 3 — 1 = 47, where the 
additional —1 comes from the fact that the norm of the 
wavefunction has no physical consequence. 


Appendix C: Projective representation, group 
extension and second cohomology 


If uja and iOb are both factor sets, then ujab = 
uJaUJb is also a factor set, where {ujaUJb)(,gi, 92 ) = 

Wq((/i, 52 )^ 6 ( 31 , 52 )- The product of factor sets is associ¬ 
ated with tensor product of projective representations: if 
uja,ujb are factor sets of Fa,rb, respectively, then tUab is 
factor set of the tensor product representation Fq 0 F;,. 

Now, let us define the equivalent class for factor sets. 
Suppose we allow a redefinition of the F’s by 

F'(g) = A(g)F(5) (C4) 

where A defined as A : G —>■ A. This induces transforma¬ 
tion of the factor set: 

^'( 91 , 92 ) = X{9i) ^^Xi92)M9i92)~^i^i9i, 92 ) (C5) 

where lo' is also a factor set. Two factor sets ui and uj' 
are said to be equivalent if they are related by the above 
equation for some A, and we write as w ~ w'. We group 
all equivalent w as a class, and define equivalence class 
by c(a;). Then, one can easily verify that the equivalent 
classes form an Abelian group with product defined by 

c(wi)c(w2) = c(wia;2) (C6) 

The Abelian group of factor set equivalence classes is 
isomorphic to the cohomology group H^{G,A). We can 
view it as definition of the second group cohomology. 
Any factor set w is named as cocycle, which is classi¬ 
fied by Z‘^{G,A), while A is called coboundary, classified 
by B^{G,A). Then, we have 


In this appendix, we will introduce mathematical tools 
for symmetry fractionalization, including projective rep¬ 
resentation, group extension as well as the second coho¬ 
mology. Readers may refer Ref l37l for more details. 

Consider a group G with elements 9 G G. We call r(g) 
a projective representation of G with coefficient A, where 
A is an Abelian group, if 

r(gi)r(g2) = w(gi,g2)r(5ig2) (Cl) 


H^{G,A) = Z^{G,A)/B^{G,A) (C7) 

We point out here, the definition of H‘^{G, A) depends 
on the action of G on A. In mathematics language, A 
is a G-module, which is equivalent to say that Vg S G 
may have nontrivial action on A with the action to be 
automorphism of A, rather than just an Abelian group. 
For example, a trivial module just means A is invariant 
under G: 


Here w is a map, which is defined as w : G x G —>■ A. 
According to associativity of matrix product, we get 

r(ffi)r(52)r(g3) = ^( 51 , 52 )^( 5152 ,5'3)r(5ig253) 

= w(5i,5253)®^w(52,53)r(gi525'3)- 

(C 2 ) 

where appearance of ^\o{g 2 , 93 ) comes from commuta¬ 
tion of gi and uj{g 2 ,g 3 ), which indicates action of G on 
coefficient A may be nontrivial. Further, we require the 
action oiVg G G on Abelian group A should be an auto¬ 
morphism of A. Then the associativity constraint for w 
is 

‘*^(51.52)^(3152,53) =^{91,9293) ^^^{92,93)- (C 3 ) 

Any function uj satisfy the associativity constraint is 
called a factor set. 


% = a, V 3 e G, a G A (C 8 ) 

One should always fix a G-module A, and then classify 
projective representation with coefficient A. However, 
in our case, IGG = Z 2 , the automorphism of Z 2 only 
contains trivial one. 

We now put projective representation aside and turn 
to discussion about group extension. Assume group E 
has a normal subgroup A. Then, we can define G as 
quotient group 

G = E/A (C9) 

with associated homomorphism tt : E ^ G. Then, it is 
natural to define G-module A: given g G G, the action 
of 3 on A is characterized by 

= gag~^ (CIO) 
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where we choose g so that ttQ) = g. E is called an 
extension of a group G by G-module A. In particular, A 
is central in E if and only if the G-action is trivial. In 
this case, the extension is called a central extension. 


Now, let us discuss about the relation between group 
extension and projective representation. Roughly speak¬ 
ing, the equivalent class of projective representation has 
one-to-one correspondance with group extension of G by 
G-module A. Namely, group extension E is also classified 
by 2-cohomology H^{G,A). Projective representation P 
can be viewed as a map T : G ^ E such that ttoP = id^. 
Then, factor set w is naturally induced by P, and auto¬ 
matically satisfies associativity constraint Eq.(C5|. No¬ 
tice that the choice of P is far from unique, and we can 
always redefine P as shown in Eq.(C4). 


In the following, we will develop a general method to 
solve the inequivalent projective representations for dis¬ 
crete group G with G-module A, and the corresponding 
extended group is E. Particularly, we will focus on G 
as the symmetry group of kagome PEPS defined in Ap¬ 
pendix]^ while A = IGG = Z 2 is a trivial G-module. 
Spin-rotation symmetry will also be discussed. 


Let us first set up the general framework. G is defined 
by generators {Ti, T 2 , Gg, cr, T} as well a s th e relation 
between these generators, as shown in Eq.( |A3| ). In other 
words, Vg G G, there is one integer set {ni, 712 ,nce,rio-} 
such that 


follows 

p(T2)-ip(ri)-ip(T2)p(ri) = 7712, 

P(a)-ip(Ti)-'P(a)P(T2)=77,T., 

P(a)-ip(T2)-ip(a)P(Ti)=7;,T,, 

P(G6)-'P(T2)-ip(G6)P(Ti) = 7?c«T,, 

P(G6)-'P(T2)-lp(ri)P(G6)P(r2) = 7?CeT., 

P(a)-ip(G6)P(a)P(G6)=7?,Ca, 

r(G6)'* = T^Ce, 

r(cr)2 = ga, 

nr? = 77 r, 

P(g)-ip(r)-ip(g)P(r) = ggT, V 5 = Ti, 2 ,a,G 6 (CM) 

where we get 13 77 ’s. One may expect the number of 
cohomology classes should be 2^^, however, as we will 
see later, there is lots redundancy. 

Now, let’s try to eliminate those redundant parameters 
by choosing gauge of P(g). By doing gauge transforma¬ 
tion P(Ti) gce,T 2 ^iTi) and T 2 gCeT 2 V<yTiT 2 , we are 
able to set = VrrTi = I- Notice, other 77 ’s may also 

change, however, we can always absorb the change by 
redefining other 77 ’s. 

Then, we have 

P(T2) = P(77)P(ri)P(7T)-l (C15) 


5 = T”iT2"’^Gg‘^V”-' (CIl) 

As discussed before, projective representation can be 
constructed from E by map P. Let us first choose the 
gauge such that P(l) = I. Then, this implies w(l,I) = 
a;(g, I) = ui{l,g) = 1,V(7 G G. Let us consider a particu¬ 
lar relation between generators as 

T-^t-^T 2 Ti = e (C12) 

Lifting this relation to E by P, we have 

p(T 2 )-ip(Ti)-ip(T 2 )r(ri) = 7712 (CIS) 

where 7712 G Z 2 . Similarly, for all relations, we obtain a 
set of 77 ’s. 

These 77 ’s are closely related to classification of pro¬ 
jective representation. However, there are some issues 
arise. First, the 77 ’s are not, in general, one-to-one corre¬ 
spondance with cohomology classes. There may be some 
redundancy in this description. Second, some choices of 
77 ’s may be inconsistent and not give a legitimate factor 
set. To see this, let us solve the classification of pro¬ 
jective representations for the kagome lattice symmetry 
group completely. Conditions for group relations are as 


By applying the above equation, we get 

77,r, = P(a)-ip(ri)-ip(a)P(r2) 

= P(a)-ip(ri)-ip(a)P(a)P(Ti)P(a)-i = 1. 

(C16) 

as well as 

7?T.r = p(T2)-ip(r)-ip(T2)p(r) 

= p(a)p(ri)-ip(a)-ip(r)-'p(a)p(ri)p(a)-ip(r) 
= 77.rr(a)p(Ti)p(r)-ip(Ti)p(a)-ip(r) 

= (77<,r)"r(r)p(Ti)-ip(r)-ip(ri) = gr.r- 

(C17) 

After above calculation, we are left with 9 free tunable 
Z 2 parameters, 

{ 7712 , VCeTi , VaCe , Va, VCe , VT, VTiTVaT, ICeT} (C18) 

So, we expect H'^{G, Z 2 ) = 2®. 

Now, let’s adding spin rotation symmetry Rg{9n). We 
have 

P(R5(27r)) = g0=2-K, 

T{Rg{9n))V{g) = 77g.,«P(5)P(i?,(0n)) (C19) 

where g G Ti, T 2 , Gg, cr, T. As argued in Refl37l one can 
always set = 1. Thus, by including spin rotation 

symmetry, we get an extra parameter ge= 2 Tr- So the num¬ 
ber of cohomology classes becomes 2^°. 
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Appendix D: Distinguishing different classes by 
lattice quantum numbers 

It has been shown that for system defined on a torus, 
lattice quantum numbers can be served as useful tools 
to distinguish different phase^^UllIlIIinil, Here, we will 
show lattice quantum numbers are also very useful to 
distinguish different classes. 

First, let us set the framework to extract quantum 
numbers of symmetric PEPS wavefunction I'k) defined 
on a torus. Similar to the PEPS on a infinite plane, ten¬ 
sors of the torus PEPS wavefunction |'I') satisfy 

^{x,y,b) ~ ^^rR ® R(x,y,b) 

where R is the global symmetry operator, while Wr acts 
on virtual legs as a gauge transformation. Then, the 
global quantum number of R is simply product of all 
Qr. Namely, we get 

i?|4/)= n (D 2 ) 

x^y,s 

Let us focus on two particular classes: 

1. Class-I is labeled by 7712 = rjCg = b Vo- = ^ and 
Xct = XT = 1- Qi = Q 2 spin liquid phas^^ belongs 
to this class. 

2. Class-II is labeled by r]i 2 = I, vCe = Vo = J and 
Xo = XT = 1- Qi = —Q 2 spin liquid phas^ 
belongs to this class. 

To distinguish quantum numbers of these two classes, 
let us consider systems on a torus with (4n -|- 2) unit 
cells with even number of unit cells in Ti direction and 
odd number of unit cells in T 2 direction. Notice that 
wavefunctions defined on this system explicitly break Cq 
rotation symmetry. However, it still preserves inversion 
symmetry R^^ = Cq. In the following, we will show that 
the ground state manifold of these two classes defined on 
systems with (4n -|- 2) unit cells form distinct representa¬ 
tions of symmetry group generated by translation Ti, T 2 
and inversion i?^. 

Let us first list symmetry transformation rules for infi¬ 
nite PEPS. For translation, we can choose a proper gauge 
such that Wti (a:, y, s, i) = IFt 2 (a:, y, s, i) = I as well as 
0Ti(a:,y, s) = 0 T 2 (a:,y, s) = 1 for both two classes. The 
symmetry transformation rule of reflection Rt^ can be 
generated by Cq rotation as 

Wr^ (a:, y, s, i) =Wce (x, y, s, i)Wce y, s, i))- 

^Ce(CQ^(T,y,s,i)) 

(a:, y, s) =Oce (t, y, s)Qce y, s))- 

QcACQ\x,y,s)) (D3) 

Thus, for infinite PEPS in Class-I, we have 


TABLE I. Translation and inversion quantum numbers for 
topological degenerate ground states on (4n -|- 2)-uc lattice 
(even by odd) samples 


(a): Quantum number for Class-I 


Sym. 

1 ^ 0 , 0 ) 


l^'O.vr) 

|^7r,7r) 

Ti 

1 

-1 

1 

-1 

T 2 

1 

1 

1 

1 

Rtt 

1 

-1 

1 

-1 


(b): Quantum number for Class-II 


Sym. 

o' 


l^'O.vr) 

|^7r,7r) 

Ti 

1 

-1 

1 

-1 

Tz 

1 

1 

1 

1 

Rtt 

-1 

1 

-1 

1 


And for Class-II, we have 


WR^ix,y,s,a/c) = I, 
WR^{x,y,s,b/d) = J, 
‘dR^{x,y,s) = i. 


(D5) 


Now, let us turn to PEPS on a torus with (4n-|-2) unit 
cells. For these two classes, one can construct a symmet¬ 
ric wavefunction Idio^o) with the symmetry transforma¬ 
tion rules defined the same as infinite PEPS, since the 
transformation rules of Ti, T 2 and Rt^ on virtual legs of 
both classes are compatible with the system size. Sym¬ 
metry quantum numbers of these states can be calcu¬ 
lated using Eq.( |D2 ), where the result is listed in the first 
columns in Table [I 

As discussed in Sec|V) other bases of ground state man¬ 
ifold is obtained by inserting non-contractible flux loops, 
labeled by |'I',r,o)) and |'I',r, 7 i-)- Wr and Qr change 

their values after the loop insertion, but it is easy to 
extract the quantum numbers of these states. For ex¬ 
ample, let us consider Ti quantum number of |4',r,o)- R 
will move the non-contractible g loop with one lattice 
spacing, leading to a new wavefunction q). One can 
easily figure out that the PEPS wavefunction |4'o_„.) is 
related to llko,,!-) by a Z 2 gauge transformation on the 
column sandwiched by g loops of |'I'o, 7 r) and |'I'(r,g), plus 
the Ti symmetry transformation rule of |'I'o,o)- Since 
there are odd number of sites per column, and site tensors 
are Z 2 odd, the single column Z 2 gauge transformation 
contributes an extra —1 to Ti quantum number. Follow¬ 
ing similar strategy, one can obtain the representation of 
symmetry group on the whole ground state manifold for 
both classes. We list the result in Table HI 

As seen in Table [J ground state manifolds of the 
two classes have distinct representations. So, these two 
classes can be distinguished by lattice quantum numbers. 


Appendix E: Valence bond solid phase 


WR^{x,y,s,i) = I, 
QR„{x,y,s) = 1 . 


(D4) 


As we mentioned before, a single class includes many 
different phases, which can only be distinguished by finite 
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FIG. 15. Visulization of the 12-site valence bond solid state 
as coverings of spin singlets (red thick bonds) on the kagome 
lattice. 
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size scaling. One may ask, is it possible that some par¬ 
ticular phase can be described by different classes? We 
think that the answer is yes, and we will give an example 
in the following. 

Let us focus on Class-I and Class-II discussed in Ap¬ 
pendix]^ As we argued before, the ground state spaces of 
these two corresponding QSL have different lattice quan¬ 
tum numbers on systems with (4n -|- 2) unit cells. 

Now, let us consider the valence bond solid (VBS) 
phases. The VBS pattern can be obtained by Landau- 
Ginzburg theory of visons. The effective Lagrangian of 
visons is constrained by visons’ transformation rules un¬ 
der symmetries, which are calculated in SecjV] It turns 
out that the transformation rules only depend on spinon 
distributions, which are the same for all classes. Thus, 
we expect both Class-I and Class-II give the same VBS 
order pattern after vison condensation. This seems to 
contradict with the quantum number discrepancy men¬ 
tioned before. 

The first observation is that for those samples where 
two classes have different lattice quantum numbers, it 
is impossible to write a compatible VBS order with the 
lattice size. In other words, there are always domain wall 
configurations on those samples. 

To see this, let us consider a particular 12-site VBS or¬ 
der pattern shown in Fig.(15) which is compatible with 
vison symmetry transformation ruleJ^. As shown in Ap¬ 
pendix]^ Class-I and Class-II have different lattice quan¬ 
tum numbers on a torus with (4n-|-2) unit cells with even 
number of unit cells in T) direction and odd in T 2 direc¬ 
tion. It is straight forward to see that on those samples, 
one can never avoid domain walls. In the following , we 
will show the different lattice quantum numbers are ac¬ 
tually caused by different quantum fluctuations along the 
domain wall. 

Let us focus on the simplest case where the sample is 
a chain along the Ti direction with even number of unit 
cells, but with only one unit cell along the T 2 direction, 
as shown in Fig.(16). Notice that the periodic boundary 
condition is imposed. Then, every site is connected with 
four bonds. Arrows on Fig.(16j,,b) denotes the direction 
of singlet bonds of Class-I and Class-II respectively. The 
direction of single bonds (red ones) on Fig.( 16 :-f) follows 


FIG. 16. (a),(b): ID chain in Ti direction. The periodic 

boundary condition is imposed, so every site shares four 
bonds. If two sites connected by a bond form a spin sin¬ 
glet, then we require the direction of the singlet is along the 
arrow on the bond. Arrows on (a) is consistent with Class- 
1, while arrows on (b) is consistent with Class-II. (c),(d),(e); 
Four possible domain wall configurations connected by sym¬ 
metry T 2 and Rtt- The direction of spin singlets (red thick 
bonds) follows arrows in (a). 


the convention in Fig.(|I^. We label these four domain 
wall states as \4>i), 1^2)7103) |(/) 4 ) respectively. In the 

following, we will show how these four simple VBS con¬ 
figurations give states with different quantum numbers. 

On the Hilbert space spanned by these four basis, the 
representation of Ti , T 2 and reads 



1 

0 


/I 

0 

0 


1 

0 

0 


0 

1 

0 

0 

0 

0 

0 

1 ,72 = 

0 

0 

1 

0 

Vo 

0 

0 

l ) 

Vo 

0 

0 

ij 


/O 0 1 0\ 
[ 0 0 0 1 I 
10 0 0 
\0 1 0 0 / 
(El) 


Then, eigenstates and eigenvalues for these three matrix 
are as following 


Eigenstates 

Ti 

72 

Rtt 

|<7’i) + \h) + \h) + l</>4) 

1 

1 

1 

ki) + ^ 2 ) - \h) - 04) 

1 

1 

-1 

0l) - 02) + 03) - 04) 

-1 

1 

1 

0l) - 02) - 03) + 04) 

-1 

1 

-1 


So different superpositions of VBS configurations give 
different quantum numbers. The above picture is sim¬ 
ilar to spin liquid phases (RVB states). Different spin 
liquid phases are distinguished by relative phases of dif¬ 
ferent configurations of bond coverings. These relative 
phase factors may result in different quantum numbers 
on some finite size sample. However, unlike spin liquid 
phases, fluctuation along the VBS domain wall is essen¬ 
tially ID physics. This remains true after considering 
samples with more unit cells along the T 2 direction. In 
one dimension, in the thermodynamic limit, the system 
will be pinned to a particular VBS configuration, and 
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a 



-c 


b 


FIG. 17. A site on the square lattice with four virtual legs 
labeled as a, b, c, d. 


the information of phase factors are lost. So we believe 
different ways of fluctuations along the VBS domain wall 
will not give new phases. 

Now, let us check the correctness of the fluctuation- 
along-domain-wall picture by studying the possible su¬ 
perposition of these four states in Class-I and Class-II 
respectively. By comparing directions of singlet bonds, 
We conclude that for Class-I, the possible superposition 
of these four states is -|- \(j)2) + -I- \4'4))- Fur¬ 
ther, we observe |^i) — 1^2) — — \4'4) is also consis¬ 

tent with Class-I with a non-contractible flux loops in T2 
direction. While for Class-II, there are two cases. For 
chains with (4n + 2) uc, the possible superposition of 
these four states is (—|0i) — |<^2) + \(l>3) + \4'4)) as well as 
\4’i) ~ \4’2) + I'/'a) — \4’4) (with no-ncontractible flux loop 
in T2 direction). For chains with 4n uc, we get the same 
result as Class-I. Thus, for systems with (4n -|- 2) sites, 
the two classes always have different quantum numbers. 
The above observation is consistent with the quantum 
numbers of Class-I and Class-II obtained in the previous 
appendix. 


Appendix F: An example on the square lattice 

As a pedagogical example, here we present the classi¬ 
fication of symmetric PEPS with IGG = Z2 for systems 
on the square lattice with a half-integer spin per site, in 
the presence of lattice translation, lattice C4 rotation and 
spin SU{2) rotation symmetries. 

The lattice symmetry group is generated by Ti, T2, (74, 
which transform the virtual leg labeled by (x, y, i) as: 

Ti{x,y,i) = (a; -b l,y,i) 

T2ix,y,i) = {x,y + l,i) 

G4{x, y, i) = {-y, x, G^ii)), (FI) 

where i = a,b,c,d labels the four virtual legs on a site 
tensor at (x, y) in a counter-clockwise fashion, as shown 
in Fig. (Hii- The counter-clockwise (74 rotates the legs as: 

(74(a) = 6, (74(6) = c, G4{c)=d, (74(d) = a. (F2) 

These generators satisfy the following identities which 


define the space group: 

[1]: 

[ 2 ]: 

[3] : 

[4] : 


TiT2T^^T^^ = e 
G4 ^TiG 4T2 = e 
C4^T2C4T^^ = e 
Gi=e 


And the on-site physical spin rotation by an angle 0 
around the spin axis n: Ugfi^ which forms a half-integer 
spin irrep of the SU{2), commutes with all lattice sym¬ 
metries: 


g-^ -Uen-g-U^n V5 = Ti,r2,(74 (F3) 


Below we solve the implementations of these symmetries 
on PEPS with IGG = Z2 = {I, J}, step by step. As 
discussed in Sec jlll A] each symmetry element R is asso¬ 
ciated with its own ry-ambiguity and e-ambiguity. And 
the PEPS itself has an ^-ambiguity and a E-ambiguity, 
a statement unrelated to specific symmetry elements. 

( 1 ) Choose the virtual basis such that i{x,yG) = J = 
Idi 0 (—lua) is a site and leg independent diagonal ma¬ 
trix, according to the discussion in |IIIB This determines 
J up to an overall ±I sign, which we will use in step-(IO). 
(Note that we have not attached any physical meanings 
for the Di and D2 sectors yet.) All the remaining V- 
ambiguity matrices (could be site and leg dependent) as 
well as all the Wr matrices (could be site and leg depen¬ 
dent) commute with J, so they are block diagonal and 
act within the Di and D2 subspaces. 

(2) Consider identity-[2]. Applying Eq.([^ to this 
identity: 


2;, (74(i))IFTi {-y, X, G4{i))WcA-y - 1, X, Gi^i)) 
■ {x,y + l,i)= il[2]X[2] {x, y, i) (F4) 


Since Wti appears in this equation only once, we can 
always use the rjTi ambiguity: Wt^ JFFti to tune 
r][ 2 ] = F This fixes the relative 77 -ambiguity for Wti,Wt 2 , 
and still leaves an overall 77 -ambiguity Wt^ , Wtj —>■ 

(3) Using the remaining E-ambiguity to transform 
Wt 2 , according to Eq. ( [^ : 

WT2{x,y,i) V{x,y,i)WT2{x,y,i)V~^ix,y - l,i)- 

(F5) 


So we can set Wt2 (x, y, i) = I. There is no cts ambiguity 
left. The remaining E-ambiguity satisfies V{x,y,i) = 
V{x,0,i). 

Next, we use this remaining E-ambiguity to transform 
Wti along the row of sites at y = 0: 

WTi(a:,0,7) ^ E(a:,0,7)IETi(a:,0,z)E-^(a:- 1,0,7), 

(F6) 


so that IEti(x, 0 ,z) = I. Now the remaining 
E-ambiguity is site independent but could be leg- 
dependent: V{x,y,i) = y^)- The remaining am¬ 
biguity satisfy cti (a:, 0, i) = I. 
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Consider identity-[1]. 
tity: 


Applying Eq.(28) to this iden- 


WTA^,V,i)WT^{x- - l,i) 

■W;^^{x,y,i) =yi2Xi2{x,y,i) (F7) 


The remaining ry-ambiguity for Wt^ , cannot tune 
away 7712 . Using results obtained previously in the step, 
we have 


ambiguities left. There is a remaining overall ry-ambiguity 
for Wti , Wt 2 ■ 

(4) Consider identity-[2] again. Based on the discus¬ 
sion in step-(2), we have: 

2 :, C4{i})Wt, i-y, X, C4{{})WcA-y - 1 . x, C4{i)) 
■Wr^ix.y+ l,i) = Xl 2 ]ix,y,i). (F14) 

Plugging in the results in step-(3), we obtain: 


Wt, (x, y, i)Wj,^^{x, y-l,i) = m 2 Xi 2 {x, y, i) (F8) 

We then use the remaining ambiguity to transform 
X 12 , according to Fq.(36): 

Xi 2 {.x,y,i) -)> Xi 2 {x, y, i)eTdx, y, i)eTi{x, y - l,i). (F9) 


So we can set Xi 2 ix, y,i) = 1. After this Fq.(F8) leads 
to WT 2 {x,y,i) = r]\ 2 ^Ti{x, 0 ,i) = 77^21 and there is no 
cti ambiguity left. In this gauge all site tensors are the 
same: = T, while bond tensors are spatial 

dependent if the 7712 index is nontrivial. 

Next we study 0 Ti and ©Ta- First we use the $- 
ambiguity, which we have not used before, to transform 
©Ta- According to Fq.(34): 


0 Ta(a;,y) ^ QT2ix,y)<^{x,y)<^*{x,y - 1 ), 


(FIO) 


so we can set ©Ta(a^,2/) = !> and the remaining 4>- 
ambiguity satisfies ^{x,y) = <i)(a;,0). 

We then use the remaining ^-ambiguity to transform 
©Ti along the row of sites at 7 / = 0: 


©Ti(a;,0) ^ ©Ti(a:,0)$(a;,0)$*(a:- 1,0), (Fll) 


so we can set ©Ti{a:,0) = 1. The remaining 4>-ambiguity 
is only a site-independent overall phase: ^{x,y) = 
$(0,0) = $. Since we will not study time-reversal sym¬ 
metry here, it turns out that this overall phase ambiguity 
is not useful. 

Applying Fq.(29) to the identity-[1], we have 


0Ti(a:,7/)©r2(a:- 1, y)©^^ (a:, yy - 1)©^^ (x, yy) 

= tJ-i2Y[x*i2ix,y,i), (F12) 

i 


where yi 2 is the site-independent phase factor obtained 
when applying 7712 on a site tensor, as mentioned in |III B| 
Using the results obtained so far, we find QTi{x,y) — 
yif 20 Ti(a:,O) = 

We can summarize the results obtained in step-(3): 

Wt, {x, y, i) = 

WT2{x,y,i) = I 
^TAx,y) = Mi 2 

©Ta(a:,y) = l. (F13) 

The remaining $-ambiguity is an overall phase. The re¬ 
maining V(x, y,i) = V(i). There is no remaining cti , CTa 


^C4(-y^x,C4(i))Wc,(-y - l,x,C4(i)) = y42Xi2](x,y,i)- 

(F15) 

Now we use the ambiguity. Applying Eq.(|M|, this 
transforms X[ 2 ] as: 


Xi 2 ](x, y,i) -t 

X[ 2 ] (x, y, i)cci (-y - 1, X, C 4 {i))€c^ (-y, x, C4{i)). (F16) 


So we can set X[ 2 ](a:,y,*) = 1, and the remaining eq^- 
am biguit y satisfies cCi{x,y,i) = ec4(0, y,7). After this, 
Eq.( |F15| ) leads to Wcja;,y,7) = TyiflUcjO, y, i). 

(5) Consider identity-[3]. Applying Eq.(| 


^C4^(-y>2:,C4(z)) - l-WcS-y^x- l,C4{i))Wj.^{x,y,i) 

= iCiTXCiTix, y, i) (F17) 


Using results in step-(3) and step-(4), we have 

Wc^i0,x,C4{i))WcA0,x - 1,C'4(7)) =VCiTXCiT{x,y,i), 

(F18) 


so we know XCiT{x,y,i) = XCiTix,0,i), independent of 

y- 

we then can use the remaining -ambiguity from 
step-(4) to transform xCaT'- 


XCiT{x,Q,i) 

XCiT{.x,Q,i)ec^{Q,x,C4{i))ecA^^x- l,G4{i)). (F19) 

So we can set XC 4 T(a^, 0 , 7 ) = 1 and the remaining eq^- 
ambiguity is site independent: ec'^(x,y, 7 ) = cCii})- 
After this, the site-dependence of Wc^ is solved: 
1 ^ 04 ( 0 , y, 7 ) = 7 y^^ylUc^( 0 , 0, i). Together with results 
in step-(4): 


Wci{x,y,i) = ’nlWc^T^cA'i)^ 


(F20) 


where we defined lUc4(7) = lUc^(0,0,7). 

So far we have not used the ry-ambiguity for Wc^ ■ The 
remaining eci{x,y,i) = eC 4 ( 0 ) and there is still a re¬ 
maining V(x,y,i) = V(i) ambiguity. 

(6) Consider ©c^. Applying Eq.(29) to identity-[2], 
together with results in step-(2,4): 


0 C4(-2/>37)©Ti(-y,a:)©c^(-y- 1, x)©Ta(a:, y) = 1. 

(F21) 


Pluggin in results in step-(3), this leads to: 


QcAx,y) = yi|©C4(0>y)- 


(F22) 
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Similarly, we apply Eq.( 29 ) to identity-[ 3 ]: 

0C4(-y>a;)0C4(-J/,a^- 1 )Mi 2 = MC 4 T H ^C4T(a;: 2 /> 

(F 23 ) 

where we used results in step-( 3 ). Plugging in xc^ t = 1 , 
which has been obtained in step-( 5 ), and use Eq.(F 22 ), 
the site-dependence of 004 is solved: 

QcAx,y) = (F 24 ) 

where we introduced 0C4 = 0C4(O, 0). 

( 7 ) Consider identity-[ 4 ]. Applying Eq.( 28 ) and the 
site-dependence of obtained in step-( 5 ), we have: 


VEc4 {i)Wc, (C!(t))Wc, {CUt))Wc, (C'4(*)) 
= VCiXcAi), 


(F 25 ) 


where XCii'i) = XcAx,]),!,) since the above relation dic¬ 
tates xc'4 to be site independent. 

Applying Eq. ( 291 and the site-dependence of 004 ob¬ 
tained in step(6), we find: 

(F 26 ) 


But due to the definition of the x ~ g^oup, we know that 

XcM) = XcS^) and XcA^) =XcM)^ so 114X04(*) = 1 - 

Consequently 0O4 = ■ Naively there would be 

four allowed roots for a given choice of 7704- However, 
we have not used the ?7-ambiguity for VEo4- Under the 
ry-ambiguity transformation: TU04 —t J-W04, clearly 0O4 
transforms as: 0O4 —t MJ0C4- 

We will prove that = —1 in step-(lO); i.e., every site 
tensor is Z2-odd. Thus, we can use the 77-ambiguity for 
fUo4 to tune away the sign in 0O4, and there are only two 
independent values for the root. We made the following 
choice: 


0C4 = 1 or 7, if 77C4 = I; 

0C4 = or , if rja = J- (F 27 ) 


After this, there is no remaining 77-ambiguity for Wc^- 
Next, we can use the remaining V(i) ambiguity to 
transform 1^04(7): 

WcA^ ^ V{i)WcAi)V-\C^\i)). (F 28 ) 


So we can set Wc’4(5) = WcAA — WCi{d) = I. This 
leaves no remaining £04-ambiguity, because eCiA) = 
ep^(c) as required in the definition of the y — group. 
The remaining U-ambiguity is site and leg independent: 
V{x,y,i) = V. 

Coming back to Eq.(F 25 ), in the current gauge we 
have: 


WcAa) =^04X04^), yi = a,b,c,d. (F 29 ) 


One may worry that if 77C4 = J, we simply have 
Uc4(a) = XC4 and the ±1 sign here may be tuned away 
by redefining the J element (recall that there is such sign 
freedom as mentioned in step-(l)). But for the moment 
let us not use this sign freedom in the definition of J 
because it will be used later in step-( 10 ). Consequently 
after step-(lO), the XC4 index here cannot be tuned away. 

(8) Consider the on-site SU{ 2 ) symmetry. We can 
apply Eq.( | 28 [ ) for a group identity in the multiplication 
table of SU{ 2 )'- 


[ 6 * 2772 ] • [6»i71i] = [ 6 * 3713 ], (F30) 

and obtain: 

Let us focus on a single virtual leg (a;, 77,7), and consider 
all the possible [02672]j [^ifii]- Oi^e can then immediately 
see that both g[0.^nPi\e4n4\ and X[e2S2],[eini](a:,7/,7) must 
satisfy 2-cocycle conditions as a function of [82^2] and 
[0i77i] (see Appendixj^for detailed discussions). Namely, 
for a fixed (x, y, 7): 

0[e2n2],[eini] G H^{SU{2),Z2) = Zi, 
X[eA].[ 04 n 4 ]A,y,i) G HASU{2),Uil)) = Z, (F32) 


Because both 2 -cohomology groups are trivial, we 
find both Vi92n2],lSini] and Xl02fi2],l0ini]{x,y,i) are 2- 
coboundaries. Consequently, one can use the 77- 
ambiguities for Wgfi and the £@^(0:, 77,7)-ambiguities to 
set v\0-,n2 ].\0,nP = I and xie2n2],lSini] A^yA) = 1 - After 
this, Eq.(F 31 ) simply means that W0ft{x,y,i) forms a 


representation of 56 /( 2 ), \/(x,y,i)- 


^03U3i^^yA)W02fi2{x,y,iW0^fiAx,y,i) = F (F 33 ) 


( 9 ) Study the site and leg dependence of Wes(a;, 77,7). 
For any space group symmetry element g, we have the 
group identity: 


g ^ ■ [ 077 ] • g ■ [ 077 ] ^ = e. (F 34 ) 

Applying Eq.(|^ to this identity, we find: 

^A^igA, y, i))W0fi{g{x, y, i))Wg{g{x, y, i))W^}{x, y, i) 
= Vg,lOfi]Xg,[0n] A, y, i) (F 35 ) 


Because we have already determined the form of 
Wg(a;,77,7) in step-( 3 , 5 , 7 ), and Wg(a;,77,7) can only be 
a power of J up to a factor, we conclude that Wg{x, y, i) 
commutes with W0fi{x,y,i)- So the above equation re¬ 
duces to: 


W 0 fi{g{x,y,i))Wg.}{x,y,i) = Vg,[0n]Xg,[0n]{x,y,i)- 

(F 36 ) 


Consequently XC4(*) = XC4 is leg independent, and 
XC4 = ±1 since e.g., xC4(a) = Xc4(c)- 


Next, one can plug in [02672] • [ 0if7i] = [03673] on the RHS 
of this equation and apply Eq.(F 33 ). One then concludes 
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that, for any fixed {x,y,i), 77g,[es] iXg,[Sn]{x,y,i)) must 
be a representation of SU{2) in Z 2 (U(i)). But such 
representations must be trivial. Th erefor e, rjg^iefi] = I 
and Xg,[en]{x,y,i) = 1, '^{x,y,i). Eq.(|F36 1 dictates that 
Wgfi must be site and leg independent: 


Wen{x,y,i) = Wefi, V(a:,y,^). (F37) 

In addition, we can consider Qgfi(x,y). Due to the 
definition of 0, the following is true for any site tensor 


=QeuWenUeuoT^^'y\ 


(F38) 


However, since is the SU{2) representation on the 
physical leg and we already showed that Wgs also form 
a representation of SU(2), Qg^ must be a representa¬ 
tion of SU{2) in U(l), which again must be trivial. So 
Qefi{x,y) = 1. 

( 10 ) Finally, consider the 9 = 2tt SU{2) rotation. 
Fq.(F38l for this particular case becomes: 


rp{^,y) ^ - Wg=2^oT^^’V\ 


(F39) 


where we use the fact that the physical spin is half-integer 
and 05=271- = 1 obtained in step-(9). This means that 
Wg= 2 Tn a site and leg independent transformation, must 
be an element in IGG and thus can be written as: 


Wg=27T = 'ne=27TXe=2nix,yG), 


(F40) 


we see that xe= 2 Tr{x,y,i) = Xe= 2 -K is site and leg inde¬ 
pendent. Due to the definition of the x ~ group, this 
limi ts X 9 = 2 Tr = ±1- Plugging the above equation back in 
Eq.(F39l, one has: 


—ye= 2 TvY\_Xe= 2 Tv = —ye= 2 Tv = ^ (F41) 

i 

So /i 5 = 27 r = — 1- This dictates that 775 = 2 ^ = J and every 
site tensor must be Z 2 odd: pj = — 1 . 

We then have: 

We= 2 . = xe=2n3 (F42) 

Note that we still have a sign ambiguity in the definition 
of J, as mentioned in step-(l). We now use this sign 
ambiguity to set xe= 2 TT = 1, and H 5 = 27 r = J. 

( 11 ) We still have the remaining V{x,y,i) = V am¬ 
biguity. Now we use this ambiguity to tranform the site 
and leg independent Wga to the standard form: 

W"5n = ©ti(In. (F43) 


namely each virtual leg is a direct sum of rii number of Si 
with different spin representation Si. Here we are only 
left with the H-ambiguity that is a direct sum of the 
similarity transformations acting in the I„. spaces. 

Summary: We find in the presence of translational 
symmetry, G 4 symmetry, spin-rotational symmetry, the 
IGG = Z 2 symmetric PEPS for a half-integer spin sys¬ 
tem on the square lattice are classified by the following 
three sets of algebraic data: 

1- 112 ,'>lC 4 T,riCi £ IGG = {I, J}. 

2- XCi which can be ±1. 


3. 0Ci which can choose values as defined in 
Eq.(lF27l). 


Since each index can choose two values, there are 2® = 32 
classes. These indices completely determine the transfor¬ 
mation rules of the site and bond tensors as: 


Wgn[x,y,i) = Wgn = ©,^i(In. © 

J = Wg^2.. 


WT 2 {x,y,i) = I, 

Wc4{x,y,a) = XCiVlWciT'nCi 
WcAx,y,blc/d) = 

(F44) 


and: 


Qgn{x,y) = 1 
0Ti(a;,2/) = M 12 . 

QT2{x,y) = 1 , 

^cAx,y) = (F45) 

Here /ii 2 = 1(—1) if V 12 — I (J), and similarly pc/r = 
1(-1) if yCiT = I (J)- 0C4 = 1 or 7(0 c 4 = 6”/"^ or 
e—/4) =i (J). 

After the physical half-integer spin is specified, e.g. 
S' = 1/2 or S' = 3/2, we know the transformation rules for 
both physical and virtual legs. One can thus determine 
the generic form of the symmetric tensor network for each 
class and use it for numerical simulation as discussed in 
Sec IlHC] 
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